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ABSTRACT 


The  objective  of  the  reported  research  was  to  investigate  the 
relationship  of  radiation  heat  transfer  to  the  in-cylinder  species  distribution  in 
direct  injection- type  diesel  engine  combustion. 

In  order  to  implement  this  objective,  both  the  experimental  and 
computational  methods  were  employed  during  the  course  of  the  study.  In  the 
experiment,  in  order  to  obtain  the  in-cylinder  species  distribution  at  successive  en¬ 
gine  crank  angles,  a  rapid  sampling  valve  has  been  used  on  a  laboratory- built  single 
cylinder  direct  injection-type  diesel  engine.  Among  the  unique  features  of  the  new 
sampling  system  is  dilution  of  the  sampled  gas  immediately  after  sampling  to  mini¬ 
mize  physical  and  chemical  reactions  which  might  take  place  in  the  sample  train. 

The  sampling  system  was  tested  in  a  bench  apparatus  before  being  directly  used  in 
the  actual  engine.  At  the  same  time,  radiative  heat  transfer  through  the  engine 
combustion  chamber  wall  was  measured  by  using  thin-film  thermocouple  with 
quartz  window  in  front.  This  was  to  measure  portions  of  radiative  and  convective 
heat  transfer,  respectively.  The  engine  for  the  data  acquisition  was  constructed 
basically  by  combining  a  reciprocating  system  of  an  ASTM  gasoline  engine  and  a 
newly  constructed  engine  cylinder  head.  Although,  the  development  of  the  engine 
has  offered  an  wide  range  of  educational  experiences,  while  some  in-cylinder  sam¬ 
pling  and  radiation  heat  transfer  data  were  obtained,  we  encountered  a  series  of 
problems  in  the  engine  requiring  replacement  of  engine  components,  which  were  of¬ 
ten  difficult  to  obtain.  At  the  end,  the  engine  was  regarded  unable  to  meet  the  need 
of  operating  for  an  extended  period  of  time  for  data  gathering.  Consequently,  it 
was  decided  to  use  a  new  single  cylinder  engine  that  was  made  available  through 
DOD-URIP.  While  the  engine  is  being  prepared  for  continuation  of  the  experiment, 
more  efforts  were  concentrated  on  the  computational  activities. 

In  the  computational  work,  a  radiation  heat  transfer  model  was 
developed  to  predict  spectral  radiation  heat  transfer  through  the  engine  combustion 
chamber  wall  for  given  in-cylinder  species  distributions.  Several  unique  approaches 
were  taken  in  the  modeling.  For  reporting  the  resuits  from  this  work  eight  separate 
manuscripts  have  been  prepared  under  the  present  sponsorship:  seven  for  publica¬ 
tion  in  journals  and  one  in  an  international  symposium.  They  are  for  reporting  our 
discovery  of  a  new  integral  function,  a  new  exact  solution  of  biackbody  functions,  a 
new  empirical  equation  of  adiabatic  flame  temperature  of  fuel  air  systems,  an 
analytical  solution  of  equation  of  radiation  functions,  a  parametric  analysis  of 
diesel  radiation  heat  transfer,  etc. 
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STATEMENT  OF  PROBLEM 


The  proposed  study  was  to  conduct  research  to  obtain  a  better 
understanding  of  the  processes  of  particulate  formation  and  heat  transfer  in  direct 
injection- type  diesel  engines.  Specifically,  it  was  to  obtain  time-  and  space- resolved 
concentration  of  gaseous  and  particulate,  to  measure  components  of  instantaneous 
heat  transfer  through  convective  and  radiative  processes,  and  to  find  the  relation¬ 
ship  of  the  in-cylinder  species  concentration  distribution  to  radiative  heat  transfer 
data  obtained  in  the  combustion  chamber. 

The  research  objective  was  to  be  implemented  by  carrying  out 
the  following  work  tasks:  (1)  Construction  of  a  new  single-cylinder  direct  injection- 
type  diesel  engine;  (2)  Fabrication  of  an  in-cylinder  sampling  valve  and  probes  to 
measure  heat  transfer  through  convective  and  radiative  processes;  (3)  Compilation 
of  in-cylinder  species  data  and  heat  transfer  measurements;  and  (4)  Development  of 
a  computational  model  to  understand  the  relationship  of  the  species  concentrations 
to  radiative  heat  transfer,  consequently  to  compare  results  from  both  the  computa¬ 
tional  and  experimental  approaches. 


The  research  was  to  be  conducted  at  the  Department  of 
Mechanical  and  Aerospace  Engineering,  Rutgers,  the  State  University  of  New 
Jersey,  by  supporting  two  graduate  students  under  the  supervision  of  Dr.  K.  T.  Rhee 
One  student  was  to  be  assigned  to  experimental  part  of  the  study,  i.e.,  to  set  up  an 
engine  apparatus,  measure  in-cylinder  species  and  radiative  heat  transfer  from  the 
engine  combustion  chamber.  The  other  student  was  to  develop  a  computational 
model  of  radiation  heat  transfer  in  diesel  combustion  to  predict  radiation  heat 
transfer  through  the  wall  when  the  combustion  product  distribution  is  given. 


1. 


Although  the  current  excess  oil  supply  creating  ecstatic  pump 
prices  of  automobile  fuels  seems  to  have  us  forget  the  painful  experiences  under¬ 
gone  during  the  oil-shortage  on  the  early  70s,  no  one  expected  such  low  pricing  to 
occur  when  the  present  research  was  proposed,  The  global  oil  shortage  has  not  gone 
away  from  us  at  all,  but  is  a  temporary  disruption  in  the  supply-demand  structure. 
Actually,  domestic  petroleum  resources  are  rapidly  being  depleted.  Evidently,  im¬ 
port  of  foreign  oil  will  not  be  the  solution  for  filling  the  increasing  domestic-oil  gap 
due  to  the  rapid  bottoming  our  oil- balance.  Failure  to  implement  solutions  to  the 
problem  in  a  timely  way  may  lead  to  serious  economic,  social  and  political 
disruption. 

One  viable  approach  to  the  oil-shortage  is  the  efficient  use  of 
limited  energy  resources.  As  the  transportation  system,  especially  automobile,  is  the 
largest  consumer  of  liquid  fuels,  advanced  automobile  technology  is  needed  badly  to 
help  find  the  solution  of  the  domestic  energy  shortage.  Note  that  the  severity  as¬ 
sociated  with  the  shortage  may  be  even  stronger  than  before  if  it  comes  to  reality 
again  if  not  prepared  for.  For  the  future,  improved  technology  for  manufacture 
and  operation  should  be  identified  to  achieve  better  fuel  economy  or  to  use  alterna¬ 
tive  fuels. 

In  spite  of  the  current  pause  in  development  of  synthetic  fuels, 
sooner  or  later  the  United  States  will  be  making  synthetic  liquid  fuels  from  coal  and 
shale.  One  of  the  first  problems  to  be  solved  is  the  determination  of  the  composi¬ 
tion  of  the  new  fuels.  It  is  particularly  important  to  know  the  compositions  of  the 
new  fuels  because  existing  highway  transportation  systems  are  highly  sensitive  to 
fuel  grade.  It  is  likely  that  new  fuel  compositions  will  not  be  much  different  from 
those  of  existing  fuels.  On  the  other  hand,  there  is  a  consensus  that  the 
carbon/hydrogen  ratio  of  the  new  fuels  may  be  somewhat  higher.  The  former  con¬ 
sideration  is  based  on  the  fact  that  there  are  millions  of  engines  already  in  the  field, 
and  their  owners  expect  many  years  of  satisfactory  service  even  by  using  new  fuels 
in  the  future.  The  opinion  on  the  C/H  ratio  is  expressed  because  the  price  of  new 
fuel  will  be  very  much  dependent  upon  the  cost  of  the  hydrogen  to  be  used  in  the 
process  of  new  fuel  production. 

The  diesel  combustion  engine  produces  a  high  power  output 
per  gallon  of  fuel,  resulting  in  increased  fuel  economy  in  transportation.  In  spite  of 
the  advantage  of  high  fuel  economy,  there  are  several  main  problems  to  be  solved 
for  becoming  a  better  propulsion  system.  Among  the  problems  associated  with  diesel 
engine  powered  vehicles  is  the  emission  of  diesel  smoke  (or  particulate).  Such  smoke 
emission  is  expected  to  be  even  worse  if  new  fuels  with  high  C/H  ratios  are  used. 

The  diesel  smoke  not  only  looks  bad  but  also  contains  polynuclear  aromatic 
hydrocarbons  (PAH),  some  of  which  are  known  to  cause  mutagenic  actions  on  living 
cells.  In  an  attempt  to  reduce  the  smoke  emission  from  diesel  engines,  the  U.S.  EP.A 
has  established  an  exhaust  particulate  regulation  to  apply  to  diesel-powered  vehicles. 
Some  studies  (Lestz  et  al.,  1982;  Bradow,  et  al.,  1982)  suggested  that  pyrene,  one  of 
PAHs,  might  react  with  NOx  in  diesel  engines  to  form  nitropvrene  which  is  con¬ 
sidered  to  be  one  of  the  most  active  mutagenic  reagents.  It  is  highly  desirable  to 
know  where,  when  and  how  this  compound  is  formed  in  diesel  engines.  Ultimately, 
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one  wants  to  minimize  such  harmful  reagents  contained  in  diesel  exhaust.  In 
addition,  diesel  exhaust  produces  unpleasant  smells,  which  seem  to  be  closely  related 
to  the  composition  and  amount  of  smoke  in  the  emission. 


The  direct  injection -type  diesel  engine  provides  better  fuel 
economy  than  the  indirect  injection- type  diesel  engine.  This  is  due  mainly  to  the 
smaller  ratio  of  combustion  chamber  surface  ratio  to  its  volume,  so  that  the  heat  loss 
to  the  chamber  wall  is  less.  And  it  is  well  known  that  the  great  portion  of  overall 
heat  transfer  in  diesel  combustion  occurs  through  radiative  heat  transfer.  As  the 
radiative  heat  transfer  is  affected  by  the  content  of  carbon  particles  appearing  in 
the  combustion  flames,  the  heat  loss  may  even  be  increased  as  the  high  C/H  ratio 
fuels  are  used  in  diesel  engines.  The  extent  of  smoke  appearing  in  the  exhaust 
during  combustion  is  considered  to  relate  closely  to  the  heat  loss  in  the  combustion 
chamber.  Without  mentioning  other  complex  processes  affecting  efficiency,  the 
large  amount  of  smoke  appearing  in  the  exhaust,  therefore,  may  represent  a  red 
signal  in  diesel  economy. 

At  the  time  the  present  research  proposal  was  submitted,  there 
had  already  been  vigorous  interest  in  eliminating  the  cooling  system  of  the  diesel 
engine,  the  adiabatic  (or  low-cooled)  engine  concept.  Although  the  waste  energy 
recovery  by  the  approach  may  not  be  significant,  the  problems  associated  with  the 
elimination  of  the  cooling  system  of  the  engine  remain  as  difficult  tasks  to  be 
solved  in  the  new  engine  development.  Since  the  conventional  diesel  engine  rejects 
as  much  as  35%  of  the  total  energy  produced  by  the  combustion  processes  to  the 
cooling  and  lubrication  systems,  the  elimination  of  the  cooling  devices  poses  con¬ 
siderable  difficulty  in  various  aspects  of  engine  production:  design  modification, 
material,  tribology,  fuel  system,  etc.  Specifically,  greater  concerns  are  those  as¬ 
sociated  with  thermal  loading,  processing/development  of  insulating  materials  and 
finding  methods  of  long-term  bonding  of  such  materials  to  the  combustion  chamber, 
and  lubrication/wear  reduction  between  the  cylinder  liner  and  piston  rings.  Since 
those  issues  are  closely  dictated  by  the  greater  thermal  loading  than  the  engine  com¬ 
ponents  can  withstand,  it  is  highly  desirable  to  have  a  better  understanding  of  heat 
transfer  in  the  engine.  This  is  not  only  for  the  fundamental  issue  of  better  under¬ 
standing  the  processes  but  also  for  helping  to  develop  better  low-cooled  engines. 

For  example,  when  the  engine  performance  and  endurance  to  survive  of  the  newly 
designed  engine  concept  (e.g.,  turbocompound)  is  analyzed,  one  of  the  first  missing 
information  for  carrying  out  the  task  is  suitable  models  predicting  heat  transfer  in 
the  engine. 


In  an  attempt  to  obtain  a  better  understanding  of  diesel  com¬ 
bustion  processes,  Rutgers,  University  conducted  diesel  engine  combustion  research 
for  some  years  before  the  award  of  the  present  research  contract.  The  objective  of 
the  previous  work  was  consistent  with  the  proposed  research  goals  in  the  present 
study  enabling  efficient  continuation  of  the  engine  study  in  our  engine  laboratory. 
The  construction  of  engine  apparatus  for  the  proposed  study  was  in  good  progress 
and  a  preliminary  radiative  heat  transfer  model  was  almost  complete  when  the 
proposal  was  submitted. 

It  was  expected  to  obtain  a  better  understanding  of  the 
relationship  of  radiation  heat  transfer  to  in-cylinder  species  distribution  in  diesel 
combustion  upon  the  completion  of  the  proposed  research  work. 
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The  objectives  of  the  proposed  work  under  the  sponsorship  were 


(1)  to  obtain  time-  and  space-resoived  particulate  data  (and  CO2,  CO,  0-<, 
local  fuel/air  ratio)  in  a  direct  injection-type  diesel  engine  combustion 
chamber, 

(2)  to  measure  the  instantaneous  heat  transfer  (convective  and  radiative, 
respectively)  through  the  combustion  chamber  wail, 

(3)  to  construct  a  computational  model  of  radiation  heat  transfer  to  predict 
radiative  heat  transfer  through  the  wall  for  given  in-cylinder  species 
distributions  in  diesel  combustion,  and 

(4)  to  investigate  the  formation  processes  of  nitropyrene  in  diesel  combustion 
by  determining  its  concentration  during  the  reaction  period. 


3.  RESEARCH  METHODOLOGY 

The  following  explains  work  tasks  that  were  to  be  carried  out 
for  implementing  the  objectives. 

(1)  A  new  single-cylinder  direct  injection-type  diesel  engine  which  were 
being  built  under  the  previous  research  contract  was  to  be  completed  for 
the  study. 

(2)  A  new  sampling  system  having  dilution  for  minimizing  physical  and 
chemical  reactions  that  might  take  in  the  sample  train  was  to  be  developed 
for  determining  in-cylinder  species  distributions. 

(3)  Heat  transfer  probes  were  to  be  fabricated  for  the  measurement  of 
instantaneous  heat  transfer  through  the  combustion  chamber  wall. 

(4)  The  preliminary  2-D  radiation  heat  transfer  model  was  to  be  improved 
for  constructing  an  extensive  3-D  non-gray  radiation  heat  transfer  model. 

Upon  the  construction  of  the  engine  set-up  and  engine  diagnos¬ 
tic  tools,  in-cylinder  gaseous  species  and  particulate  concentrations  were  to  be  taken 
at  different  engine  crank  angles  and  different  locations  in  the  combustion  chamber. 

At  the  same  time,  radiation  heat  transfer  data  were  to  be  obtained  at  corresponding 
engine  conditions. 

Once  the  species  distribution  and  radiation  heat  transfer  data 
in  the  combustion  chamber  were  obtained,  predicted  radiation  heat  transfer  using 
the  new  computational  radiation  model  based  on  the  in-cylinder  species  distribution 
was  to  be  compared  with  measured  radiation  heat  transfer  data. 

The  study  of  pyrene  in  diesel  combustion  was  to  be  conducted 
in  collaboration  with  a  properly  instrumented  group  for  the  analysis  of  samples, 
such  as,  EPA  Mobile  Source  Emission  Branch:  Time-resolved  in-cylinder  samples 
taken  from  the  engine  using  the  new  sampling  system  were  sent  to  the  analytical 
group  for  finding  the  history  of  pyrene  formation  in  diesel  combustion.  This  ap¬ 
proach  was  expected  to  provide  information  as  to  when  and  where  pyrene  is  formed 
in  a  diesel  engine  combustion  chamber. 
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The  cumulative  results  achieved  in  the  study  are  explained  in 
two  main  parts,  i.e.,  experimental  and  computational  approaches.  Included  in  the 
experimental  side  are  the  construction  of  engine  apparatus,  in-cylinder  probes  and  a 
sample  processing  system.  Extensive  results  from  the  computational  work  are 
reported. 


4-1.  EXPERIMENT 

A.  CONSTRUCTION  OF  APPARATUS 

ENGINE  SET-UP.  The  construction  of  the  new  engine  was  made  by  combin¬ 
ing  a  reciprocating  unit  from  an  ASTM  gasoline  engine,  a  new  cylinder  head,  a  new 
cylinder  block,  and  as  many  stock  parts  as  possible  to  minimize  the  machining  cost. 
For  the  construction  of  the  system  using  various  engine  components  from  different 
sources,  extensive  modifications  were  needed  in  those  components.  Among  the 
modifications  made  in  the  construction  are  balancing  of  the  reciprocating  system, 
machining  of  the  cylinder  block,  special  design  of  injection  nozzle  for  the  new 
engine,  etc.  Note  the  design  details  of  individual  engine  components  are  explained 
in  a  separate  report  (DOT/RSPA/DMA-50/84-4).  Although  some  of  its  content  may 
seem  somewhat  repetitive  because  of  various  subsequent  improvement  and  redesign 
made  in  the  system  during  the  present  research  contract,  they  are  briefly  explained 
in  Appendix- 1. 

IN-CYLINDER  SAMPLING  PROBE.  The  probe  for  obtaining  in-cylinder 
species  determination  was  newly  developed  and  tested  using  our  laboratory  built 
engine.  Unlike  the  conventional  intermittant  flow-type  sampling  valve,  the  basic 
design  idea  of  the  new  sampling  system  is  based  on  a  continuous  sampling  valve 
which  was  developed  by  Rhee,  et  al.  (1978)  for  minimizing  the  potential  effect  of 
quenched  boundary  layer  formed  around  the  probe  hole  on  the  measurements.  The 
new  system  constructed  under  the  present  support  was  to  include  an  in-probe  dilu¬ 
tion  method  in  order  to  eliminate  or  minimize  the  physical  and  chemical  reaction 
which  might  take  place  in  the  sample  train.  The  disturbance  of  the  flow  field  by 
the  presence  of  the  probe  in  the  cylinder  may  be  greater  with  the  present  continuous 
flo-.v-type  system  than  the  intermittent  flow-type  as  some  limited  bench  test 
exhibited.  However,  the  disturbance,  if  noticeable  at  all,  is  not  considered  to  sig¬ 
nificantly  affect  the  sample  because  the  volume  of  the  sample  being  taken  by  the 
probe  is  from  the  uninterrupted  part  of  the  swirling  in-cylinder  charge:  at  most,  the 
disturbance  will  affect  the  composition  of  the  combustion  products  downstream  of 
the  probe.  The  details  of  the  probe  are  explained  in  Appendix-II. 

HEAT  TRANSFER  PROBES.  In  spite  of  strong  interest  by  many  auto  en¬ 
gineers  in  measurement  of  heat  transfer  through  the  combustion  chamber  wall,  the 
measurement  in  diesel  engines  has  not  been  properly  made  due  to  the  inherent  na¬ 
ture  of  diesel  combustion,  namely,  high  soot  formations.  The  high  soot  content  in 
the  cylinder  combustion  products  rapidly  fouls  the  windows  installed  for  optical 
access  and  the  thermocouple  probe  placed  on  the  wall.  Even  neglecting  the  problems 
associated  with  the  soot  fouling  formation  on  the  probes,  the  separation  of 


the  portions  of  heat  transfer  by  convective  and  radiative  processes  remains  as  a  dif¬ 
ficult  task  to  be  solved  in  the  diesel  heat  transfer  analysis.  Discussing  the 
methodology  for  the  goal,  there  are  two  basic  approaches  that  one  could  choose  for 
separating  the  individual  effects:  one  is  to  isolate  the  portion  of  convection  heat 
transfer  in  a  measurement,  which  is  basically  a  measurement  of  radiation  heat 
transfer;  the  other  is  to  exclude  the  radiation  effect  in  a  measurement,  which  is, 
then,  a  measurement  of  convective  heat  transfer.  Although  both  approaches  have 
carefully  been  studied  in  the  present  research,  since  the  former  method  has  been 
studies  more  intensively  its  details  are  only  discussed  here.  (Presently,  the  latter 
method  is  at  the  evaluation  stage  and  its  discussion  will  be  made  upon  conducting 
its  bench  test.)  The  surface  thermocouple  method  (Overby  et  al.,  1961)  has  been 
chosen  to  obtain  the  instantaneous  heat  transfer  through  the  combustion  chamber. 
Two  probes  were  placed  in  a  mounting  plug,  one  installed  flush  with  the  wall  for 
measurement  of  total  (convection  plus  radiation)  heat  transfer  and  the  other 
recessed  with  a  quartz  window  in  front  for  measurement  of  radiation  heat  transfer. 
More  details  of  the  probes  are  discussed  in  Appendix-Ill. 

SAMPLE  PROCESSING  SYSTEM.  In  order  to  comply  with  the  definition  of 
soot  by  EPA  in  our  measurement  of  in-cylinder  combustion  products,  several  sample 
processing  methods  are  chosen  in  the  present  work:  As  mentioned  previously,  dilu¬ 
tion  of  samples  in  the  sampling  probe  is  important  for  proper  quantification  of  the 
in-cylinder  soot  distribution.  Since  soot  is  defined  as  "any  material  (excluding  con¬ 
densed  water)  collected  by  filtering  a  sample  of  diluted  exhaust,  and  the  sample 
train  be  maintained  at  125°  F,"  the  sample  taken  out  of  our  probe  was  to  be  intro¬ 
duced  into  a  filter  unit  at  the  recommended  temperature.  The  measurement  of  soot 
was  to  be  made  by  weighting  a  filter  before  and  after  a  known  volume  of  sample 
was  passed  through  it.  A  new  sample  processing  unit  was  constructed  to  stabilize 
the  filter  (before  and  after  sampling)  where  the  environment  is  controlled  at  a 
designated  condition.  Because  of  the  relatively  unstable  nature  of  particulate 
samples,  it  was  necessary  to  use  a  special  caution  in  the  collection  and  processing  of 
data.  The  chamber  built  for  this  purpose  exhibites  a  very  satisfactory  result  by 
maintaining  a  typical  condition  of  relative  humidity  of  55%  and  a  constant  tem¬ 
perature  of  82°  F.  In  addition,  in  order  to  compare  the  in-cylinder  soot  formation 
to  the  exhaust  soot  emission,  a  new  exhaust  dilution  tunnel  was  constructed  in  our 
laboratory  by  duplicating  the  GMR  dilution  tunnel  (McDonald  et  al.,  1980). 


B.  ENGINE  DATA 

ENGINE  FAILURE.  While  we  were  conducting  some  preliminary  engine 
tests  by  using  the  heat  transfer  probe  and  the  sampling  valve,  there  were  a  series  of 
engine  problems  occurring  in  the  laboratory  built-engine:  some  were  simple  and 
others  were  seemingly  serious.  For  example,  we  had  to  find  three  new  sets  of  timing 
gear  and  oil  pump  driving  gear  (located  in  the  same  space  in  the  engine)  due  to 
repeated  chaotic  failure  of  the  components.  This  was  considered  to  occur  due  to 
the  greater  load  that  the  timing  gear  had  undergone  (by  driving  the  fuel  injection 
pump)  than  originally  designed  for  use  in  an  ASTM-CFR  gasoline  engine.  Another 
persisting  problem  encountered  in  the  set-up  was  the  control  of  leakage  in  the  newly 
constructed  camshaft  installed  for  driving  the  fuel  injection  pump.  During  the 
process  of  the  engine  construction  and  the  preliminary  measurements,  we  had  to 
have  a  major  overhaul  of  the  entire  engine  many  times.  The  engine  upon  its  con 
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struction  had  more  than  fifty  hours  of  compiled  operation  time  before  the  most 
serious  problem  occurred  lately.  The  problem  in  the  engine  found  in  the  recent 
overhaul  was  that  the  crankshaft  was  noticeably  deformed  and  needed  proper  repair 
or  replacement.  It  was  judged  that  the  engine  might  not  exhibit  satisfactory 
functioning  for  the  in-cylinder  data  acquisition  and  would  take  a  considerably  long 
period  of  operation.  In  the  meantime,  we  were  awarded  a  grant  under  the  DOD- 
URIP  to  develop  a  new  single  cylinder  direct  injection- type  diesel  engine  that  can 
be  used  as  both  a  conventional  water-cooled  mode  and  a  new  uncooled  operation. 
Although  the  development  of  the  new  engine  was  not  originally  proposed  for  replac¬ 
ing  our  laboratory-built  engine,  it  was  concluded  that  use  of  the  new  engine  would 
be  the  best  solution  for  future  data  acquisition.  Consequently,  any  systematic  data 
acquisition  has  not  yet  been  made  for  discussing  meaningful  engine  data  in  this 
report. 


Note  that,  in  addition  to  the  DOD-URIP  sponsored  engine, 
another  new  unit  of  Cummins  903-block  base  single  cylinder  direct  injection-type 
diesel  engine  was  made  available  to  us  (provided  by  Cummins)  for  fundamental 
studies  of  diesel  combustion.  This  additional  engine  has  been  lined  up  for  the 
measurement  of  instantaneous  temperature  variation  on  the  engine  combustion 
chamber  walls.  For  this,  a  new  grass-hopper  linkage  that  was  constructed  in  our 
laboratory  is  being  placed  in  the  engine. 


4-2.  COMPUTATIONAL  WORK 

While  various  difficult  problems  were  encountered  in  the  ex¬ 
perimental  work,  which  hampered  the  scheduled  data  acquisition,  some  very  sig¬ 
nificant  progress  has  been  made  in  the  computational  approach.  The  main  goal  in 
the  computational  work  was  to  develop  a  computer  model  to  predict  radiation  heat 
transfer  through  the  combustion  chamber  wall.  The  model  is  explained  in  the  fol¬ 
lowing  and  several  new  techniques  employed  in  the  modeling  are  separately 
explained. 


A.  RADIATION  HEAT  TRANSFER  MODEL 

In  the  construction  of  a  radiation  heat  transfer  model  for  com¬ 
bustion  systems,  several  serious  problems  are  faced,  mainly  due  to  the  directionality 
and  spectral  nature  of  radiation  heat  transmission.  Among  the  serious  problems  to 
be  solved  in  modeling,  as  discussed  next,  may  be 

(1)  A  viable  computation  of  spectral  volume  absorptance  of  a  mixture  of 
combustion  products  along  with  determination  of  emission  characteristics 
of  species, 

(2)  A  more  accurate  computational  method  for  solving  the  equation  of 
radiation  heat  transfer  along  individual  optical  paths  in  the 
hemispherical  volume  faced  by  locations  of  the  combustion  chamber  wall. 

(3)  A  detailed  description  of  in-reactor  distribution  of  combustion  products 
and  temperature  in  terms  of  a  suitable  coordinate  system  centered  at  an 
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individual  locations  where  heat  transfer  is  computed  and  its  coupling  with 
the  above  mentioned  equation  of  radiation  heat  transfer. 


SPECTRAL  VOLUME  ABSORPTANCE  OF  SOOT  GAS  MIXTURE.  When  a 
source  of  radiation  heat  transmission,  particularly  a  mixture  of  combustion  products 
including  soot,  is  evaluated,  it  is  necessary  to  find  the  volume  absorptance  and  tem¬ 
perature  of  the  mixture.  The  most  widely  accepted  method  of  computing  the 
spectral  volume  absorptance  of  such  a  mixture  seems  that  the  absorptance  by  soot 
cloud, and  that  by  gaseous  components,*^  are  found  to  find  the  summed 
absorptance,*^  .  In  the  computation  of  soot  radiation,  since  most  of  the  practical 
combustion  systems  produce  soot  having  much  smaller  thermal  wavelength,  the  scat¬ 
tering  effect  on  the  radiation  by  soot  cloud  is  minimal  and  its  computation  may  be 
made  by  using  the  Rayleigh  limit  expression.  The  expression  is  independent  of  the 
soot  size,  inversely  proportional  to  the  wavelength,  and  proportional  to  the  soot 
volume  fraction,  i.e.. 


C  n2 ( 1-k2 )  +  2; 


The  optical  constant,  n,  k,  the  real  and  imaginary  parts  of  the  refractive  indices, 
were  measured  from  the  intensity  of  polarizing  light  reflected  from  the  soot  par¬ 
ticles  and  the  measurements  were  fitted  to  the  dispersion  equation  (Dalzell  and 
Sarofim,  1969;  Lee  and  Tien,  1981).  More  details  of  the  computation  are  shown  in 
Appendix-IV,  where  results  from  our  preliminary  radiation  heat  transfer  analysis  in 
diesel  combustion  are  reported. 

For  computation  of  emission  from  gaseous  mixture,  the  empiri¬ 
cal  wide  band  model  proposed  by  Edwards  and  Balarkrishnan  (1973)  has  often  been 
used  by  others  in  the  past,  and  the  same  techniques  were  employed  in  the  present 
modeling.  The  details  of  this  method  for  our  modeling  is  extensively  explained  in 
Appendix-IV. 

The  remaining  question  for  the  computation  of  spectral  radia¬ 
tion  of  a  soot  laden  mixture  may  be  the  summation  of  both  the  soot  radiation  and 
gas  radiation.  In  view  of  a  previous  method  of  summation  of  gas  emittance.tg  and 
soot  emittance,«s  in  order  to  find  the  total  radiation,  t  for  gray  body  analysis,  i.e., 
t:  l-(l-«g)(i-«s)  (Siegel  and  Howell,  1981),  the  following  similar  approach  to  the 
method  was  employed  in  the  computation  of  spectral  volume  absorptance:  The 
spectral  volume  absorptances  for  gas,  x^g  and  soot  cloud,  *^s  were  computed  as 

g  =  C/X  (2-a) 

where,  C  is  found  from 

\  =  l-£xp( -Ci A )d\ 

and  in  Eq.  (1),  one  finds 
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The  summation,  then,  was  found  by  the  following: 

rx  =  (C+C')A  (3). 

The  validity  of  the  new  method  for  the  spectral  property,  Eq.  (3),  although  similar 
to  that  widely  used  for  gray  body  computations,  remains  to  be  studied. 


EQUATION  OF  RADIATION  HEAT  TRANSFER.  Although  the  equation  of 
radiation  heat  transfer  in  radiatively  participating  fluid  is  well  known,  its  solution 
for  practical  combustions  is  not  readily  solvable,  due  to  complexities  of  the 
equation.  Elaborating  further  on  the  abovementioned  problem,  Fig.  1  is  shown  for 
writing  the  governing  equation  of  the  process  as 
r 


I 


where, 


■  / 


vr 


e  ( r ' 
bX 


£xp(-/  «•  [ 

r  \ 


r" ) dr" )dr ' 


(4) 


i  C 

2xnc 


I  [ Exp( hc/XKT) -1 ] 

is  die  spectral  volume  absorptance  at  (r,0,f), 
x  10”^"  erg-sec,  and  k  *  1.380  x  lO'1”  erg/K.. 
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c  =  2.998  x  10JU  cm/sec,  h  =  6.625 
The  heat  flux  along  an  optical  path 
of  a  solid  angle  in  the  radiation  range  may  be  computed  from. 
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Among  the  main  difficulties  in  finding  a  solution  of  Eq.  (4)  are  that  is  a  function 
of  both  species  concentrations  and  temperature,  that  both  variables  are  a  f unction 
of  the  location  along  the  individual  optical  path,  r,  and  that  the  mathematical  com¬ 
plexity  of  the  equation  is  severe  so  as  to  employ  simplified  methods  (as  in  previous 
studies).  Since  the  species  distribution  along  the  individual  optical  paths  is  dis¬ 
cussed  in  the  following  section,  the  new  mathematical  approaches  taken  for  Eq.  (4) 
are  mentioned  here. 


The  equation  of  radiation  heat  transfer,  Eq.  (4),  due  to  its 
mathematical  complexity,  has  been  solved  by  using  various  simplified  methods  in 
the  past,  e.g.,  the  Monte-Carlo  method  (Steward  and  Cannon,  1971),  the  zonal  method 
(Hotel  and  Sarofim,  1967;  Menguc  et  al.,  1985)  and  the  use  of  geometric  factors 
(Chang  and  Rhee,  1983).  Since  we  reported  the  abovementioned  results  obtained 
from- using  the  geometric  factors  in  2-D  model,  a  new  analytical  method  of  non-gray 
3-D  radiation  heat  transfer  has  been  developed.  Equation  (4)  has  been  processed  by 
using  several  unique  analytical  methods  for  finding  analytical  solutions  for  both  op¬ 
tically  thin  and  thick  media  and  for  using  a  relatively  simple  numerical  method  for 
optically  medially  thin  medium.  The  details  of  the  solution  are  explained  in 
Appendix-VI. 

SPECIES  DISTRIBUTION  ALONG  INDIVIDUAL  OPTICAL  PATHS.  In  con¬ 
ducting  radiation  heat  transfer  computation  for  a  container  having  soot-laden  com¬ 
bustion  products,  one  of  the  most  important  pieces  of  information  needed  for  it  is  to 
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find  the  species  distribution  along  individual  optical  paths  in  the  hemispherical 
volume  faced  by  the  location  where  heat  transmission  is  considered.  Two  tasks 
have  to  be  carried  out  for  the  objective  (refer  to  Fig.  1):  The  combustion  product 
distribution  has  to  be  given  with  respect  to  a  reference  location  in  the  combustion 
chamber,  most  probably  the  injection  nozzle  (C  in  Fig.  1);  a  convenient  method  of 
finding  the  species  distribution  along  optical  paths  with  respect  to  D.  In  conducting 
these  tasks,  it  is  necessary  to  find  distribution  functions  that  can  be  coupled  with 
Eq.  (4)  for  attempting  to  find  its  solution. 

When  Eq.  (4)  was  studied,  with  reference  to  experimental  data 
and  computational  results  as  to  in-cylinder  combustion  product  distribution  reported 
by  others,  it  was  found  that  several  needs  for  the  goal  would  be  met,  if  the  distribu¬ 
tion  with  respect  to  the  nozzle  is  expressed  in  the  following  form  of  new  equation, 

F  -  f  exp  (-ap-b08-cz2)  (6). 

where  F  represents  either  fuel/air  ratio,  or  soot  concentration,  or  CO 2  of  H?0 
concentrations;  f,  a,  b  and  c  are  constants  to  be  determined  by  using  experimental 
method  or  computational  means.  This  new  form  of  equation  is  considered  to 
reasonably  well  describe  the  species  distribution  of  nonaxisymmetric  plumes  and  of¬ 
fer  various  advantages,  e.g.,  that  the  expression  can  be  used  for  describing  the  spray 
plume  in  swirl  motions  by  including  its  effect  in  -direction,  and  again,  that  it 
facilitate  solutions  of  Eq.  (4)  as  mentioned  in  the  previous  section.  For  this, 
however,  with  a  given  species  distribution  with  respect  to  the  nozzle,  Eq.  (6),  it  was 
necessary  to  find  the  distribution  along  individual  optical  paths  (refer  to  Fig.  1). 

In  order  to  meet  this  need,  a  new  coordinate  transformation 
method  was  developed  in  the  present  study.  Since  the  details  of  derivation  are 
shown  in  Appendix- VII,  our  new  results  are  only  shown  in  the  following. 

F  -  f  ,exp(-[(r-r  )/r  ]2  (7)- 

0  d  o  w 

where,  fd,  rQ  and  rw  are  functionally  related  to0,{ ,  pd,  i^,  zd,  f,  a,  b,  c  and  c. 
Basically,  when  the  distribution  of  combustion  products  in  a  chamber  with  respect 
to  a  reference  location,  the  distribution  can  be  calculated  along  any  chosen  optical 
path  in  the  chamber  so  that  the  the  distribution  of  volume  absorptance  along  an  op¬ 
tical  path  is  computed,  according  to  Eq.  (3),  as 

<Q  -  <1  exp  (-c1  (r-r1  )2)  (8). 

This  relation  is  used  for  implementing  Eq.  (4)  for  predicting  radiation  heat  transfer 
through  a  combustion  chamber  by  integrating  Eq.  (5)  over  the  entire  hemispherical 
volume  faced  by  the  location  of  heat  transmission. 


B.  NEW  TECHNIQUES  FOR  EQUATION  OF  RADIATION  HEAT  TRANSFER 

Several  new  techniques  have  been  developed  during  the  process 
of  developing  a  new  model  of  radiative  heat  transfer,  Eq.  (4).  Among  those  worthy 
of  mention  are  a  new  integral  function,  a  new  blackbody  function,  3nd  a  new  em¬ 
pirical  equation  for  adiabatic  flame  temperature  of  fuel/air  systems. 
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A  NEW  INTEGRAL  FUNCTION  FOR  RADIATION  CALCULATIONS. 
During  the  process  of  finding  a  solution  of  Eq.  (4),  we  encountered  a  new  integral 
function  in  the  following  form. 


c 

s 


X 

«x-l 


dx 


(9). 


Since  an  analytical  solution  for  the  equation  in  not  known,  it  was  considered  to  be 
very  difficult  to  solve  Eq.  (4).  Fortunately,  an  exact  solution  of  Eq.  (9)  was  dis¬ 
covered  in  a  more  general  form  as, 

/*2  -ax  a 

j  dx  -  Ptm.a.x^  -FCa.a.Xj).  (10) 

m  a  \ 
where,  F  (m,a,x)  *  t  /.-(«♦•>*  r 

nii  i  iio  r("-l+l)  1 


Note  that  a  special  form  of  Eq.  (9)  is  Debye  functions  of  which 
values  are  included  in  most  of  the  mathematics  tables.  That  is,  we  can  find  a  solu¬ 
tion  in  closed  form  of  Debye  function  as  follows. 


<7.1 


H  tr<3>  , 

4  u,i 


)  -  F(3,  0,x1(i>l 


(11). 


Details  of  the  derivation  of  this  new  integral  function  are  explained  in  Appendix- 
VIII. 


BLACKBODY  RADIATION  FUNCTIONS.  Among  the  more  significant  of 
our  new  integral  functions  was  another  new  discovery  of  blackbody  radiation  func¬ 


tions  by  using  Eq.  (10),  as  defined  in  the  following. 

fi 
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(12) 


Note  that  the  functions  are  tabulated  in  most  basic  heat  transfer  textbooks.  In 
addition,  for  the  computer  solution  of  various  types  of  radiation  problems,  previous 
researchers  employed  polynomial  curves  fitted  to  the  function  (Pivovonsky  and 
Nagel  (1961).  With  the  new  exact  solution  of  the  abovementioned  integral  equation, 
it  was  found  that  the  blackbody  fuctions  can  be  solved  in  a  closed  form  as  shown 
below. 


o-XT 


(XT 


err 


d (XT)  -  l  [  -2—  (x3  +  -2lC  +  ii  +  ± 


n-1 


)  1 


(13) 


The  further  details  of  the  solution  derivation  are  shown  in  Appendix-X. 


EMPIRICAL  EQUATIONS  FOR  ADIABATIC  FI.AME  TEMPERATURES. 

One  of  the  most  time-consuming  processes  in  computation  of  the  new  model  was  to 

compute  the  adiabatic  flame  temperature  (AFT),  which  was  made  based  on  fuel/air 
ratios.  Primarily,  in  an  attempt  to  reduce  the  computer  time  in  the  beginning,  a 


[»M>] 


In  the  analysis,  a  new  equation  was  developed  to  estimate  the  AFT  for  fuel-air  sys¬ 
tems  of  several  families  of  fuel  under  varied  reaction  pressures  as  shown  in  the 
following. 

T„  =  ti[\  +  h  ■  In +  c  ■  (In </>)*]  ( ' K)  (14) 

where  a  *  a\  +  «_•  •  In (P:,P») 
h  m  />|  +  />-i  •  ln(  />//*n) 

<•  =  n  +  ■  ln(  PI  Pa) 

Fuel  volume  fraction 

(Fuel  volume  IntctionlsmirH 
P0  m  I  atm 

and  aj,  bj,  ...  are  coefficients  of  each  family  of  fuels  separately  given  as  shown  in 
Appendix- XI. 

This  formula,  however,  does  not  include  the  variation  of  ini¬ 
tial  temperature  of  fuel/air  mixture.  Upon  extensive  study  of  the  AFT  of  fuel/air 
systems,  a  better  formula  for  computing  AFTs  was  developed  by  including  all  the 
major  fuels  and  variables  as  shown  below. 

7"n  =  /ti(l  +  A*  In  4>  +  ,4j(ln 

A,  »  An  +  An  In  Z’  +  ztoOn  P)1  (t-  1.2,3)  (15) 

A„  -  A„\  +  A„i  InT+^oallnfF  ( j-  1.2.3) 

A(,k  m  Aim*  Aim  In  Me  +  -ti/vaO0  'Ve)2  (k  =*  '.2.3) 

where,  A were  found  for  individual  fuel  families  and  also  fuels.  Considerations 
given  in  its  development  and  the  abovementioned  constants  for  fuels  are  given  in 
Appendix- XII. 

When  the  AFT  is  computed  for  a  mixture  of  the  initial  tem¬ 
perature  of  298°  C,  use  of  Eq.  (12)  may  be  more  convenient;  and  for  mixtures  at 
temperature  other  than  298°  C,  use  of  Eq.  (13)  is  recommended. 


C.  RELATIONSHIP  OF  GASEOUS  AND  PARTICULATE  CONCENTRA¬ 
TIONS  TO  RADIATION  HEAT  TRANSFER  IN  DIESEL  COMBUSTION 

The  objective  of  the  proposed  study  has  been  implemented  by 
conducting  parametric  analysis  of  radiative  heat  transfer  in  diesel  combustion  using 
our  new  computational  model.  Granted  that  the  new  model  remains  to  be  verified, 
some  results  from  the  present  analysis  appear  to  be  very  reasonable. 

Several  major  variables  were  considered  for  the  analysis  as 
listed  in  the  following  (An  in-depth  discussion  is  found  in  Appendix-XII).  They  are. 

(1)  The  location  in  the  combustion  chamber, 

(2)  The  surface  emissivity. 
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(3)  The  soot  concentration, 

(4)  The  soot  and  gas  radiation, 

(5)  The  reaction  pressure, 

(6)  The  combustion  chamber  wall  temperature,  etc. 

Some  of  the  more  significant  conclusions  drawn  from  the 
analysis  are  listed  in  the  following: 

o.  The  most  portion  of  radiation  heat  flux  incident  on  a  location  in  the 
chamber  wall  comes  from  the  radiation  source  right  near  the  location. 

o.  The  effect  of  surface  emissivity  on  diesel  radiation  may  not  be  negligible, 
in  particular,  when  the  surface  temperature  is  high  as  in  a  new  uncooled 
diesel  engine. 

o.  The  diesel  radiation  heat  transfer  is  highly  dependent  on  the  soot 

concentration  and  weakly  sensitive  to  the  gaseous  species  concentration. 

o.  The  pressure  effect  on  the  radiation  is  almost  negligible  in  diesel 
combustion. 

o.  The  surface  temperature  of  the  chamber  wall  is  not  a  sensitive  parameter 
on  the  diesel  radiation  heat  transfer. 


5.  SUMMARY  AND  RECOMMENDATION 

The  relationship  of  radiation  heat  transfer  to  gaseous  and  par 
ticulate  concentrations  in  diesel  combustion  has  been  studied  by  using  both  ex¬ 
perimental  and  computational  methods. 

In  the  experimental  work,  because  of  unfortunate  and  serious 
engine  failure  the  in-cylinder  measurement  was  deferred  until  a  new  engine  that  is 
being  set  up  becomes  fully  functional.  Among  the  work  tasks  which  have  been 
carried  out  in  the  experiment  are  (l)  construction  of  a  single-cylinder  direct 
injection-type  diesel  engine,  (2)  development  of  an  in-cylinder  sampling  probe,  (3) 
fabrication  of  heat  transfer  probes,  and  (4)  construction  of  sample  processing 
systems,  including  an  exhaust  dilution  unit  and  a  sample  stabilization  chamber. 
These  engine  diagnostic  tools  are  being  prepared  using  data  acquired  from  the  new 
engine. 


The  computational  work  modeling  radiation  heat  transfer  in 
diesel  combustion  has  been  greatly  successful,  to  the  present  researchers’  view,  who 
have  prepared  seven  separate  refereed  publications  on  the  subject,  as  listed  in 
section-7  and  included  in  the  appendices. 

Due  to  the  engine  problem  encountered  in  the  experimental 
approach,  the  objective  of  verifying  our  new  computatiun.il  model  has  not  ac¬ 
complished  before  the  contract  period  was  over.  It  is  highly  desirable  to  continue 
the  present  work  for  achieving  this  remaining  objective  and  a  continuation  would 
make  it  possible  to  achieve  the  objective  in  the  most  economical  way  in  time  and 
cost  because  we  have  trained  ourselves  and  acquired  tools  needed  for  the  study. 
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Eight  separate  papers  have  been  prepared  under  the  present 
sponsorship  during  the  three- year  contract  period.  They  are  included  in  the 
Appendices.  The  title  of  the  papers  are  listed  in  the  following. 

1.  "Computation  of  Radiation  Heat  Transfer  in  Diesel  Combustion", 

SAE  Transactions,  1984,  pp.  3.1002-3.1016. 

2.  "Adiabatic  Flame  Temperature  Estimates  of  Lean  Fuel/Air  Mixtures", 
Combustion  Science  and  Technology,  1983,  vol.  35,  pp.  203-206. 

3.  "A  Useful  Integral  Function  and  Its  Application  in  Thermal  Radiation 
Computations,"  Int.  Comm.  Heat  and  Mass  Transfer,  1983,  vol.  10,  pp.  329-33 

4.  "Blackbody  Radiation  Functions," 

Int.  Comm.  Heat  and  Mass  Transfer,  1984,  vol.  1 1,  pp.  451-455. 

5.  "Empirical  Equations  for  Adiabatic  Flame  Temperatures  for  Some  Fuel-Air 
Combustion  Systems,"  Combustion  Science  and  Technology,  1985, 

vol.  44,  pp.  75-88. 

6.  "Coordinate  Transformation  Method  for  Radiation  Heat  Transfer  Prediction  in 
Soot  Laden  Combustion  Products,"  being  reviewed  for  publication. 

7.  "An  Analytical  Method  of  Non-Gray  Three-dimensional  Radiation  Heat 
Transfer  in  Spray  Combustion,"  being  reviewed  for  publication. 

8.  "A  Parametric  Analysis  of  Radiation  Heat  Transfer  in  Diesel  Injection  Diesel 
Combustion,"  Proceedings  of  Int’l  Symposium  on  Diagnostics  and  Modeling  of 
Combustion  in  Reciprocating  Engines,  Sept.  4-6,  1985,  Tokyo,  pp.  167-174. 
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APPENDIX -I 


CONSTRUCTION  OF  A  SINGLE  CYLINDER  LABORATORY  DIESEL  ENGINE 

Construction  of  the  system  was  a  continuing  endeavor  for  our 
engine  research.  Several  main  engine  components  were  newly  designed/constructed 
for  the  construction  of  the  system.  Details  of  the  components  are  explained  in  the 
following. 

RECIPROCATING  UNIT  AND  CRANKCASE.  Among  the  main  design  and 
modification  details  are, 

-a  large  hole  (6  in.  diameter)  bored  at  the  top  of  the  crankcase, 

-a  new  camshaft  for  replacing  the  existing  one  in  the  ASTM  engine  to  connect 
a  fuel  injection  pump  and  an  encoder, 

-balancing  of  the  reciprocating  system  by  installing  two  balancing  weights 
(1,109  grams  each),  etc. 

The  sealing  at  both  ends  of  the  camshaft  was  constructed  a  few 
times  for  eliminating  excessive  leakage  found  during  engine  operation.  The  instal¬ 
lation  of  balancing  weights  at  the  end  of  the  connecting  rod  and  balancing  piston 
was  rather  involved  due  to  the  space  limitation  in  the  crank  case:  a  new  oil  strainer 
was  fabricated  at  a  new  location  for  reserving  space  for  their  installation.  Of 
greater  concern  in  the  reciprocating  unit  was  the  endurance  of  the  timing  gear  be¬ 
cause  the  gear  took  a  heavier  load  to  the  drive  injection  pump  unlike  the  simple 
load  undergone  when  operated  as  a  ASTM  gasoline  engine:  the  gear  was  repeatedly 
broken  three  times  requiring  replacement  of  both  the  timing  gear  and  the  oil-pump 
driving  gear.  It  was  found  that  the  original  crank  shaft  might  not  been  suitable  for 
being  used  in  the  new  diesel  engine  as  we  discovered  that  the  unit  was  bent  to  one 
direction  creating  excessive  vibrations  of  the  engine. 

CYLINDER  BLOCK.  A  new  cylinder  block  was  constructed  where  the 
cylinder  liner  was  installed  and  a  water  jacket  was  attached  for  cooling  the  liner 
(the  cooling  of  the  cylinder  head  was  made  through  a  separate  coolant  flow).  The 
construction  details  are, 

-a  cylinder  block  that  houses  a  cylinder  liner  from  GM-53  series, 
-accommodation  of  corresponding  engine  components  including  liner,  piston, 
piston-ring,  con-rod,  bearings  with  proper  modifications  where  needed. 

-a  newly  fabricated  water  jacket  package  installed  on  two  cavities  of  the  new 
cylinder  block, 

-new  intake  and  exhaust  manifolds  connected  to  the  engine  block,  etc. 

The  construction  of  the  cylinder  block  was  a  difficult  task 
since  the  proper  scavenging  of  the  engine  was  required  in  the  system.  An  approach 
taken  for  this  goal  was  to  install  several  different  newly  fabricated  vane  blocks 
before  the  intake  ports  in  order  to  identify  the  most  suitable  one  for  the  engine: 
this  was  basically  to  minimize  the  soot  emission  from  the  engine.  In  this  experience, 
a  swirl  meter  was  installed  in  the  cylinder  under  steady  flow  conditions  to  find  the 
highest  swirl  motion  of  the  intake  air,  since  the  swirl  of  the  engine  configuration 
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was  very  low  due  to  the  simple  design  of  the  engine  components.  An  aluminum 
block  carved  for  air  flow  seemed  the  most  suitable  of  all  we  tried.  Upon  extensive 
engine  operation  periods,  we  found,  however,  a  severe  deposit  of  wet  soot  on  the  in¬ 
take  port  indicating  significant  back  flow  of  air  during  the  early  part  of  the  ex¬ 
haust  period.  Since  the  fuel  was  injected  every  two  cycles  providing  an  extra 
scavenging,  the  construction  of  vane  blocks  was  not  further  pursued.  However,  this 
approach  might  have  been  responsible  for  the  deposit  of  wet  soot  due  to  over¬ 
cooling  of  the  engine  components. 

CYLINDER  HEAD.  The  construction  of  the  engine  cylinder  head  was  another 
unique  experience.  The  main  head  was  constructed  out  of  a  cold-roll  block  includ¬ 
ing  the  following  machinings: 

-a  cylindrical  cavity  for  the  combustion  chamber  by  matching  with  the 
cylinder  liner  and  head  gasket  (Fig  A-l), 

-a  water  jacket  around  the  cavity  connecting  the  coolant  flow  separate  from 
that  in  the  cylinder  block  (Fig  A-2), 

-a  mounting  plug  and  its  cavity  in  the  cylinder  head  for  in-cylinder  probes 
enabling  selection  of  measurement  locations  during  the  engine  operation, 
-installation  of  the  injection  nozzle  with  minimum  use  of  space  in  the  head 
for  reserving  maximum  room  for  the  in-cylinder  probe  mounting  as 
mentioned  above, 

-three  injection  nozzles  with  different  sac  volumes  which  were  specially 
manufactured  for  the  new  engine  experiment  (Fig  A-3),  and 
-brackets  for  mounting  thermostat  and  coolant  connections. 

Two  of  the  main  concerns  in  its  construction  were  the  cooling 
of  the  cylinder  head  and  selection  of  the  injection  nozzle  suitable  to  the  engine.  It 
was  found  that  the  cylinder  head  was  excessively  cooled  in  the  early  tests.  A  new 
thermostat  was  installed  for  eliminating  the  problem.  The  selection  of  the  injection 
nozzle  still  remains  as  an  unresolved  experience  due  to  the  insufficient  detailed 
measurements  and  relatively  short  period  of  experiment  with  the  engine;  we  have 
been  using  one  type  of  nozzle  of  three  selections  for  eliminating  various  immediate 
engine  design  issues. 

ENGINE  ACCESSORIES.  Several  new  engine  accessories  were  required  in  the 
engine  apparatus  for  facilitating  the  engine  operation  and  control. 

-a  mounting  plate  (24-24-0.5  inches)  installed  in  the  front  side  of  the  engine 
for  mounting  various  engine  components:  a  fuel  injection  pump  (American 
Bosch,  APE  1B50P5625B);  a  water  pump;  a  fuel  filter  and  water  separator 
assembly;  a  phase  shifter  (for  the  timing  control  of  in-cylinder  sampling); 
installation  of  an  encorder;  etc., 

-a  new  phase  shifter  which  is  to  vary  the  relative  phase  angle  between  the 
engine  camshaft  and  the  sampling  system  enabling  the  shift  of  timing 
during  the  engine  operation, 

-an  air  induction  unit  for  controlling  the  inlet  air  pressure  and  f'ow  rate. 

-a  closed-loop  coolant  unit  including  devices  for  control  of  its  flow  rate  and 
temperature:  control  thermostat;  magnetic  valve;  water  pump;  electric  heater, 
etc.,  and 
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-an  exhaust  gas  system  including  a  surge  tank  where  exhaust  samples  were 
taken. 


The  new  coolant  system  needs  more  explanation.  It  was 
required  to  operate  the  engine  for  an  extensive  period  before  achieving  its  warm-up 
condition.  It  seemed  undesirable  to  operate  the  laboratory- built  engine  for  such 
repeated  periods  of  time  particularly  when  the  system  had  some  unresolved  design 
problems  including  various  stc  '  items  from  various  sources  of  which  mutual  com¬ 
patibility  was  not  established.  In  order  to  minimize  the  warm-up  period,  a  separate 
heating  unit  was  employed  to  heat  the  coolant  and  to  circulate  it  through  the  engine 
so  that  the  entire  engine  block  attained  an  elevated  temperature  before  the  firing  of 
the  engine. 


APPENDIX -II 

DEVELOPMENT  OF  SAMPLING  PROBE 

When  the  research  proposai  for  the  present  work  was 
submitted,  the  principal  investigator  found  more  than  thirty  major  technical  papers 
prepared  (as  early  as  the  1920s)  reporting  results  obtained  by  using  in-cvlinder  sam¬ 
pling  probes  for  engine  studies.  Some  more  papers  appeared  since  then  to  add  to  the 
archives  of  the  sampling  method.  The  sampling  valves  used  in  those  studies  were 
mostly  of  the  intermittent  flow-type  where  the  sampling  hole  is  kept  normally 
closed  and  rapidly  opened  for  a  short  period  of  time  to  take  a  small  amount  of  in¬ 
cylinder  gas  volume  at  the  desired  time  of  sampling.  This  type  of  valves  were  con¬ 
sidered  to  have  some  problems  in  obtaining  reliable  sample  compositions  (Zvonov 
et  al.,  1968;  Whitehouse  et  al.,  1977)  due  to  the  quenched  layer  which  might  exist  in 
the  close  vicinity  of  the  sampling  hole.  In  an  attempt  to  minimize  such  a  potentially 
severe  problem,  Rhee  et  al.  (1978)  developed  a  new  continuous  sampling  system. 

The  new  system  employs  a  continuous  stream  of  in-cylinder  gas 
flow  through  the  sampling  hole  (Fig  A-4).  The  flow  is  rapidly  diverted  for  a  short 
period  of  time  to  sample  volume  for  the  subsequent  analysis.  Briefly  mentioning  the 
operation  of  the  new  system,  a  flow  passing  through  the  orifice  in  the  probe  (Fig  A- 
4  (A))  is  normally  connected  to  a  large  (waste  gas)  section  of  the  oscillation  tube 
(Fig  A-4  (C))  which  is  placed  in  the  probe  under  spring  load.  When  the  sample  is  to 
be  obtained,  the  sampling  tube  is  rapidly  rotated  through  a  sectoral  range  as  indi¬ 
cated  by  the  broken  lines  in  the  figure.  During  this  period  the  sampling  hole  is 
briefly  swept  by  the  sample  tube  by  taking  a  small  amount  of  gas  out  of  the  con¬ 
tinuous  flow  stream.  The  sampling  in  the  next  cycle  is  made  by  moving  the  oscilla¬ 
tion  tube  in  the  opposite  direction  for  a  new  sampling.  This  actuation  motion  of  the 
tube  was  made  by  using  a  hydraulic  piston-cylinder  unit  controlled  by  a  four-gate 
valve  that  was  operated  in  phase  of  the  phase  shifter  connected  to  the  engine  (Fig 
A-5).  The  hydraulic  force  was  provided  from  a  separate  high-pressure  accumulator. 

Of  more  significance,  the  new  system  constructed  in  the  cur¬ 
rent  study  was  to  include  the  dilution  tube  in  the  sample  tube  so  that  a  dilution  gas 
having  a  non-combustion  product  trace-gas  is  injected  into  the  sample  tube  for 
quenching  any  physical  and  chemical  reactions  that  potentially  take  place  in  the 
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system.  Due  to  the  limited  resources,  only  one  student  was  assigned  to  the  ex¬ 
perimental  work  (further,  there  were  unexpectedly  large  amounts  of  work  needed  in 
the  engine  construction)  so  that  the  usefulness  of  the  new  sampling  system  has  not 
yet  been  evaluated.  It  will  be  very  interesting  to  systematically  evaluate  the  new 
system  in  an  engine  in  our  continuing  study. 


APPENDIX-111 

HEAT  TRANSFER  PROBES 

There  have  been  two  views  of  heat  transfer,  in  particular 
radiation  heat  transfer,  that  potentially  determines  iiesel  economy:  "Since  the  com¬ 
bustion  flame  in  the  diesel  engine  is  luminous  (but  not  in  the  SI  engine),  it  has  been 
speculated  for  many  years  that  radiation  heat  loss  may  be  an  important  factor  of 
diesel  economy"  (Obert,  1973).  This  view  has  the  same  emphasis  of  the  portion  of 
radiation  out  of  the  total  in-cylinder  heat  radiation  as  that  by  Ebersole  et  al.  (1963) 
who  reported  that  the  portion  by  radiative  process  was  as  high  as  40%  of  the  total 
measured.  On  the  other  hand,  Amann  (1980)  reported  his  computation  suggesting 
that  the  improvement  in  thermal  efficiency  might  be  only  a  few  percent  with  up  to 
85%  of  recovery  in  heat  loss  through  the  water  jacket  and  lubricating  system  of 
diesel  engines.  This  view  is  supported  by  a  second  law  analysis  of  an  engine  by 
Primus,  et  al.  (1984).  Although  more  data  on  radiation  heat  transfer  is  neededfor  a 
better  evaluation  of  its  impact  on  engine  economy  since  the  loss  through  the  radia¬ 
tive  process  when  the  compression  ratio  is  very  high,  the  importance  of  the  heat 
transfer  clearly  resides  in  the  thermal  loading  of  engine  components  of  the  combus¬ 
tion  chamber.  Particularly  since  the  radiative  heat  transfer  is  highly  direction- 
dependent  and  since  the  fate  of  critical  components  of  the  advanced  engines,  are 
affected  by  the  thermal  load,  the  further  study  of  heat  transfer  analysis  is 
warranted. 

The  basic  details  of  the  present  probes  are  the  same  as  the 
original  design  as  shown  in  Fig  A-6.  The  Figure  shows  the  probes  mounted  in  the 
cylinder  head:  one  is  to  measure  the  total  (convective  plus  radiative)  heat  transfer 
and  the  other  is  to  measure  the  portion  of  radiative  heat  transfer  through  the  com¬ 
bustion  chamber  wall.  The  measurement  of  the  radiative  portion  of  heat  transfer 
was  to  be  obtained  by  eliminating  the  convective  portion  out  of  the  total  heat 
transfer.  The  detector  is  ASEA  double  thermocouples  which  enable  measurement  of 
unsteady  heat  flow  with  rising  time  of  up  to  2  microseconds.  The  main  weakness  of 
this  method  are  that  there  are  soot  fouling  formation  on  the  window  and  that  the 
recessed  mount  of  the  probe  needs  a  proper  evaluation  of  its  geometric  factor.  The 
evaluation  of  the  geometric  factor  was  attempted  by  using  the  flash  bulb  technique. 
The  method  is  to  employ  flash  bulbs  used  in  photography  that  produces  a  relatively 
strong  radiation  with  designated  intensities  by  the  manufacturers.  After  several 
trials  by  using  a  bench  rig,  this  method  was  found  unsuitable  due  to  the  excessive 
variation  of  radiation  intensity  of  the  flash  bulbs  overtaking  the  effect  of  the 
recessed  mount. 

Consequently,  we  decided  to  directly  evaluate  the  effects  in 
the  engine  as  Ebersole  et  al.  (1963).  Some  engine  measurements  were  made  by  using 


the  probe  in  the  engine,  which  need  further  processing  prior  to  analysis.  This  task 
has  not  been  made  because  of  the  following  reason.  While  the  construction  of  the 
engine  setup  and  the  sampling  probe  was  carried  out  by  a  graduate  student  (J. 
Matonick)  (note  that  the  remaining  graduate  assistantship  was  given  to  S.L.  Chang 
for  the  computational  work),  the  investigation  of  heat  transfer  was  conducted  by 
another  student  (C.  Sheft)  on  a  Departmental  support.  At  the  time  of  writing  the 
thesis,  C.  Sheft  accepted  a  job  offer  from  industry  and  decided  to  complete  the 
work,  which  is  being  carried  out. 
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Fig.  A-4  The  Continuous  Flow-Type  Sampling  System  with  Dilution 


Fig.  A-5  The  Sampling  Valve  Actuation  System 


Appendix-IV.  Computation  of  Rad.  Heat  Transfer  in  Diesel  Combustion 


831332 

Computation  of  Radiation  Heat 
Transfer  in  Diesel  Combustion 

S.  L.  Chang  and  K.  T.  Rhee 

Dept,  of  Mechanical  and  Aerospace  Engrg. 

Rutgers  Univ. 


ABSTRACT 

A  theoretical  model  of  radiation  heat  trans¬ 
fer  has  been  developed.  A  computation  of  radia- 
cion  heat  flux  at  a  particular  location  in  the 
combustion  chamber  by  using  the  present  model 
requires  in-cylinder  time-  and  space-resolved 
species  data  and  cylinder  pressure.  From  the 
species  data,  the  burned  fuel/air  ratio  distri¬ 
bution  is  inferred  to  compute  space-resolved 
adiabatic  flame  temperature.  For  the  computation 
of  the  spectral  emissivity  of  an  isothermal  vol¬ 
ume  of  adiabatic  temperature  containing  soot,  the 
Rayleigh-limit  expression  is  used.  The  refrac¬ 
tion  indices  in  the  expression  are  obtained  by 
using  the  dispersion  equations  based  on  the  elec¬ 
tronic  theory  encompassing  both  free  and  bound 
electrons.  For  che  spectral  emissivity  from  the 
gaseous  component  in  che  volume,  the  semi-empiri¬ 
cal  band  model  is  used.  A  parametric  analysis  of 
radiation  heat  transfer  in  diesel  combustion  is 
made  by  using  Che  present  model;  a  prediction  by 
Che  model  qualitatively  compares  with  some  of  che 
reported  experimental  data. 


INTRODUCTION 

The  diesel  combustion  engine  produces  high 
power  output  per  unit  voLume  of  fuel.  Because  of 
the  ever  growing  shortage  of  petroleum,  che  die¬ 
sel  engine  has  been  recognized  to  be  a  more  im¬ 
portant  means  for  transportation,  which  is  well 
reflected  in  che  recent  increased  number  of  die¬ 
sel  engine-powered  vehicles. 

One  of  che  concerns  in  the  wide  use  of  die¬ 
sel  engines  is  the  emission  of  diesel  smoke  (or 
particulate).  The  diesel  smoke  not  only  looks 
bad,  buc  also  contains  polynuclear  aromatic  hy¬ 
drocarbon,  some  of  which  are  known  co  cause  muta¬ 
genic  actions  on  living  cells  (U*.  In  addition, 

♦Numbers  in  parenthesis  designate  references  at 
end  of  paper. 


in  spite  of  the  high  fuel  economy  of  diesel  pow¬ 
ered  vehicles,  there  is  a  consensus  that  it  can 
be  further  improved.  Since  the  conventional  die¬ 
sel  engines  reject  as  much  as  ->0%  of  che  total 
energy  produced  by  the  combustion  reaction  to  che 
cooling  and  lubrication  system,  one  is  interested 
in  reducing  che  heat  loss  for  improving  the  en¬ 
gine  system.  Unlike  che  premixed  spark-ignition 
engines,  the  flames  in  diesel  combustion  are 
highly  luminous.  Although  it  varies  from  engine 
to  engine,  che  portion  of  radiation  heac  loss  out 
of  the  total  heat  transfer  in  diesel  engines  has 
been  reported  co  range  up  co  30%  or  more  (2-5). 
Owing  to  the  fact  that  the  concentration  of  the 
particulates  appearing  in  che  combuscion  flame  is 
closely  related  co  its  radiation  heat  emission, 
without  mentioning  other  complex  processes  af¬ 
fecting  the  efficiency,  the  large  amount  of  smoke 
appearing  in  the  exhaust,  therefore,  mav  suggest 
a  negative  sign  in  diesel  economy. 

Regarding  che  engine  improvement,  it  was 
suggested  (6)  the  thermal  efficiency  of  diesel 
engines  would  gain  a  mere  few  percenc  with  up  to 
85%  of  recovery  in  heac  loss  through  the  water 
jackec  and  lubrication  system  of  diesel  engines 
by  insulating  the  combustion  chamber.  The  insul¬ 
ation  was  proposed  to  achieve,  bv  using  ceramic 
components  in  che  engine,  an  adiabatic  engine 
concept  (7).  A  successful  development  of  an  ici- 
abatic  diesel  engine  will  eliminate  the  cumber¬ 
some  engine  cooling  svstem,  e.g.,  che  water  jack¬ 
et  and  the  radiator  making  the  svstem  more  com¬ 
pact  and  tidy,  an  additional  factor  in  improving 
fuel  economy  (8).  A  prototype  engine  of  the  adi¬ 
abatic  concepc  has  exhibited  a  specific  fuel  con¬ 
sumption  race  of  0.28  Ib/hp-hr  (3).  Furthermore, 
when  a  turbocharger  is  combined  with  the  adiabat¬ 
ic  engine  (as  a  "compound  engine’"' ,  the  overall 
efficiency  of  the  engine  will  be  increased  i i ,  1 '  . 

It  appears  to  be  clear  chat  the  detailed 
understanding  of  heat  transfer  in  diesel  combus¬ 
tion,  a  precondition  in  the  improvement  ot  diesel 
engines  in  various  aspects,  requires  extensive 
information  as  co  the  in-cvl inder  species  distri¬ 
bution  during  che  reaction  period.  In  irder  to 
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help  meet  such  needs,  research  of  diesel  combus¬ 
tion  has  been  conducted  at  Rutgers,  The  State 
University  of  New  Jersey.  The  study  is  to  inves¬ 
tigate  the  processes  of  diesel  particulate  for¬ 
mation  and  radiation  heat  transfer  in  direct  in¬ 
jection-type  diesel  engines  by  achieving  the 
following  objectives:  (1)  to  obtain  time-  and 
space-resolved  gaseous  and  particulate  concen¬ 
tration  in  direct  injection-type  diesel  engine 
combustion,  (2)  to  measure  instantaneous  heat 
transfer  (convective  and  radiative,  respective¬ 
ly)  through  the  combuscion  chamber  wall,  and  (3) 
to  identify  the  relationship  of  the  in-cylinder 
species  concentrations  to  the  radiative  heat 
transfer  by  constructing  a  theoretical  radiation 
heac  transfer  model.  The  present  paper  reports 
the  theoretical  model  of  radiation  heac  cransfer 
developed  during  the  course  of  che  above  men¬ 
tioned  diesel  combuscion  research.  A  computation 
of  radiation  heat  cransfer  in  diesel  combuscion 
by  using  che  present  model  requires  che  in-cylin- 
der  species  distribution  and  cylinder  pressure. 

APPROACH  OF  THEORETICAL  MODEL 

Among  in-cylinder  information  needed  in  che 
computation  of  the  radiative  energy  incident  over 
a  particular  location  of  the  diesel  combustion 
chamber  wall  are  che  distribution  of  species  and 
temperature,  pressure,  and  geometric  factor.  In 
che  computation,  it  is  necessary  to  find  che 
time-resolved  spectral  absorption  coefficient  of 
each  volume  of  combustion  products  throughout  the 
chamber,  since  che  reaction  process  is  highly 
nonuniform  in  nature  with  both  locacion  and  time. 

Since  Kunitomo  et  al.  (10)  reported  a  model 
to  predict  radiation  heat  flux  in  diesel  combus¬ 
cion  chat  employed  a  gray  body  of  gaseous  and 
particulate  mixtures  with  a  uniform  concentration 
so  as  to  exclude  che  geometric  factor,  much  re¬ 
search  has  been  carried  out  on  radiation  heac 
cransfer  of  flames  enabling  a  more  extensive  com¬ 
putation:  che  in-cylinder  sampling  of  diesel 

combustion  shed  a  light  to  offer  in-depth  species 
information  in  the  combustion  chamber  (11,12);  a 
comprehensive  technique  was  proposed  to  compute 
thermal  radiation  by  gaseous  mixture  (13);  the 
method  of  computing  che  optical  constants  of  soot 
mixture  has  greatly  been  improved  (14,15);  che 
computation  of  chemical  equilibrium  properties  of 
combustion  becomes  feasible  (16,17,18). 

The  present  model  is  an  attempt  to  combine 
Chose  improved  techniques  for  the  computation  of 
radiation  energy  incident  over  individual  loca¬ 
tions  of  Che  combustion  chamber  wall.  In  che 
computation,  the  radiation  from  a  small  isother¬ 
mal  volume  of  combustion  products  is  assumed  to 
be  Che  superimposed  result  of  the  emissions  from 
sooc  and  gaseous  mixtures  without  mutual  inter¬ 
ference.  The  emissivity  of  soot,  -;g,  may  be  ex¬ 
pressed  in  terms  of  che  volumetric  absor^tance 
of  soot,  <s,  and  che  optical  length,  L,  deter¬ 
mined  bv  the  size  of  che  volume,  as  follows. 


Similarly,  the  emissivity  of  gases,  £g,  can 
be  shown  by  using  che  volumetric  absorptance  of 
gases,  «rg. 


Since  it  is  assumed  that  che  total  absorp¬ 
tion  coefficient,  < ,  is  the  sum  of  <s  and  <g, 
one  may  write  the  emissivity  of  the  combined  ra¬ 
diation,  e,  as. 


e  -  1  -  (1-e  )(l-e  )  (1) 

s  g 

Therefore,  the  radiation  energy  streaming 
away  from  the  volume,  B,  by  knowing  its  tempera¬ 
ture,  T,  can  be  found  as. 


B  =  -;oT4  (2) 

where,  j  is  the  Stefan-Bolczman  constant. 

The  incident  radiation  flux,  H,  on  a  partic¬ 
ular  detecting  area,  Aj ,  originated  from  the 
emitting  source  area,  Ae ,  can  be  expressed  as. 


A  —A 
e  d 


F  -  oT 

A  -A. 
e  d 


(3) 


where,  F^  is  the  geometric  factor  from  A  to 
e  d 

Ad‘ 

Integrating  over  the  entire  emitting  source, 
one  can  find  che  tocal  radiation  energy  incident 
on  the  detecting  area,  chat  is, 


V  -  i.A  FA  -A  -T  dA 

d  e  e  d 


(4) 


The  result  of  che  present  compucacion  is 
reported  by  using  the  "apparent  emissivity"  de¬ 
fined  as. 


H,  :T 
Ad  r 


(5) 


where,  T  is  arbitrarily  chosen  to  be  2400  K. 
r 


COMPUTATIONAL  DETAILS 


The  details  of  che  computational  steps  and 
assumptions  of  the  methodology  outlined  above  are 
explained  for  che  following  items:  (1)  the  tem¬ 
perature  of  each  mixture  volume.  (2)  the  soot 
radiation,  (3)  the  gas  radiation  and  (-*)  a  compu¬ 
cacion  of  radiation  heat  transfer  in  diesel  com¬ 
bust  ion  . 

Adiabatic  Flame  Temperature  -  The  tempera¬ 

ture  of  a  designated  volume  of  in-ev  1  inder  mix¬ 
ture  was  assumed  to  attain  the  adiabatic  flame 
temperature  which  can  be  computed  bv  using  the 
species  distribution  which  mav  be  obtained  either 
theorec ical  1 v  or  experimental Iv .  In  iddition, 

the  encitv  in  the  volume  is  issumed  to  be  it  in 
eauilibrium  of  "he  temperature  ind  pressure. 
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Besides  fuels,  eleven  gaseous  species  were  in¬ 
cluded  in  Che  computations  of  the  adiabatic  flame 
temperature  and  che  equilibrium  composition  by 
using  the  burned  fuel/air  ratio  inferred  from  the 
species  distribution.  Although  such  computations 
have  been  widely  employed  by  others  for  flame 
studies  in  the  past  (16,17),  Che  present  paper 
offers  a  new  technique  of  their  computations  that 
greatly  reduce  che  computer  time.  The  details  of 
the  computation  and  che  new  iteration  technique 
are  shown  in  Appendix  I. 

Soot  Radiation  -  The  thermal  radiacion  of 

soot  in  diesel  combuscion  is  the  combined  result 
of  the  emission,  absorption,  and  scattering  of 
soot  particles.  Since,  according  to  Planck's 
function,  the  main  portion  of  the  thermal  emis¬ 
sion  from  the  soot  is  infrared,  i.e.  ,  in  the 
range  of  wavelength,  X  *  l-10um,  the  present  com¬ 
putation  considers  che  corresponding  range  of 
wavelength.  Regarding  the  absorption,  if  the 
designated  volume  of  mixture  attains  a  thermal 
equilibrium,  the  absorption  of  the  radiation  can 
be  determined  by  the  Kirchhoff  principle  which 
states  chat  the  spectral  volume  absorption  coef¬ 
ficient  is  exactly  equal  to  the  spectral  volume 
eraissivity  for  a  system  in  a  thermal  equilibrium 
with  its  surroundings.  In  the  present  paper,  for 
convenience,  che  absorption  coefficient  is  used 
inscead  of  the  volume  emissivity  or  absorptivity. 
Discussing  the  scattering,  Kunitomo  (19)  reported 
that  che  great  majority  of  diesel  exhaust  par¬ 
ticulate  was  so  small  that  their  scattering  af¬ 
fected  less  Chan  17.  at  a  photon  wavelength  of  lum 
and  considerably  less  at  a  larger  wavelength. 
Assuming  that  the  in-cylinder  soot  size  distribu¬ 
tion  is  in  che  similar  range  of  chat  in  the  ex¬ 
haust,  which  was  supported  by  a  recent  experi¬ 
mental  study  (20),  the  scattering  of  soot  is 
neglected  in  the  thermal  radiacion  of  diesel 
combustion. 

With  the  above  assumptions,  the  soot  radia¬ 
tion  is  computed  in  the  following  sequence:  (1) 
to  formulate  che  spectral  emissivity  of  the  soot 
-;SA,  (2)  to  find  the  spectral  absorption  coef¬ 
ficient,  <s;,  ,  which  is  relaced  to  che  refraction 
index  of  the  soot,  and  (3)  to  compute  the  soot 
emissivity . 

As  explained  earlier,  che  spectral  absorp¬ 
tion  coefficients  can  be  found  from  che  following 
relationship , 


of  the  soot,  and  fv  is  che  soot  volume  fraction. 

The  above  expression  to  employ  dispersion 
equations  based  on  the  classical  electronic 
theory  (14,15)  was  used  in  the  present  computa¬ 
tion.  According  to  the  experimental  results, 
the  complex  optical  constants  (n-ink)  in  Eq .  (7) 
vary  with  the  wavelength  as  shown  in  Table  I. 

Table  I:  Optical  Constants  of  Soot  ( li ) 
Acetylene  Soot  Propane  Soot 


(H/C- 

1/14.7) 

(H/C* 

1/4.6) 

l(um) 

n 

nk 

n 

nk 

0.4358 

1.56 

0.46 

1.57 

0.46 

0.45 

1.56 

0.48 

1.56 

0.50 

0.55 

1.56 

0.46 

1.57 

0.53 

0.65 

1.57 

0.44 

1.56 

0.52 

0.8065 

1.57 

0.46 

1.57 

0.49 

2.5 

2.31 

1.26 

2.04 

1.15 

3.0 

2.62 

1.62 

2  21 

1.23 

4.0 

2.74 

1.64 

2. 38 

1 . 44 

6.0 

3.22 

1.34 

2.62 

1.67 

10.0 

4.80 

3.32 

3.48 

2.46 

Dalzell  and  Sarafim  (14)  suggested  chat  Che 
optical  constants  of  different  sooc  mav  depend, 
to  a  small  extent,  on  difference  in  molecular 
structure  and,  to  a  large  extent,  on  che  differ¬ 
ences  in  the  ratio  of  hydrogen  to  carbon  mole¬ 
cules  (H/C).  Note  that  che  H/C  of  che  soot  in 
most  combustion  systems  is  found  to  fall  between 
the  cwo  extremes  as  shown  in  the  table.  As  for 
che  temperacure  dependency,  it  was  found  that  the 
optical  constants  were  insensitive  to  temperacure 
(14,15).  In  providing  a  basis  for  interpolating 
the  results  of  the  table,  a  dispersion  model  was 
found  to  fic  the  data  bv  choosing  the  appropriate 
parametric  constants  (14).  More  details  of  che 
use  of  the  dispersion  equation  for  finding  the 
optical  constants  are  shown  in  Appendix  II. 

Upon  the  given  values  of  fuel  air  ratio  ind 
fv,  the  spectral  emissivitv  of  soot,  ■:  s  .  ,  or  a 
volume  can  be  determined  bv  using  che  above  men¬ 
tioned  optical  constants  of  the  soot.  Fr>m  che 
definition,  the  soot  emissivity,  -:s,  is  expressed 
as , 

i  =  I  -t  .  £,  (T) d X /  (3) 

s  o  s-i  b ' 


si 


(6) 


Since  the  median  size  of  soot  is  in  the 
range  of  O.lum  (21,22),  which  is  approximately 
one  tench  of  infrared  wavelength,  <s,  mav  be 
found  using  che  following  Ray le igh-limit  expres- 
s ion  ( 21 ) . 


where,  E^.(T)  is  the  Planck's  radiation  function 


E  . (T)  -  — 
b  •  .  o 


2-hc 


[  exp  (he  /  '■  KT )  —  1  ] 


*  2.998  .<  10 10  cm/ sec. 

=  1.380  x  !0~^  erg/’K,  inc 


36n^’<  (  - /  ">  )  f 

<  *  - r-L -  <  7 ) 

sa  2  2  -  a  2 

[ n“ (  l-k“ ) +2 ]  +4n  k“ 

where,  n  and  nk  are  che  complex  refraction  indices 


h  -  P.P2  5  x  10  erg.  sec. 


Supstituting  Eq  .  'hi  into  Eq .  '  3 )  ,  one  mtaans, 
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OD  —  ^  >  L  f 

es  -  1  -  /  e  E^CndA/aT4  (9) 

Many  researchers  (23-26)  attempted  to  find 
the  solution  to  Eq.  (9)  by  assuming  a  simple  ex¬ 
pression  of  the  spectral  absorption  coefficient 
as, 

<.  -  <  f  /Ad  (10) 

A  O  V 

where,  <0  is  a  constant  independent  of  tempera¬ 
ture  and  d  is  chosen  to  be  close  Co  1. 

Although  this  enabled  one  to  further  simpli¬ 
fy  Eq.  (9)  to  obtain  an  integral  equation  con¬ 
taining  the  hexagarama  function  (24,27),  the  final 
result  was  calculated  by  numerical  means. 

Even  though  their  method  was  relatively 
simple  and  widely  used,  a  better  technique  to 
find  an  exact  solution  of  Eq.  (9)  has  been  de¬ 
veloped  in  the  present  work.  Since  the  refrac¬ 
tion  indices  are  almost  constant  in  each  sub¬ 
interval  of  the  entire  integration  range,  the 
solution  to  Eq.  (9)  can  be  found  as  a  closed-form 
by  using  the  following  variable  transformation 
method.  In  order  to  replace  the  temperature  and 
wavelength  variables  in  Eq.  (9),  a  new  dimension¬ 
less  variable  x  is  defined  as  follows. 


x  -  hc/AKT  (11) 

To  obtain  a  tidy  expression,  another  parameter, 
a,  is  given  as. 


_ 36Trn  k _  KT 

[n2-n2k2+2]+4n4k2  V  ch 


Equation  (9)  ,  then,  can  be  expressed  in 
terms  of  x  and  a  as, 


,  to  r  -ax  x 

1  -  — r  e  -  dx 

4  '  x  , 

-r  o  e  -1 


where , 


S ( 3 ,a ,x)  *  l  e 
n»l 


■  ( n+a ) x 


r  ,  f(3+l)  ,  ,-i-l  3-i, 

,;0  [“(3-i+l)(n+a)  x  !  ••• 

and  T  is  the  gamma  function  defined  as  f(m  +  1)  = 
m  •  T(m)  and  T(l)  -  1. 

3y  using  the  above  solution,  the  soot  emis- 
sivity  of  an  isothermal  volume  at  a  given  temper¬ 
ature  can  be  found  as. 


1  -  =t  L  [S ( 3  ,a  .  ,x  . )  -  S ( 3 ,a  .  ,x  . ) ]  (17) 
4  .  *•  it  t  t+1 

T  1=1 


Gas  Radiation  -  Although  chere  are  many 

gaseous  species  in  flames,  carbon  dioxide  and 
water  vapor  are  reported  to  be  mainly  responsible 
for  the  thermal  radiation  (10,24,29).  Therefore, 
these  two  species  were  only  considered  in  the 
present  model.  The  gas  radiation,  depending  upon 
the  molecular  structure,  has  several  blocks  of 
che  lines  in  the  spectra  of  soot.  By  averaging 
over  che  wave  number  of  each  block,  a  continuous 
distribution  of  the  emission  (or  absorption)  co¬ 
efficient,  which  is  also  called  a  band,  can  be 
computed  by  assuming  that  waves  outside  these 
bands  neither  emitted  or  absorbed. 

The  major  bands  of  both  gases  occur,  for 
water  vapor,  of  the  principal  vibration-rotational 
bands  at  1.38um,  1.87um,  2.7um,  and  6.3um  and  of 
a  weak  rotational  band  at  20um;  for  carbon  diox¬ 
ide,  of  the  principal  bands  at  l.9um,  2.7cm, 

4.3ura,  and  15um,  and  of  the  two  minor  bands  at 
9.4um,  10.4um.  There  are  two  overlapping  bands 
at  2.7um;  one  is  from  water  vapor  and  the  other 
from  carbon  dioxide. 

By  definition,  the  spectral  band  absorptance, 
x  ,  of  an  isothermal  volume  is  written  in  a  simi¬ 
lar  form  of  Eq.  (6)  as, 


or  in  the  following  numerical  integration  form. 


a  =  1 


exp(-«.  L) 

J 


(13) 


»  1 
s 


v 

15  - 

i-1 


j-xi+l 


.  x 
i 


3 


(14) 


where,  M  is  the  number  of  che  intervals  in  the 
whole  integral  range  and  ai  is  the  parameter  a 
of  the  ich  interval. 

The  solution  of  che  integral  in  Eq .  (14)  has 
been  found  by  using  a  unique  method  developed  by 
the  authors  (28).  Without  showing  che  procedure, 
che  exact  solution  is  given  as  follows. 


where,  j  is  the  wave  number,  L  is  the  optical 
length,  and  <kJ  is  the  volume  absorption  coeffi¬ 
cient  of  the  gas. 

The  experimental  results  (30-33)  exhibited 
the  distribution  of  the  spectral  absorption  coef¬ 
ficient  to  be  rather  irregular  so  that  a  general 
expression  of  the  radiation  bands  is  unavailable. 
For  practical  applications,  several  models  have 
been  proposed  co  fit  the  experimental  data.  Of 
those,  the  wide  band  model  proposed  by  Edwards 
and  Balakrishnan  (13)  was  emploved  in  the  present 
work.  The  model  assumes  a  block  approximat ion 
that  calls  for  black-bodv  in  the  bands  and  is 
transparent  otherwise.  The  cotal  band  absorbence. 
A^,  which  is  band  width  of  the  ith  band,  is  de¬ 
fined  as. 
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where,  Is  Che  wave  number  ac  the  center  of 
the  1th  band. 

The  model  provides,  for  each  band  of  the 
gas  radiation,  a  semiempirical  correlation  of 
the  integrated  band  intensity,  3,  Che  band  width 
parameter,  A  ,  the  line  width  parameter,  3,  and 
the  effective  pressure,  Pe.  The  use  of  these 
parameters  for  finding  A  is  shown  later  in  Table 
II,  and  computation  of  each  parameter  is  shown 
in  the  following. 

The  effective  pressure,  Pe,  which  is  re¬ 
sponsible  for  the  pressure  broadening  effect,  is 
written  as, 

Pe  *  [P/P0  +  Cpi/PoXb“l)  1N  (20) 

where,  PQ  is  one  atmosphere,  P  is  the  cylinder 
pressure  of  the  engine,  P^  is  the  partial  pres¬ 
sure  of  the  absorbing  gas,  and  b,N  are  the  self¬ 
broadening  coefficients  given  in  Appendix  III. 

The  computation  of  several  band  parameters 
(e.gi,  integrated  band  intensity  and  band  width 
parameters)  needed  in  the  computation  of  the  band 
absorptance  is  explained  in  Appendix  III. 

By  knowing  the  gas  density,  0,  optical 
lengch,  L,  and  the  above  mentioned  band  para¬ 
meters,  a  set  of  universal  band  absorption  cor¬ 
relations  can  be  found  for  the  computation  of  the 
total  absorpcance, A^ ,  of  the  ith  gas  radiation 
band.  Table  II  shows  the  correlations,  which 
relate  the  dimensionless  integrated  band  absoro- 
cance.A*  *  A/A0,  to  the  optical  depth  parameter, 
r  *  aoL/A0,  and  the  effective  pressure  parameter, 
1  »  3Pe. 

Table  II:  Band  Absorption  Correlations  (13) 


Table  III: 

Ranges  of  H2O  amd  CO?  Emission  Band 


Bands 

V1 

'Jc 

(u) 

(cm-1) 

(cm-1 

20 

0 

6.3 

1600 

2.7 

3760 

1.87 

5350 

1.38 

7250 

15 

667 

10.4 

960 

9.4 

1060 

4.3 

2.7 

3660 

2.0 

5200 

The  range  of  integration  for  each  band  is  shown 
in  Table  III.  For  the  case  of  a  nonoverlapped 
band,  if  the  lower  or  upper  bound  of  Che  band  is 
given,  the  other  bound  can  be  found  by 


If  Che  center  of  the  band,  is  given  (Table 
III) ,  both  upper  and  lower  bound  wave  number  can 
be  calculated  by 


v 

u,l 


c ,  I 


A./2 

1 


(23) 


l.i 


c.i 


A./2 

1 


(24) 


n 

<  1 


>  1 


min 

0 

n 

1/n 

0 


1/n 


A* 


1/7 

(Anx)  J"/ “  -  n 


In (nr)  +  2  -  n 


1  +  Inr 


where,  t  »  aoL/A0,  n  »  3P@,  and  A* 
chat  the  line  width  parameter,  3, 


A/A0.  Note 
and  the  band 


Appendix  A-III. 

Since  the  gas  emissivity  in  the  bands  is 
assumed  to  be  unity  while  outside  Che  band  is 
zero,  the  gas  emissivity  of  the  ith  band  may  be 
expressed  in  a  dimensionless  variable,  x. 


After  all,  the  gas  emissivity  of  ith  band 
whose  range  is  xj  ^  <  x 
as  follows. 


x,,  ,■  can  be  expressed 

a  »  X 


6  .  *^7  [S(3,0,x.  .)  -  S(3,0.x  .)’ 

g , 1  4  1,1  u » 1 


where , 


S ( 3 ,0 ,x)  »  7  ie"r 

n=L  ^ 


( :5  • 


-  x3~x  2(3+1) 
i-0  nX+1  r(3-i+l) 


The  gas  emissivitv  of  an  isothermal  gas 
volume  at  a  flame  temperature,  T,  and  pressure, 

P,  can  be  determined  by  summing  the  radiation  of 
all  bands  of  water  vapor  and  carbon  dioxide,  i.e. 


'g.i 


15 


I 

1 


3 

x 

band  i  eX-l 


dx 


(21) 


where,  x  -  hcv/KT. 

The  exact  solution  of  Eq .  (’ll)  can  be  found 
bv  using  a  similar  method  employed  for  Eq .  (13). 


£='■£.  (2b) 

g  -  g .  1 

The  emissivity  of  the  combined  radiation  is 
obtained  bv  introducing  Eqs.  (17)  ind  1 2b )  in  Eq . 
(1). 

Computation  of  Radiation  Heat  Transfer  in 
Diesel  Combustion  -  For  a  parametric  commutation 


4  ;  ■ 


9,  degree 

(A)  (B) 

Figure  1.  Normalized  Distribution  of  Burned  Fuel/Air  Ratio  and  Soot  Volume 
Fraction  in  (A)  Radial  Direction  and  (B)  Angular  Direction. 


of  radiation  heat  transfer  in  diesel  combustion, 
several  assumptions  were  taken  as  listed  below. 

(1)  The  system  of  which  size  is  small  but 
large  enough  co  represent  the  local  vol¬ 
ume  of  mixture  is  isolated  allowing  no 
heat  and  work  flow  across  the  boundaries. 

(2)  The  gaseous  species  in  the  system  attain 

an  equilibrium  state  of  the  adiabatic 
flame  temperature  determined  by  the 
local  burned  fuel/air  ratio  and  cylinder 
pressure.  ' 

(3)  The  radiation  source  is  a  2-dimensional 
surface  with  prescribed  opcical  length. 

(4)  The  fuel/air  ratio,  soot  volume  fraction, 
and  gas  species  concentration  are  avail¬ 
able  enabling  construction  of  the  2- 
dimensional  radiation  source. 

As  to  the  assumption  of  isolated  system,  the 
svstem  receives  the  piston  work  while  it  loses 
energy  chrough  radiation;  this  mutual  compensa¬ 
tion  may  partially  justify  the  assumption.  Al¬ 
though  the  processes  in  diesel  combustion  occur 
ac  highlv  non-equilibrium  states,  the  system  vol¬ 
ume  was  chosen  to  be  small  and  to  attain  an  equi¬ 
librium  for  simplicity  of  computation.  Since  the 
oDtical  length  is  varied  over  the  source,  the  2- 
dimensional  emission  source  was  considered  co 
represent  the  emission  characteristics  of  the  in- 
cvlinder  mixture  volume.  Note  a  studv  to  con¬ 
struct  a  3-dimensional  model  co  eliminate  this 
assumption  is  under  wav.  Although  in-cylinder 
species  information  will  be  obtained  either  ex¬ 
perimentally  or  theoretically,  the  present  compu¬ 
tation  was  carried  out  for  assumed  distributions 
of  species.  The  distributions  of  Che  burned  fuel, 
air  racio  and  sooc  concentration  assumed  tor  the 
present  parametric  studv  showed  chat  chev  de¬ 
creased  exponentially  along  the  radial  direction 
ind  held  a  normal  distribution  in  the  inguiar 
direction  .Fig.  1].  The  non-dimensional  form  of 
the  distribution.  : .  mav  be  expressed  is  follows. 


,  -ad  -br 

t  »  f  e  e  (. 

o 

where,  r,  J  are  Che  radius  and  angular  coordi¬ 
nates,  t0,  a,  and  b  are  to  be  determined  from 
experimental  data. 


L  Emitting  Volume 


Figure  2.  Two-Dimensional  Radiation  Emitting  .Source 

In  order  to  complete  tne  comput  it  ion ,  i  nu¬ 
meric. ll  icheme  was  emnl  'VMU.  Ficur*'  cnv, 

2  -  :  imens  un.i  1  -n  i  ;s  i  m  -ur:  uc  ;  i  v  i  : 
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sectors  according  to  the  fuel  plume  position. 

Each  sector  is  divided  into  n  rings  with  a  uni¬ 
form  width,  and  each  ring  is  again  cut  into  equi- 
size  cells.  By  setting  a  nodal  point  at  the 
center  of  each  cell,  one  obtains  n^  nodals  with 
the  same  area  in  a  sector.  In  each  nodal  point, 
the  fuel/air  ratio,  gas  species  concentration, 
soot  volume  fraction,  and  optical  length  were 
specified  according  to  Eq.  (27). 

Upon  the  implementation  of  the  above  numeri¬ 
cal  scheme,  a  computer  program  became  available 
for  computation  of  radiation  heat  flux  incident 
on  any  location  in  the  cylinder  head.  Table  IV 
shows  nineteen  subroutines  included  in  the  pro¬ 
gram  in  addition  to  the  main  program. 


Table  IV.  Subroutines  for  Various  Computations 
Subroutines  Subject  of  Computation 


AFT 

HEAT 

EQSTAT 

EQCON 

ROOT 

FACT 

RA 

SFACT 

POINTER 

PO INTER 1 

SOOT 

OPTI 

SERIES 

GAS 

CHIWHI 

BANDA 

A1 

AH20 

ACO 


Adiabatic  flame  temperature 
Enthalpies  of  species 
Equilibrium  stace  composition 
Gibbs  free  energy 

Positive  real  root  of  2nd  or  3rd  order 
algebraic  eq.  in  EQSTAT 
Gamma  function  values 
Combined  radiation  heat  flux 
Geometric  factor 
Interpolation  procedure  in  EOCON 
Interpolation  procedure  in  HEAT 
Soot  emissivity 

Optical  constants  of  soot  by  Eq .  (A15) 
Integration  of  Eq .  (16) 

Gas  emissivity 

V  and  i  defined  in  Eqs.  (A17)  &  (A20) 
Integrated  band  absorptance-Table  III 
Integrated  band  absorptance-Table  III 
Evaluate  the  emission  band  positions 
of  H-jO 

Evaluate  the  emission  band  positions 
of  COo 


PARAMETRIC  STUDY  OF  RADIATION  HEAT  TRANSFER 


The  construction  of  the  present  model  of 
radiation  heat  transfer  has  enabled  a  parametric 
analysis  of  its  processes  in  diesel  combustion. 
Several  assumptions  were  employed  in  the  analvsis. 

(1)  The  detection  surface  of  which  incident 
radiation  heat  flux  is  computed  is  paral¬ 
lel  with  the  emitting  surface  so  that  the 
shape  factor  can  be  explicitly  expressed. 

( 2)  There  are  four  plumes  of  pie  shape  with 
a  uniform  optical  length  of  10  mm  over 
the  emitting  surface. 

(3)  The  fuel  volume  fraction  of  stoichio¬ 
metric,  J,  and  soot  volume  fraction,  fv, 
have  profiles  in  the  plume  as  expressed 
in  the  following  equations. 


*  2.333  i  exp i —l.oi 
o 


exp  f -0 . 3r /R 1 


(23) 


f  »  2.333  f  ,  expl-1.5*-1/-  )']  expf-O.Sr  R)  i29) 
v  vO  o 

where.  i0  is  ~/4,  R  is  the  cylinder  radius,  and 


>0 >  fvQ  are  the  mean  values  of  the  fuel  volume 
fraction  of  stoichiometric  and  soot  volume  frac¬ 
tion  in  the  combustion  chamber,  respectively. 

Among  the  parameters  varied  in  the  computa¬ 
tion  were  the  position  of  detection  surface  on 
the  cylinder  head,  the  mean  burned  equivalence 
ratio,  the  mean  soot  volume  fraction,  and  the 
clearance  between  the  emitting  surface  and  the 
cylinder  head,  and  the  cylinder  pressure.  Some 
of  the  results  from  the  computation  are  shown  in 
the  following.  Note  chat,  for  comparison  pur¬ 
poses,  the  results  are  expressed  in  apparent 
emissivity  (this  will  be  represented  by  AE  here¬ 
after)  defined  in  Eq.  (5)  instead  of  the  heat 
flux. 

Location  of  Detector  -  The  radiation  heat 

flux  incident  on  different  locations  of  the  cyl¬ 
inder  head  was  computed.  The  AE  at  different 
locations  along  r-direction  and  over  9-direction 
are  shown  in  Figs.  3  and  4,  respectively.  The  AE 
monotonously  decreases  in  the  positive  r-direc¬ 
tion,  and  holds  a  bell-shaped  profile  over  the 
9-direction.  Although  expected,  it  is  interesting 
lo  find  such  non-uniformity  of  AE  over  the  cylin¬ 
der  head.  The  variation  is  much  greater  over  the 
9-direction  than  the  r-direction.  The  distribu¬ 
tion  of  AE  is  quite  similar  to  those  of  burned 
fuel/air  ratio,  species  concentration  projected 
on  Che  cylinder  head,  which  is  not  a  surprising 
result  since  the  radiation  heat  flux  is  strongly 
dependent  on  those  distributions.  This  suggests 
that  the  present  model  may  help  evaluate  the 
thermal  load  of  a  new  engine  head  once  the  engine 
fuel  injector  and  combustion  characteristics  are 
predicted  or  given.  The  bell-shaped  distribution 
of  AE  over  the  9-direction  indicates  chat  the  AE 
at  a  particular  location  is  almost  entirely  dic¬ 
tated  by  the  fuel  plume  right  underneath.  This 
is  mainly  due  to  the  small  clearance,  z,  between 
Che  cylinder  head  and  the  emitting  source  and  ics 
further  analysis  is  followed  in  the  next. 

Location  of  Emitcing  Plain  -  Since  the 

diesel  combustion  products  move  awav  from  the 
cylinder  head  as  che  piston  descends,  it  is  im¬ 
portant  to  evaluate  che  AE  for  different  clear¬ 
ances  between  che  cylinder  head  and  the  emitting 
plain,  z.  Figure  5  exhibits  some  interesting 
results.  For  smaller  z,  the  AE  of  locations  close 
to  che  axis  of  the  plume  is  greater  and  more  rap¬ 
idly  decreases  in  the  9-direction;  the  AE  there¬ 
fore  attains  greater  values  for  greater  z  at  the 
locations  far  from  che  axis.  This  is  explained 
bv  the  larger  value  of  the  shape  factor  for  great¬ 
er  z. 

Position  of  plumes  -  In  the  above  analy¬ 

ses,  it  was  found  chat  the  radiation  heat  flux 
incident  on  a  location  is  stronglv  affected  bv 
the  plume  right  below.  Figure  o  shows  che  evalu¬ 
ation  of  the  contribution  bv  the  sprav  plumes  on 
the  AE  at  locations  along  r/R  -  1.2  ind  z.R  - 
3  40  over  the  ‘-direction.  The  \E  due  to  the 
plume  right  under  the  detector  it  ■  -  1  is  t . ' 

while  the  AE  due  Co  ill  ocher  plumes  is  in  ;  v 
0.0026.  This  reaffirms  chat  the  AE  of  i  part  ic  i~ 
lar  location  is  almost  entirely  controlled  tv  the 
nearest  sprav  plume. 
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Figure  3.  Apparenc  Emissivity  ac  Different 
Locations  Along  the  r-direction 


Figure  4.  Apparent  Emissivity  at  Different 
Locations  over  the  3-direction 
with  z/R  «  3/40 
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Figure  3.  Apparent  Emissivity  due  to  the  Emitting 
Plain  Located  at  Different  Clearances 
with  r/R  =  1/2 


Figure  6.  Effect  of  Plume  Locations  m  \DDur>. 
Emissivitv 


Soot  Concentration  -  Since  the  representa¬ 

tive  in-cvlinder  soot  data  are  not  presently 
available,  the  effect  of  soot  volume  fraction, 
fv,  on  the  AE  was  analyzed  by  using  various 
values  of  fv  arbitrarily  chosen  but  within  5  to 
100  times  greater  than  the  reported  exhaust  fv. 
This  is  consistent  with  reports  by  Aovagi  et  al. 
(12)  who  suggested  that  the  amount  of  soot  per 
exhaust  gas  volume  is  approximately  1/40  of  the 
amount  of  soot  per  cylinder  volume  and  by 
Kittleson  and  Du  (20)  who  measured  to  find  that 
the  amount  of  soot  in  the  exhaust  is  approxi¬ 
mately  1/3  of  the  total  soot  appearing  in  the 
cylinder.  The  mutual  consistency  stands  since 
the  volume  of  the  soot  it  TDC  is  expected  Co  be 
jdd rox ima te ly  1.11  of  that  in  Che  exnaust.  Figure 
shows  che  AE  of  i  detecting  area  over  the  •- 
direction  with  r.  R  3  1,  2,  z/R  =  3/40.  The  AE 


increases  with  soot  concentration  is  expected. 
However,  it  is  interesting  to  find  that  the  emis 
sivitv  does  not  proportionally  increase  witn  the 
soot  concentration.  The  increment  of  che  AE  is 
larger  3C  lower  concentration  range  of  fv  and 
vice  versa.  The  partial  reason  for  che  conver¬ 
gence  of  AE  ac  high  soot  concentration  is  due  to 
the  exponential  function  in  the  relationship  of 
the  emissivity  and  the  absorption  coefficient  is 
shown  in  Eq .  (6)  . 

Soot  and  Gass  Radiation  -  Although  it  has 

often  been  mentioned  in  literature  '  3?d  that  the 
high  luminosity  of  soot  flames  in  diesel  combus¬ 
tion  is  responsible  far  the  diesel  radiation  lea 
loss,  onlv  i  limited  evaluation  has  neen  senar- 
itelv  made  on  che  contribution  on  AE  bv  loot  mu 
gas  of  eomDustion  nroducts  in  diesel  : /mouse  ion  . 


snows  cne  .u.s  v- 
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Figure  7.  Sooc  Volume  Fraction  Effect 
on  Apparent  Emissivitv 


Figure  3.  Apparent  Emissivities  bv  Soots  and 
Gases  in  Diesel  Combustion 
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Figure  Pressure  Etrect  on  Apparent  Emissivitv 
tor  a  High  Soot  Concentrat  ion  Gomoust ion 

computed  for  seemingly  i  typical  diesel  comous- 
tion  condition.  The  total  A£  which  results  from 
ooth  sooc  and  gas  is  found  almost  entirely  caused 
bv  the  soot  i.i  flames,  i.e.,  the  gas  radiation 
in  diesel  comouscion  is  extremelv  small  compared 
to  Che  soot  radiation.  This  mav  suggest  chat  the 
radiation  heac  loss  from  premixed  gaso line- fue led 
engine  combustion  is  negligible.  Note  that  che 
present  result  suggesting  the  gas  radiation  to  be 
a  mere  fractional  contribution  on  che  tocal  radi¬ 
ation  is  considerably  different  from  chose  com¬ 
puted  bv  Kunitomo  ec  al.  (10).  Their  results 
predicted  that  the  gas  radiacion  ranged  between 
20%  to  50%  of  che  combined  radiacion.  The  iis- 
crepancy  becween  the  present  results  and  their 
results  may  be  caused  by  che  differences  in  the 
computation . 

Pressure  Effect  — —  3ecause  of  che  inher¬ 
ently  high  pressure  reaction  condition  of  Jiesei 
combustion,  it  mav  be  important  Co  look  into  the 


Figure  10.  Pressure  Effect  >n  Apparent  Emissivit 
for  a  Low  boot  Concentration  combust:, 

pressure  effect  on  the  thermal  radiation  in  the 
combustion  champer.  The  inalysis  was  made  with¬ 
out  including  the  modification  of  fuel  sprav  ind 
atomization  affected  ov  pressure.  The  computa¬ 
tion  was  made  for  plumes  with  a  high  sooc  concen 
tracion  (  f„  -  ■4.0*10”'’)  and  with  low  concentra¬ 
tion  ( fv  =  4.0*10”')  as  shown  m  Figs.  1  and  10. 
respec c ive lv .  The  results  suggest  that  for  hign 
soot  concentration  flames  the  pressure  effect  on 
the  AE  is  almost  negligible.  Since  the  flame 
temperature  is  verv  insensitive  to  the  pressure 
variation,  the  results  indicace  no  et'ect  on  the 
sooc  radiation.  In  Fig.  10  plotting  ‘or  a  low- 
soot  radiation,  the  pressure  ei  feet  an  the  \i  .  , 
remarkable,  which  suggests  chat  there  exi.-,t  ,  t  v 
pressure  broadening  etfect  of  the  gas  riaiatun 
in  Jiesei  combustion.  Note,  however,  tint  ’•he 
impact  of  the  pressure  broadening  in  i  tvoi.  il 
lieseL  combustion  conuition  is  found  nee  i  ,  c  :  n  i  . 


jmui  nog  i  i  g  in  . 


RADIATION  HEAT  TRANSFER 


Radiation  Heat  Loss  in  Diesel  Combustion  — 
Due  to  the  lack  of  in-cvlinder  particulate  and 
gaseous  species  data,  the  presenc  parametric 
study  was  limited  for  an  engine  crank  angle,  al¬ 
though  the  new  model  can  be  employed  for  the 
evaluation  of  the  entire  period  of  combustion. 

A  computation  was  made,  however,  for  a  limited 
portion  of  the  combustion  period  by  using  a  set 
of  in-cylinder  daca  arbitrarily  chosen  buc  close 
co  real  based  on  the  following  assumptions.  (1) 
The  cylinder  pressure  increases  with  engine  crank 
angle  until  8°CA  after  top  dead  center  (TDC)  and 
decreases  thereafter.  (2)  The  soot  volume  frac¬ 
tion  increases  with  crank  angle  buc  starts  to 
decrease  16°CA  after  TDC.  The  relationship  of 
fv  to  CA  was  chosen  co  be  identical  co  that  of 
heat  release  rate  co  CA  based  on  the  experimental 
findings  reported  by  Kitrleson  and  Du  (20).  They 
found  chat  the  amount  of  in-cylinder  soot  almosc 
exactly  followed  the  heac  release  race  co  CA  his¬ 
tory.  (3)  The  burned  fuel/air  ratio  continuously 
increases  during  the  period  of  combustion.  (A) 
The  optical  length  is  assumed  to  be  10mm  over 
the  entire  emitting  surface  with  varying  species 
concentration.  The  combined  AE  (due  co  soot  and 
gases),  tc,  and  gas  AE ,  -tg,  computed  bv  using  the 
above  input  for  a  detecting  area  at  r/R  =  1/2  and 
1=0  are  shown  in  Fig.  11.  It  Is  interesting  to 
find  the  AE  is  very  low  prior  to  ISO^CA  while  fv 
is  relatively  high  and  vice  versa  after  190“2A. 
This  may  be  explained  in  terms  of  temperature 
chat  is  dictated  by  the  burned  fuel/'air  ratio. 
Again,  it  is  predicted  that  the  gas  AE  is  a  neg¬ 
ligible  portion  of  the  combined  AE  in  diesel  com¬ 
bustion  . 

Comparison  with  Other  Reports  -  The  pres¬ 

enc  model  will  not  be  compared  with  ocher  models 
(3,10)  since  thev  did  not  reporc  the  in-cvlinder 
species  daca  and  the  comoutacion  of  space-re¬ 
solved  heat  flux  in  diesel  combustion.  However, 
a  comparison  wich  some  experimental  results  is 
attempted  co  a  limiced  extent.  The  main  diffi¬ 
culty  encountered  in  the  selection  of  compatible 


experimental  data  for  the  comparison  was  me 
variety  in  measurement  techniques,  measurement 
locations  in  the  cylinder,  engine  operating  -:r.- 
dition,  types  of  engine  chosen  for  the  rrevicus 
data  acquisition,  etc.  ''2,36,33,39).  The  'ooten- 
tial"  radiation  heat  transfer  daca  react: ec 
Flynn  et  al.  (1972)  were  chosen  for  me  rre  =  e" 
discussion  since  their  data  did  not  ir.cl-:e  :-.t 
emissivity  of  the  detector.  The  same  :ata  w~ 
in  Fig.  11  were  used  co  compute  the  one.  si 
thickness,  KL,  co  compare  with  Flynn's  data . 
that  since  the  apparent  optical  thicxness  re¬ 
ported  by  them  were  compuced  oased  on  me  meas¬ 
ured  apparent  radiation  temperature,  an  ic 
ment  was  made  on  two  sets  of  data  of  di: recent 
injection  timing  co  achieve  a  more  consistin' 
comparison  wich  Che  present  results  Fie. 

Even  wich  the  above  mencionea  difficulties,  -  -»• 
interesting  notes  may  be  drawn  from  erne  - 
son.  The  prediction  is  lower  than  me  measure¬ 
ments  by  a  factor  of  2—*.  Among  proDaole  reason 
for  the  prediction  of  low  optical  chioxness  man 
be  the  substantially  greater  opcical  len.ctn. 
in  their  measurement  than  the  presenc  one.  How¬ 
ever,  ic  is  encouraging  to  find  chat  me  trend 
predicted  KL  relationsnip  to  CA  compares  well 
with  chat  from  the  experiment.  Vhen  the  in¬ 
cylinder  radiation  heac  flux  daca  become  lvail- 
able  from  the  present  ongoing  study,  che  mocel 
will  be  further  analyzed  for  the  verification. 

SUMMARY 

A  new  radiation  heac  transfer  model  or  die¬ 
sel  combustion  has  been  developed.  Use  of  the 
presenc  model  requires  in-cylinder  information  is 
co  time-  and  space-resolved  particulace  and  cas¬ 
eous  species  and  pressure-time  data. 

The  model  employs,  for  the  emissivity  of  i 
volume  gases  containing  sooc,  Ray ie igh- 1 imi t  ex¬ 
pression  to  compute  the  spectral  absorption  :oef- 
ficient  of  sooc,  the  dispersion  eauations  to  ccm- 
oute  che  refraction  indices  included  in  the 
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Rayleigh-limit  expression,  the  band  model  to  ob¬ 
tain  the  spectral  band  absorptance  of  gaseous 
mixture,  the  adiabatic  flame  temperature  obtained 
based  on  the  burned  fuel/air  ratio,  and  for  the 
computation  including  the  geometric  consideration 
of  the  combustion  system,  a  numerical  scheme  in¬ 
cluding  the  shape  factor  to  integrate  over  the 
entire  non-uniform  emitting  source. 

A  parametric  analysis  of  the  diesel  radia¬ 
tion  heat  transfer  was  carried  out  by  using  the 
present  model.  Among  the  findings  from  the  anal¬ 
ysis  are: 

(1)  The  thermal  radiation  flux  incident  on  a 
particular  location  of  the  cylinder  is 
mainly  originated  from  the  nearest  fuel 
plume  and  very  little  from  other  plumes; 

(2)  The  radiation  flux  greatly  decreases  as 
the  clearance  between  the  detecting  lo¬ 
cation  and  the  plain  of  emitting  source 
increases ; 

(3)  The  radiation  flux  is  generated  mostly 
due  to  the  presence  of  soot  and  negli¬ 
gibly  due  to  the  presence  of  gaseous 
combustion  products; 

(4)  The  soot  radiation  increases  with  the 
soot  concentration  greatly  when  the 
average  soot  concentration  is  low  and 
vice  versa; 

(5)  There  is  an  obvious  pressure  broadening 
effect  on  the  gas  radiation  in  diesel 
combustion.  However,  its  impact  on  the 
combined  radiation  in  diesel  combustion 
is  negligible; 

(6)  The  thermal  radiation  in  the  very  begin¬ 
ning  of  the  combustion  period  may  be  low 
due  co  the  low  temperature  even  though 
the  soot  concentration  is  potentially 
high  at  the  time;  and, 

(7)  A  prediction  by  che  present  model  qual¬ 
itatively  compares  with  some  of  the  re¬ 
ported  experimental  data. 
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APPENDIX  I 

Adiabatic  Flame  Temperature  Computation 

In  the  computation  of  adiabatic  flame  tem¬ 
perature,  eleven  gaseous  species  were  considered 
besides  the  fuel,  including  seven  dissociation 
reactions  with  known  individual  equilibrium  con¬ 
stants,  a  stoichiometric  equation  and  an  energy 
equation. 

The  stoichiometric  equation  for  the  reaction 
is 

C  H  +  w(I+J/4)  02  +  3.764v(I+J/4)  N, 

- -  (w-1) (I+J/4)  02  +  3 . 764w( I+J/4 )  N, 

+  I  CO,  +  J/2  H,0  (Al) 

In  addition,  the  equations  showing  che  dis¬ 
sociation  reactions  involving  eleven  species  are. 


C°2 

—  CO  +  1/2  0, 

(v 

(A2-1) 

H2° 
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—  H,  +  1/2  0, 

(e2) 

(A2-2 ) 

1/2 

H2 

+  1/2  0,  - - 

OH  (e3) 

(A2-3) 

1/2 

N2 

+  1/2  02  - - 

NO  (e ,  ) 

(A2-4) 

1/2 
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- -  H 

(e.) 

( A  2  —  5 ) 

4. 

1/2 

o2 

- -  0 

fV 

( A2-6 ) 

1/2 

N, 

-  N 

(e_3 
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(A2-7) 

w  : 

•  (F/A)  / ( F/A) 

=  l/> , 

(A3) 

F/.A  is  che  molar  fuel-air  racio,  (F/A)s  is  che 
stoichiometric  molar  fuel-air  racio,  J  is,  there¬ 
fore,  the  fuel  volume  fraction  of  stoichiometric, 
and  e^  is  the  extent  of  the  corresponding  reaction 
<A2-i)  which  is  defined  as  the  ratio  of  the  change 
of  each  species  to  its  stoichiometric  coefficient 
(18)  . 

In  the  computation,  seven  extents  of  reaction 
i  -  1,  ....  7)  were  used  instead  of  the  elev¬ 
en  species  mass-terms  in  che  equations  by  identi¬ 
fying  the  number  of  moles  of  each  species,  namelv, 
ZO->,  ZN-i,  ....  etc.,  and  che  total  number  of  moles 
in  the  product,  Za^. 

Z,  “  l)  <  I*  I  /  4 )  +  ( e ,  +e  ,  -e  ,-e  .  -e  ,  1  2  <  A —  —  1 3 
a  1  2  3  ->  6 


3 .  76uv (  I  +■  i  .  4  > 


e_  1 


A-.-21 


VW.V-V 
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ZOH  ■  e3 


So  ’  e4 
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(A4-11) 


Z  ,  =  4.764w(I+J/4)  +  J/4 

ail 


+  (e1+e2+e5+e6+e7) /2 


Seven  equilibrium  constant  equations  for  Eqs. 
(A2)  can  be  written  in  those  terms  defined  in  Eqs. 
(A4)  as  follows 


ZC0  V 


(P/Zall> 


Z  Z  1/2 
H2  °2 


Kp3  =  1/2  1/2 

h2  °2 


„  OH 

P4  *  1/2  1/2 
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where,  P  is  the  cylinder  pressure. 

Substituting  Eqs.  (A4)  and  (A5)  into  the 
above  equilibrium  constant  equations,  one  obtains 
seven  simultaneous  algebraic  equations.  In  addi¬ 
tion,  the  energy  equation  of  the  adiabatic  system 
is  written  as  follows. 

I  Z  [ <H  -H  )  +  (AH  )  ]  * 

prod.  1  2  f  f  1 

>  Z  f(Hx  -H  ,  +  (AHf)  ]  (A7) 

reac .  J  If  f  J 

where  is  enthalpy  at  initial  temperature,  Hf-, 
is  the  enthalpy  at  flame  temperature,  AHf  is  the" 
heat  of  formation  at  reference  temperature  Tf, 
and  Tf  is  273°K. 

A  new  iteration  scheme  was  developed  to  find 
the  adiabatic  flame  temperature  and  the  equilib¬ 
rium  composition  by  using  the  above  equations  and 
reported  Kp's  and  iHf's. 

The  key  of  the  method  is  the  simplified  forms 
of  equation  after  introducing  Eqs.  (A4)  into  Eqs. 
(A6) .  They  are 

0. 5(P-Kp2)e2  +  (C1-C2+IKp2)e2 

+  (2IC,-0.5I2Kp2)e1-I2C2  =  0  <Ad) 

^Kp2el+(I_el)Kpl'e2  *  0-3[Je1Kp2 
+  (e3+e5> (I-e1)Kp1 ]  «  0  'A9> 

(l-l/4KpZ)e2  +  [C3/2-e5/4+e,/2]Kp2e, 

+  (ej/2-e7)C2Kp2  =0  A 1 j 

( 1-Kp2 / 4 )  e2  +  0. 5(C.+C_Kp2'ie  -  :  *  J  \ 

U  -*  434-*  -4  J 

(ZC0,Kpl)e5  +  °‘  5'2iZ0,Kp5e3 


eiVe2-°-5VKp5  =  ° 


[  (I-e^JKp^+e^Kp^/Zle^-  hw-lMI*: 

+  0.  Sfe^+e^-e^-e^  ]  )e1<pfj  *  ) 

(zco2KpI)e7  +  °-52co2o2Kp7e7  *  2‘:o2), 


-  3.764w(I+J/4) ] Kp “  =  0 
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-  [ (w-1) (I+J/4)  +  (a^+e^-e^-eg) /2 ] 

C4  -  [(w-1) (I+J/4)  +  (e1+e2-e3-e6)/2]Kp^ 

and 

C5  -  3. 764w( I+J/4)  -  e  / 2. 


[ 1-exp (-  l  ±xk5k1  -CT) 

a(T)  =  aQ  - ^ -  (A16) 

[l-exp(-  T  ±x  i  ]  V(Tq) 
k=l 

where , 


A  brief  inspection  of  the  abpve  equations 

suggests  that  if  the  values  of  e2»  e3 . sj 

are  fixed,  the  value  of  can  be  easily  found  as 
an  exact  solution  of  a  simple  third  order  alge¬ 
braic  equation.  Similarly,  others  can  be  found 
in  sequence.  This  is  repeated  to  find  a  set  of 
converged  values  of  e^,  which  are  used  to  find 
che  equilibrium  composition  in  Eqs.  (A2) .  The 
above  process  is  iterated  to  find  a  convergence 
of  che  adiabacic  flame  cemperacure  of  the  compo¬ 
sitions  . 


APPENDIX  II 

Optical  Constants  of  Soot 
3y  Using  Dispersion  Equations 


The  complex  refraction  indices  of  che  soot 
are  calculated  by  using  che  dispersions  equations. 
The  dispersion  equations  chat  are  known  to  include 
che  contributions  of  both  free  and  bound  electrons 
are  shown  as 
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where,  Fj  is  the  number  of  effective  electrons 
per  unit  volume;  hj  is  the  electron  damping  con¬ 
stant;  u  is  the  frequency  of  the  radiation; 
is  Che  natural  frequency  of  the  j th  electron;  e 
is  the  electron  charge,  1.602*10'^C,  m  is  che 
mass  of  the  electron,  9 .  l*10“-^Kg ;  and  Eg  is  che 
electric  permittivity,  8. 854*10-l-2c2sec2/Kg-m3. 
Note  chat  the  numerical  value  of  each  constant  is 
shown  in  Table  AI . 


APPENDIX  III 

Computation  of  Band  Parameters  of  Cases 

The  integrated  band  intensity  of  each  band 
is  approximated  by  (34,35), 


x,  =  hcv,  /KT 

k  k 

X0,k  =  hcvk/KTo 

T-  =  100°K 
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In  che  above  equations,  the  statistical 
weight  gk  is  unity  for  the  degenerate  vibration 
and  the  vg  .  is  zero  if  a  plus  sign  is  associated 

with  Sk  in’ the  upper  state  (vk+o^ . vk+; k , . . . ) 

and  is  5k  when  a  minus  sign  appears.  The  quan¬ 
tity  1q  had  to  be  determined  experimentally  (30, 
31,33)  for  each  species  of  the  combustion  prod¬ 
ucts. 

The  band  width  parameter,  Ag,  is  approximated 
by  (13), 

Ao  =  Aoo(T/To)1/2  (A18) 


where,  Agg  is  experimentally  determined  for  each 
band . 

The  line  width  parameter,  3,  is  calculated 
by  (35), 


1/2  XT) 
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(A19) 


where , 
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Two  weak  bands  of  HoO  require  a  special 
method  of  computation  3S  explained  in  the  follow¬ 
ing. 

For  the  inclusion  of  two  weak  bands  of  H  -.0 , 
the  following  method  was  taken  in  the  present  com¬ 
putation.  In  che  computation  of  i  ind  ,  for  the 
pure  rotacional  band  (  "k  =0,  K  -  1,2,3),  the 
following  equations  were  used. 


tlT)  =  i 
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and 

B(T)  -  30(T/T0)'1/2  (A22) 

For  Che  2.7um  H20  band  where  chree  overlap¬ 
ping  bands  occur,  a  was  found  by  summing  a-Cj  » 
1,2,3)  by  using  Eqs.  (A16)  and  (A17).  Values  of 
3j  were  compuced  by  using  Eqs.  (A19)  and  (A20) 
and  averaged  by 


3 


3 


(a. 3.) 
J  J 


(A23) 


Table  AIII  shows  the  recommended  values  of 
3q,  3q,  and  Aqq  for  che  major  absorption  bands  of 
HjO  and  C02  (13). 


Table  AI. 

The  Constants  of 

the  Dispersion  Equations  (14) 

Electron 

Fj  (electrons/m^) 

(rad . /sec . ) 

(rad. /sec. ) 

free  c 

4.06  x  1027 

6.00  x  1015 

— 

bound  1 

27 

2.69  x  10 

6.00  x  1015 

1.25  x  1015 

bound  2 

2.86  x  1026 

7.25  x  1015 

7.25  x  1015 
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ABSTRACT 


A  computational  study  of  radiation  heat  transfer  in  a  combustor  having 
fuel  spray  combustion  is  reported  in  the  present  paper.  When  the  distribution 
of  radiatively  participating  species  is  known  in  either  a  normal  distribution 
or  an  exponential  distribution  in  any  direction  of  a  cylindrical  coordinate 
with  respect  to  a  referenced  location,  e.g.,  in  a  form,  F  *  f  exp  (-ap-b02-cz2) , 
where  f,  a,  b,  and  c  are  distribution  constants  for  individual  species,  the 
governing  equation  of  radiation  heat  transfer  without  scattering  in  such  a 
system  can  readily  be  solved  by  using  the  herein  reported  method.  Several 
results  from  the  present  authors’  work  are  explained  that  are  employed  for 
implementing  the  solution:  a  new  coordinate  transformat  ion  method  to  fine  the 
species  distribution  along  individual  optical  paths  centered  at  any  chosen 
location  in  the  combustor;  a  new  formula  for  the  adiabatic  flame  temperature  in 
a  logarithmic  function  expressed  in  terms  of  relevant  variables;  a  new  integral 
function  facilitating  solutions  for  various  radiation  equations;  a  new  inverse 
error  function  enabling  some  exact  solution  of  the  governing  equation  of 
radiation  heat  transfer,  etc. 

An  exact  solution  is  found  for  both  optically  thin  and  thick  media.  And  a 
numerical  solution  is  offered  for  optically  medially  thin  medium.  The  present 
analytical  method  is  applied  to  radiation  heat  transfer  in  a  direct  injection- 
type  diesel  engine  combustion  chamber.  Some  representat i ve  results  are 
discussed  that  have  been  obtained  from  a  parametric  analysis  of  its  radiation 


heat  transfer. 


INTRODUCTION 


When  an  equation  of  radiative  heat  transfer  is  solved  for  predicting  the 
heat  flux  incident  upon  various  locations  over  the  surface  of  a  combustor 
having  flame  plumes,  several  difficult  problems  are  encountered.  In  order  to 
briefly  review  the  difficulties,  a  combustor  having  a  jet  flame  is  considered 
as  shown  in  Fig.  1.  The  heat  flux  through  an  optical  path  (r,  9,t)  centered  at 
location,  D  may  be  calculated  by  solving  the  governing  equation  of  radiation 
neglecting  scattering  effect, 


qeo 


I  1.(0)  cose  di 
0  A 


when  the  spectral  radiation  intensity,  I  (r)  is  represented  as 
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r  <^(r ' )  r’ 
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r 

where  <x  is  the  spectral  volume  apsorptance  at  (r.e.t)  and  efeA  is  the  spectral 
blackbody  radiation.  The  total  radiation  flux  on  location  D,  then,  is 
calculated  by  integrating  over  the  entire  hemispherical  volume  faced  by  the 
location.  Among  the  main  difficulties  in  implementing  the  above  computation 
are  that  can  be  determined  only  when  the  distributions  of  radiatively 
participating  species  and  temperature  are  known  along  the  individual  optical 
paths  and  that  I.  is  spectrum-dependent. 

In  order  to  obtain  a  better  solution  of  Eq,  (1),  some  new  methods  were 
proposed  for  overcoming  the  abovementioned  difficulties  when  the  authors 
computed  radiation  heat  transfer  in  a  direct  injection-type  diesel  engine:  a 
coordinate  transformation  technique  to  find  the  species  distribution  along 
individual  optical  paths  when  the  distribution  is  given  with  respect  to  a 
chosen  location,  e.g.,  the  injection  nozzle  hole  [1];  a  computation  method  of 


finding  the  spectral  volume  absorptances  of  soot  and  gas  for  their  mixtures 
[1,2]  by  using  reported  results  [3,^.5]  and  new  analytical  and  numerical 
methods  of  solving  Eq.  (1),  which  is  considered  in  the  present  paper.  This 
paper,  therefore,  may  serve  as  a  summary  report  of  the  authors'  modeling  study 
of  radiation  heat  transfer  in  a  combustor  having  flame  jets. 

Elaborating  the  species  distribution  in  the  combustion  chamber,  it  was 
found  to  conveniently  be  described  by  the  following  equation  in  a  cylindrical 
coordinate  with  respect  to  the  injection  nozzle  [1]  when  distribution  is  <nown 
in  either  a  normal  (or  a  skewed  normal)  distribution  or  an  exponential 
distribution  in  any  direction  of  a  cylindrical  coordinate,  e.g.,  in  a  form; 

F  *  f  exp  (-ap-b<?J-cz2)  (3) 

where  F  represents  either  fuel/air  ratio  or  soot  concentrat ion  or  or  H_3 

<;  2 

concentration;  f,  a,  b  and  c  are  constants  to  be  separately  determined  for  eacn 
medium  by  either  experimental  or  theoretical  means.  The  expression  is 
considered  to  reasonably  well  describe  the  species  distribution  of  non- 
axisynmetric  plumes  in  view  of  some  reported  results  from  both  theoretical  and 
experimental  studies  [6,7],  and  the  new  expression  offers  various  advantages, 
e.g.,  that  the  expression  can  be  used  to  describe  the  spray  plume  in  swirl 
motions  by  including  its  effect  in  q-direction  and  that  it  facilitates  the 
sought  solution  of  Eq.  (1)  [1],  When  the  species  distribution  with  rescect  to 
referenced  position  in  the  combustor  is  known,  the  species  distribution  along 
the  optical  path,  r,  of  any  chosen  direction  (r,6,^)  and  location,  D 
(p.,'>.  ,z.)  may  be  found  as  [1], 

Qua 

F  -  f  exp(-[(r-r  )/r  ]2  (u) 

o  d  o  w 

wh4re  f  ,  r  and  r  are  functionally  related  to  9,  f, ,  ojt  ,  z  .  f,  a,  b,  and 
u  o  w  1  q  Q  Q 

c.  This  distribution  equation,  then,  may  be  used  for  determining  the  volume 
absorptance  according  to  the  reported  methods  [3,p,5]  as 


2 

<  -  exp  (-c, (r-r.)  )  (5). 

o  1  r  1  i 

The  temperature  distribution,  T,  needed  in  the  above  equation  and  e  for  Eq. 

D  ^ 

(1)  may  be  determined  by  using  the  distribution  of  burned  fuei/air  ratio,  w, 
from  Eq.  (3)-  This  can  be  achieved  by  using  our  new  formula  for  the  adiabatic 
flame  temperature  [3]; 

T  -  T  *  T,  2,nu  ♦  T_  (Slnu)2  (6) 

o  1  d 

where  T1  ,  and  T^  are  functionally  related  to  the  reaction  pressure,  the 
initial  mixture  temperature  and  the  number  of  carbon  atoms  of  the  fuel 
molecule.  Even  with  the  above  rearrangement  and  simplif ication  of  terms  in  Eq. 
(1),  its  solution  was  not  readily  found  until  some  new  mathematical  techniques 
were  employed  as  explained  in  the  following. 


ANALYTICAL  AND  NUMERICAL  SOLUTIONS 

The  solution  of  Eq.  (1)  may  be  considered  for  the  following  three  cases: 
(1)  the  optically  thin  medium;  (2)  optically  thick  medium  and  (3)  optically 
medially  thick  medium:  the  solution  for  each  case  will  be  explained  upon 
further  simplifying  the  governing  equation  as  follows.  Introducing  new  terms, 

,  kT  r 

<  *  <,.A  [3.9],  x  »  hc/kXT  and  a  -  —  /  <  dr, 

o  a  he  .  o 


where  h  is  Planck's  constant;  c  the  speed  of  light  and  *  the  3oltzmann 
constant,  Eq.  (1)  can  be  rewritten  for  directional  radiation  intensity  as, 


15  ok  ,  o. 

~5 — —  J  <  T 

m  he  0  o 


®  ^  -ax 

f  — - -  dx  dr. 

0  ex-1 


During  the  process  of  simplifying  the  above  equation,  t.ne  authors  discovered  i 


new  useful  integration  function  [10],  i.e.. 
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x-  -ax  m 

f  - - —  dx  -  Ftm.a.x, )  -  F(m,a,x  ), 

x  ,  1  2 

x,  e  ~i 


F(m,a,x)  *  I  {e 
n-1 


-(n+a)x  j.  j-  ) —  (n+a)  *  1  x^  1  ]  > ,  where  r(m*l)  -  ■  ,  2, 

i-0  r(m~i+1) 


Mote  that  Eq.  (7)  can  also  be  used  for  determining  a  solution  in  closed 
form  of  the  black  body  radiation  function,  Fq_^.j.  C11], 

XT  S  (?) 

Fo-x-  *  j.  ^r~d(XT) 

0  o  i 


4  *  c£2!tx3  .  aaf  *  fit .  &^,3  . 

it  ,  n  n  2  _3 


From  Eq,  (7),  since  one  can  find 

®  u  -ax  « 

/  -  dx  =  £  - ■?  , 

0  eX-1  n»l  ( n+a ) 

the  governing  equation  is  simplified  as 

™  r  < 

360  ok  .  1  f  o  o 

105  '  "5  hc  n-1  n5  0  (1*a/n)5 

r 

With  the  definition  of  the  optical  depth,  t  -  /  dr,  the  equation  may  be 

0 

rewritten  in  a  tidier  form, 

t  360  ak  "  1  ,To  T5 

l  .  r  »  - 5 — : -  l  — =■  J  - - - c  d  T  '  d  J 

71  nc  n-1  n5  0  d*^l)5 

nhc 

The  solution  of  this  equation  is  considered  for  the  previously  mentioned  three 
cases  as  follows. 


OPTICALLY  THIN  MEDIUM:  When  *  approaches  zero,  i.e.,  an  optically  thin  limit, 
an  exact  solution  for  Eq.  (S)  was  found.  In  this  case,  Eq.  '3;  becomes 


,  ~  36C  ak  _  1  ,‘o  „5,,  5kTr, 

xee:  *  ^h^~  -."si  T  vl"  dT 

n- 1  n  0 


"or  finding  the  solution  of  the  aoove  equation,  since  ?  is  functionally  related 

to  the  optical  depth,  t  and  the  burned  fuel/air  ratio,  ,,  and  since  nay  be 

r  2 

found  as  u  ■  u1 e  from  Eq.  (4),  Eq .  (6)  may  be  written  as 

7  »  ?  ♦  T1  (In  ’jl)1  -  r2)  +  7 2  (In  -  r2)2  -  7^  -  r2  *  7^  r  (1C) 

2 

where,  T,  -  T  *  T,  In  u,  *  T_(  2-n  u,  )  and 
3  o  1  12  1 

Tu  -  T1  +  2T2  Zn(u^ ) . 

In  addition,  for  obtaining  the  solution  of  Eq.  (3),.  it  was  necessary  to  find 
the  relationship  between  t  and  r.  From  definition  and  Eq.  (5),  the  optical 

depth,  t  is  rewritten  as, 

,  ,2 
/  ~c  ^ ( r-r 1 ) 

t  «  <  f  e  dr 

0 

<1  /—  .  _  _ 

-  —  /  —  {erf  L'C,  ( i — r , )  ]  ♦  erf  (/c.  r , ) } 

2C,  11  11 


/Q 

or,  erf  [/c~  (r-r,)]  -  —  /—  -  erf  (/c~  r , ) . 

11  <  IT  11 

In  view  of  Eq.  (10)  and  for  further  simplifying  Eq .  (11),  a  variable 
transformation  was  employed  to  obtain 


(r-r,). 


t  *  —  /— -  [t  -  rr-  S'-?-  erf  (/c~  r,)3 
'  2c,  II 


fin 


<.  IT 


*  d,  Ct-t,] 

Chen  Eq.  (11)  may  be  written  as 
erf  (r)  •  t 


or , 


r  -  erf  ( r ) 


-1 


(  1  2) 


7  ) 


where,  erf  represents  the  inverse  error  function,  whose  exact  solution  is 


found  as  (note  the  details  of  its  derivation  are  explained  in  Appendix 


A,  -  m  -3  7  tt  2  — 5  , 

—  (T  *  —  T  +  ~  T 


‘text,  the  above  equation  is  introduced  into  Eq.  (10)  to  find 
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T  “2  *  -  “U 
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where , 


4  ?4  and  ?6 


16  ‘2' 


From  Eq.  (15),  one  also  may  find  the  following  in  order  to  use  then  i 

5  5  4-2  4-4 

7  '  Y  -  5  ~3  T5  7  *  5  T3  T6  ’  • 

T6  -  7  6  -  5  7  5  r2  *  6  7  5  7,  7\ 

3  3  3  « 

Finally,  the  above  equations  and  Eq .  (12)  are  introduced  into  Eq 


obtain 
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Consequently,  Eq.  (16)  can  be  used  for  finding  an  exact  solution  of  the 
equation  of  radiation  for  an  optically  thin  medium  in  a  combustor  with 


plumes  as, 


l  .  360 _ok  !  i_  ((T  5  .  T  6  d  T  _  [d  2(t  _t  }2 

3 C  tt  he  5  3  nhc  3  1  0  2nnc  1  0  T1 

n» l  n 


S  14 
-  (-  T  7 
3  3  5 


10  k  r 


15kT3T3r.  »,  ,«  »  «, 

2nhc  "1  T0  T1  di  T1  1 


(T„V - -  -3— 6 )  (d.5(x  -t.)5  -  d.5  x.5) 


v-  ir»  i.-v  v*  VA  V  L.'  ^ 


Optically  thick  medium:  An  exact  solution  can  also  be  found  in  an  optically 

thick  medium:  assuming  the  volume  absorptance  <  is  0  for  r  <  r  ,  and  <Q  for  rQ 

<  r  <  r  *c ,  one  can  rewrite  Eq.  (2)  as 
o 

r  +e  r" 

I,_  •  f  °  —  e.  exp  (-  /  <  d r")dr’ 

o  •  b  0 

This  equation  is  further  solved  to  find  the  result, 

r  *e  <  r' 

I,.  »  /  0  —  e,  exD  (-  /  <  dr")dr' 

9£  it  b  o 

r  r 

o  o 

<  e  .  r  +e 

-  -  /  exp(-<  (r’-r  ))dr* 

it  oo 

r 

o 


oT  r.  -7, 
— —  [  1  — e  J 

TT 


(13) 


Optically  medially  thin  medium:  When  the  system  holds  an  optically 
medially  thin  (or  thick)  medium,  Eq.  (2)  cannot  analytically  be  solved  as  the 
above  discussed  cases:  rather,  a  numerical  method  was  employed  for  it. 

Without  going  into  the  details  of  its  derivation  (shown  elsewhere  [11]),  but  by 
knowing  that  the  volume  absorptance  <(r)  can  be  found  from  Eq .  (b),  the 
following  result  was  used  for  the  solution. 
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15  o 


at 


T  .  J 
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I  /  I  <  .  .  i  [F(4,a,x .  , 

*i  t.  ,  j-i  01J  J_1 

1-1  J 


)  -  F(U,a,x  ) ]d7 


(19) 


where,  <Q  is  the  volume  absorptance  defined  in  the  ith  subrange  of  r-domain 
and  the  j th  subrange  of  i-domain.  The  numerical  results  of  Eq.  (19)  were 


obtained  using  Caussian  integration  method. 


RESULTS  AND  CONCLUSION 


The  present  computational  method  enables  an  extensive  parametric  study  of 

radiation  heat  transfer  in  a  combustor  with  flame  jets  [2].  Since  this  summary 

paper  is  to  primarily  report  the  methodology  of  computation  by  using  the 

authors'  several  new  findings  reported  elsewhere  [1  ,  2,  3,  10],  only  some 

representati ve  results  considered  for  a  direct  inject  ion-type  diesel  engine  are 

presented  here.  For  this,  in  order  to  be  consistent  and  to  have  better  mutual 

comparisons,  the  same  computational  conditions  employed  for  the  previous 

results  are  considered  here.  The  numerical  values  of  f  in  Eq.  (3)  were,  for 

soot,  HgO,  CCj  and  fuel/air  ratio,  3  x  10  °,  0.01,  0.01  and  1.0,  respectively; 

the  distribution  constants,  a,  b,  and  c  were,  for  those  species,  0.6,  2. 36  and 

2.C,  accordingly.  Note  again  that  those  constants  and  distribution  made  in 

each  coordinate  of  Eq.  (3)  can  be  determined  by  computational  or  experimental 

method,  or  arbitrarily  chosen  for  parametric  studies.  The  engine  details  for 

the  computation  were:  the  chamber  surface  temperature,  500°K;  the  number  of 

spray  plumes,  the  piston  radius,  R  *  4.92  cm;  the  combustion  chamber  bowl 

radius,  3-44  cm;  the  compression  ratio,  23;  the  engine  speed,  1000  rpm;  the 

surface  spectral  emissivity,  0.95;  the  fuel  composition,  ;  the  time  of 

lo  3 u 

computation  in  cycle,  40  degrees  after  top  dead  center.  The  results  are  shown 

U 

in  normalized  forms:  normalized  radiation  heat  transfer,  Q/jT  ,  where  T  is 

2400°K  and  o  is  the  Stefan-Boltzmann  constant;  normalized  directional 

4  a 

intensity,  I  /irc7  ;  normalized  spectral  radiation  heat  transfer,  Q^/oT  . 

The  radiation  heat  transfer  through  a  solid  angle  incident  on  the  cylinder 

head,  are  shown  in  Fig.  2.  Where  there  are  four  flame  jets  in  the  combustion 

chamber,  radiation  heat  transfer  flux  incident  on  the  locations  along  the 

radial  direction  varies  depending  upon  the  direction  in  the  hemispheri cal 

volume  of  integration.  The  schematic  drawing  in  the  figure  shows  the  spherical 


coordinate  identifying  the  direction  of  solid  angles.  Note  that,  for  cases  9  * 

60°  and  76. 7°,  heat  transfer  on  r/R  <  0.1  is  in  large  from  plume  It  1  and  that 

oeyond  r/R  -  0.1  from  both  plumes  #1  and  2.  The  directional  radiation 

intensity,  Ia,  was  combined  with  the  spectral  emissivity  of  the  chamber 

ot, 

wall,  e,,  »  0.95,  to  obtain  radiation  heat  transfer  through  the  cylinder  head  by 

integrating  as,  Q  =  /  I  cos0  sin  0  d8  d£  (Fig.  3).  The  net  heat  transfer 

at, 

computed  with  inclusion  of  wall  emission  at  500°K,  reflection  from  the 

opposite-side  walls  are  represented  by  "Net";  those  obtained  excluding  the 

reflection  is  denoted  by  "No  reflection”;  and  the  heat  flux  incident  on  the 

surface  without  including  wall  effect  are  shown  by  "Spray  plume".  Explaining 

the  results,  the  rapid  reduction  in  heat  transfer  around  r/R  *  0.7  is  caused  by 

the  presence  of  the  piston  bowl.  The  double  hump  in  the  plot  may  be 

surprising.  In  order  to  clarify  this,  an  additional  computation  made  for  solid 

angle  of  zineith  angle,  9  -  39. 7°,  at  varied  azimuthal  angles  and  locations  on 

the  cylinder  head,  is  3hown  in  Fig.  U.  From  the  results,  one  finds  that  while 

the  radiation  heat  flux  incident  on  locations  around  the  nozzle  (r/R  »  0)  is 

from  all  of  the  four  plumes,  the  flux  on  those  near  r/R  »  0.1  is  from  one 

plume,  i.e.,  the  nearest  plume  to  the  location.  If  t.ne  comcutation  is  made  for 

a  single  plume,  the  heat  transfer  should  have  asymptotically  been  approaching 

zero  at  r/R  -  0.  This  trend  is  qualitatively  compared  with  the  experimental 

results  obtained  by  using  a  single  plume  [12]. 

Since  the  distribution  constants  being  used  for  Eq.  (3)  in  the  present 

analysis  were  chosen  based  on  an  in-cylinder  soot  measurement  [13],  it  seems  of 

interest  to  compute  radiative  heat  transfer  contributed  by  in-cylinder  soot  and 

gaseous  species,  i.e.,  mainly  C02  and  H^O.  Fig.  5  shows  their  spectrum- 

resolved  radiation  heat  transfer  smoothed  over  individual  wave  bonds  m 

U 

normalized  form,  Q^/oT  ,  with  respect  to  normalized  wave  number,  5  jm/^.  It 
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appears  clear  that  thermal  radiation  transmission  caused  by  the  presence  of 
soot  almost  dominates  the  entire  radiation  process  and  the  radiation  .by  gaseous 
species  is  very  small.  The  great  portion  of  radiation  by  gas  found  around 
6um/A  is  explained  by  the  strong  emission  bands  of  H^O  at  A  =  2.7  and  CO^  at 
A  -  2.7  and  U.3  um.  Since  the  spectral  absorptance  of  surface  can  be  ' 
incorporated  with  results  like  those  in  Fig.  5,  the  analysis  may  be  made  for 
combustors  having  various  surface  coatings  with  known  spectral  emissivity. 

The  main  issue  associated  with  the  use  of  the  present  analysis  method  of 
radiation  heat  transfer  may  be  justification  of  the  species  distribution 
described  by  Eq.  (3).  Since  the  distribution  widely  varies  depending  upon  the 
fuel  injector  and  combustor  condition,  it  is  difficult  to  exactly  describe  the 
variations  by  a  single  equation.  Reported  results  suggest,  however,  the 
individual  species  distribution  may  be  simplified  by  using  exponential  or 
normal  (or  skewed  normal)  distribution  in  each  coordinate  component  of  a 
cylindrical  system,  e.g.,  Eq.  (3).  It  is  found  in  the  present  study  that  the 
radiation  heat  analysis  is  greatly  simplified  when  the  distribution  is  made  ir. 
sucn  a  form  without  resorting  to  the  rather  complex  zonal  method  or  the  use  of 
the  geometric  factors.  Further,  the  use  of  the  present  method  readily  enables 
parametric  analysis  of  radiation  heat  transfer,  e.g.,  in  order  to  identify  a 
desirable  injector  suitable  for  a  chosen  combustor  system. 
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AM  INVERSE  ERROR  FUNCTION 


The  error  function,  erf(x),  although  defined  in  several  ways,  is 
essentially  equivalent  to  the  following  expression: 


erf(x)  »  —  /  e 
/rr  0 


X  -t2 


dt 


—  z  (-nn  x 


2n+1 


/t?  n-0 


n ! ( 2n+1  ) 


(A-l  ) 


The  equation  is  the  integral  of  the  so-called  Gaussian  or  normal  function  and 
occurs  frequently  in  the  study  of  the  general  theory  of  probability  [1U],  the 
analysis  of  transient  heat  flow  in  a  semi-infinite  solid  [15],  and  the 
computation  of  radiation  heat  transfer  with  radiatively  participating  media 
[10],  etc.  For  meeting  the  needs  of  such  studies,  the  error  function,  its 
derivatives  and  integrals  have  been  tabulated  [16].  As  found  in  the  text,  it 
was  needed  to  obtain  the  inverse  form  of  the  error  function,  i.e.,  in  Eq.  (A- 
1),  the  value  of  x  for  a  corresponding  erf(x).  The  following  derivation  shows 
a  new  inverse  error  function  in  closed  form: 

When  the  error  function,  erf(x),  is  set  equal  to  y,  i.e.,  y  =  erf(x),  its 


inverse  form  may  be  expressed,  for  the  convenience  of  discussion, 


x  *  erf  (y) 


(  A  — 2 ) 


or, 


x  -  fre(y).  ( a-3 ) 

By  Taylor's  expansion  method,  the  inverse  error  function  can  be  expressed 


as 


,  :  *<n,<0)  n 

x  *  L  -  Y 

n-0  n! 


(  A-u ) 


(n) 


where,  x  (0)  is  the  nth  derivative  of  x  with  respect  to  y  at  y 


.  n 


order  to  find  x 


\  » w 


(0),  the  following  steps  are  taken. 


x.(y)  -  **  «  (^f1 

x  dy  dx; 


1%-L—  /  e't2  dt]}'1 

dx  A  o 


/t t  X 

r  • 


x"(y)  -  Cx’(y)]  -  ^  [x'(y)][g]  1 


u  2  x 
■  2  X  6 


x(3)(y)  -  ^  [x ' '  (y )  ] 


2  /— 

r  /-  i,  2x  2x  i  r  *  ^  X  ^ 
-  [^  0  *  )  e  JL— ^  e  -1 


*3/2  il  2i  3x2 

r  (1  +  ^x  )  eJ  , 


x(U)  (y  )  -  ^ —  ( 7  x  *  1 2x3)  eU*  , 


5/2  2 

x(5)(y)  *  |  (7  *  92x2  *  96xS  e5x 


3  2 
x(6)(y)  .  ( i 27x  *  652x3  ♦  U80x5)  e6x 


7 /2 

x(7)(y)  ,  (127  ♦  3^80x2  *  1022«xJ 


Consequently,  one  can  rewrite  Eq.  ( H)  as 


+  5760x°)  e 


" 


A  close  investigation  of  the  above  results  indicates  that  a  term  L^-j  is 
repeated  in  each  derivative  and  that  the  even  derivatives  become  zero  at  y  • 
This  leads  to  writing  Eq.  (5),  the  inverse  error  function,  in  a  tidier  form, 

fre(y)  -  I  fn(0) y)2n+1 /(2n>l ) !  |y|  <  1  <A-€ 

n»0 

where,  f  (0)  is  obtained  from  the  reciprocally  defined  function,  f  (t), 
n  n 

expressed  as 

f  ( t )  ■  1 ,  and 
o 

f  , (t)  -  f "  (t)  *  2(Un  ♦  3)  x  f ’ (t) 
n+l  n  n 

+  2[2n-1  *  4(2n+1 ) (n+1 )x2]  f  (t).  (n  -  0,1,2,.. 

n 

Among  the  analytical  properties  of  the  new  inverse  error  function  are 
fre(0)  -  0 

fre(-y)  «  -fre(y),  and 
fre( 1 )  »  *. 

In  addition,  the  convergence  of  Eq.  (6)  is  demonstrated  in  Figure  A-1. 
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Nomenclature 


r 

a  Species  distribution  constant  or  (kT/hc)  /  <Q  dr 

o 

b  Species  distribution  constant 

c  Species  distribution  constant  or  speed  of  light 

f  Species  distribution  constant 

h  Plank's  constant 

I  Radiation  Intensity 

k  Boltzman  constant 

q  Radiation  heat  flux 

r  Optical  path 

R  Radius  of  cylinder 

T  Temperature 

p,  4> ,  z  Components  in  cylindrical  coordinate 

r,  5,  £.  Components  in  spherical  coordinate 

J  Volume  absorptance 

X  Wavelength 

t  Optical  depth 

a  Stefan-3oltzman  constant 

Subscripts 
a  Adiabatic 

b  Blackbody 

d  Detector 


X 


Spectral 
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COORDINATE  TRANSFORMATION  METHOO  FOR 
RAOIATION  HEAT  TRANSFER  PREDICTION  IN 
SOOT  IA0EN  COMBUSTION  PRODUCTS 

by 

S.  L.  Chang  and  K.  T.  Rhee 
Department  of  Mechanical  and  Aerospace  Engineering 
Rutgers,  The  State  University  of  New  Jersey 
New  Brunswick,  New  Jersey  08903 

INTRODUCTION 

Computational  modeling  of  radiative  heat  transfer  has  recently  become  a 
more  important  means  for  the  analysis  of  thermal  loading  of  combustion  chambers 
with  luminous  flames.  This  seems  to  arise  due  to  both  the  renewea  recognition 
of  need  for  an  extensive  radiation  heat  transfer  analysis  to  achieve  a  better 
design  of  many  practical  systems  (e.g.,  furnace,  turbine  combustor  and  diesel 
engines)  and  promising  new  contributions  from  recent  studies  enabling  more 
comprehensive  modeling.  In  the  modeling,  however,  several  difficult  problems 
are  faced,  mainly  due  to  the  directional ity  and  spectral  nature  of  radiative 
processes.  Most  serious  problems  to  be  solved  in  improving  its  modeling  may  ce 
1 isted  as  fol lows: 

(1)  A  viable  computation  of  spectral  volume  absorptance  of  a  mixture  of 
combustion  products; 

(2)  A  more  accurate  computational  method  for  integrating  the  eyuat’on 
radiation  heat  transfer  along  individual  optical  paths  in  the  hemi spher < ca : 
volume  faced  by  each  location  of  the  reactor  wall  (hereafter  called  a  detector) 
and  over  the  wave  length  range  of  thermal  radiation;  and 

(3)  A  detailed  description  of  in-reactor  distributions  of  combustion 
products  and  temperature  in  terms  of  a  suitable  coordinate  system  centered  at 
an  individual  detector  and  its  coupling  with  the  above  radiative  transfer 


eouati on . 


In  the  present  paper,  being  mainly  concerned  with  the  last  item  of  the 
list,  a  brief  discussion  is  given  for  others  in  order  to  review  the  overall 


2. 


■problem. 

spectral  volume  absorptance  of  soot/gas  mixture 

When  the  spectral  volume  absorptance  of  combustion  products  are  calculated 

for  the  range  of  thermal  radiation,  the  respective  absorptances  due  to  soot 

cloud  and  gas  mixture  are  computed  as  well  as  their  summation.  The  most  widely 

accepted  method  of  computing  the  absorptance  of  soot  r  ,  by  using  either 

\  •  s 

Rayleigh  limit  expression  [1]  for  median  size  of  particles  smaller  than 
0.1  micrometer,  or  Mie  theory  [2]  for  particles  with  greater  sizes  where  it  is 
required  to  know  the  size  distribution  to  take  into  account  scattering.  Since 
practical  combustion  systems,  in  general,  produce  small-size  soots,  the  Rayleigh 
expression  is  only  considered  here.  The  expression  is  independent  of  the  size 
of  the  particles,  inversely  proportional  to  the  wavelength  l,  and  propor t i ona 1 
to  the  soot  volume  fraction,  i.e., 

36  n2  k(*A)  f 

_  _  v _  ,  i  , 

rX~?  ?  2  <1  ?  (  *■  /' 

[n  (  1-k  )  +  2]  +  4  nV 

The  optical  constants,  n,  k,  the  real  and  imaginary  parts  of  the  refractive 
indices,  were  measured  from  the  intensity  of  polarizing  light  reflected  from  the 
soot  particles  and  the  measurements  were  fitted  to  the  dispersion  equation  [3,4] 
It  was  reported  that  *-  may  slightly  increase  with  wavelength  and 
insignificantly  varies  with  kind  of  soot  and  temperature.  For  computation  of 
emissions  from  gaseous  mixtures,  the  empirical  wide  band  model  proposed  by 
Edwards  and  Bal arkri shnan  [5]  has  often  been  used  [6-8].  The  significance  of 
the  method  is  to  find,  instead  of  finding  gray  body  absorptance  of  gas  volumes. 


the  width  of  individual  emission  bands,  A  ,  compared  to  blackbody  spectrum  to 

approximate  the  strength  of  original  emission  spectrum.  This  unique  technique 

facilitates  an  evaluation  of  the  spectral  dependence  of  gas  radiation,  when 

emission  bands  overlap  each  other  at  a  single  spectrum,  a  separate  correction  is 

given  [5].  With  regard  to  the  surrmatlon  of  both  emissions  due  to  soots  and 

gases,  even  though  any  method  has  not  fully  been  verified,  several  techniques 

have  been  used  in  the  past:  for  an  assumed  homogeneous  and  gray  body  emission, 

by  adding  the  gas  emittance  times  the  soot  transmittance  to  soot  emittance  [9]; 

by  adding  the  soot  and  gas  emittances  es  and  to  find  total  emittance,  e ^ , 

i.e.,  et  =  1- ( 1— e  )(l-c$)  [10];  by  introducing  a  new  approximation  using  a 

pentagarmia  functional  form  [11],  etc.  Unlike  those  for  computation  of  gray  body 

emittance  of  gas/soot  mixture  volume,  the  surrmation  of  spectral  emissivity  of 

such  mixtures  needs  special  consideration  in  order  to  evaluate  the  spectral 

dependence  of  radiation  heat  transfer.  A  new  method  of  summing  spectral  volume 

> 

absorptances  has  been  used  by  the  authors  to  find  the  spectral  emittance  of 
combustion  products  comprising  both  soot  and  gases  for  their  recent  studies  [121. 
In  the  computation  of  the  spectral  volume  absorptance  of  a  soot/gas  mixture, 
the  spectral  volume  absorptance  of  gaseous  mixture,  is  expressed,  for 

convenience  of  computation,  in  a  form 

Vg  =  CA-  :21 

Here,  C  was  found  from  the  relationship.  A,  =  l-Exp(-Cl/\)d\  upon  the  computa 
tlon  of  A^  for  each  emission  band  by  using  the  Edwards'  wide  band  non-gray  model. 
On  t'he  other  hand,  the  volume  absorptance  due  to  soot  suspending  in  the  mixture, 
r\  s’  ^  1s  use<^  to  C1  from  the  expression  s  =  C'/\.  The  spectral 

absorptance  of  soot/gas  mixture,  «- ^ ,  is  then  computed  from  the  following: 


^  =  (OC  )  A 


This  method  for  quantifying  spectral  dependence  of  radiative  processes  or 
a  combustion  gas  mixture  with  soot  suspension  is  basically  similar  to  one  of 
those  listed  above  for  gray  body  radiation  analysis.  That  is,  a  surmiation  for 
individual  spectral  emissivity  of  such  mixtures  by  using  the  new  method,  without 
exceeding  blackbody  emissivity  at  a  spectrum  upon  the  addition  of  both  emissions 
is  quite  equivalent  to  the  following  relationship: 

£t.\  =  (4) 

The  validity  of  the  new  summation  method,  however  similar  to  one  often  employed 
for  gray  body  analysis,  is  not  fully  accepted  at  present,  so  that  a  future 
systematic  experimental  evaluation  of  summation  is  needed.  Further,  since  the 
calculation  of  emittance'  as  a  thermal  equilibrium  mixture  is  made  by  computing 
the  spectral  volume  absorptance  instead  of  emittance  as  explained  above,  the 
assumption  of  Klrchhoff's  law  to  hold  for  it  in  a  system  is  more  justified  than 
a  similar  assumption  made  for  gray  body  computations. 


EQUATION  OF  RA0IAT10N  HEAT  TRANSFER  EQUATION 

The  equation  of  radiation  heat  transfer  in  radiatively  parti ci pati ng  fluids 
confined  in  a  volume  is  well  known.  The  governing  equation  of  monochromatic 
radiation  heat  transfer  Is  sought  for  the  processes  in  combustors  containing 
plumes  of  small  size  soot  laden  flames  as  shown  in  Fig.  1  where  pertinent 
details  are  shown,  e.g.,  geometric  information  and  spectral  reflectivity  of  wall 
surfaces,  8^,  etc.  Referring  to  Fig.  1,  the  governing  equation  of  local 
spectral  radiation  intensity,  I  (r),  is  written  as 

A. 


r  o  \(r'  ) 
I.  (r)  =  - - 


(  r  '  )  Exp( -/  r  i r " ) dr " ) dr '  , 


where,  e  =  — = - — --- -  ,  *■  is  the  spectral  volume  absorptance  at  (r, 9,5), 

bX  \5  [Exp<hc/\KT)-lj 

c  *  2.998  x  1010  cm/sec,  h  =  6.625  x  lO-27  erg-sec,  and  k  =  1.380  x  10‘16  erg/K. 

The  heat  flux  along  an  optical  path  of  a  solid  angle  (8,?)  in  the  thermal  radiation 
range  may  be,  then,  computed  from 


q 


8? 


1(0)  cos  8  d\ 


(6) 


The  mafn  difficulty  in  finding  the  solution  of  Ed.  (5)  resides  in  the  fact  tha 
r  is  a  function  of  both  species  concentrati ons  and  temperature  and  that  both 
variables  are  a  function  of  the  location  along  the  individual  optical  path,  r.  'he 
present  paper  considers  a  new  technique  to  describe  the  species  distribution 
along  individual  optical  paths  and  its  coupling  with  £q.  (5).  When  such 
distribution  details  become  available,  a  more  realistic  solution  of  Eq.  (6)  is 
accessible,  as  reported  elsewhere  [13]. 

The  surface  properties  are  combined  with  the  radiation  heat  fluxes  incident 
on  individual  locations  of  the  combustor  wall.  The  extinction  coefficient, 
t  =  JY  -dL,  along  each  optical  path  found  from  the  computation  as  explained 
above  and  the  spectral  reflectivities,  s  ,,  are  included  in  the  expression 

A.  *  1 

for  the  radiosity  of  surface  1,  3., 


-T  , 


(Ij  ♦  Bxl  cos  8 |  I x ( 0 ) )  +  8vl  cos  8 |  e  (I2  +  cos  I v ( 0 ) ) 


BXl  8X2  cos  91  C0S  92 


Upon  imp! ementa ti on  of  computations  by  coupling  Eq.  (5)  with  Eq.  (7),  the  net  he 
transferred  through  a  chosen  spectrum  at  a  detector  can  be  calculated. 


IN-C0M8UST0R  SPECIES  DISTRIBUTION 


Since  an  analysis  of  space-  and  spectrum-resolved  radiation  heat  transfer 
In  a  combustor,  as  considered  here,  requires  detailed  distributions  of  in-reactor 
species  and  temperature,  search  of  a  useful  technique  for  this  is  clearly 
warranted.  Unfortunately ,  at  present  neither  convenient  means  for  obtaining 
in-combustor  data  sufficient  for  use  in  radiation  heat  transfer  analysis  has 
been  available.  The  distributions  may  be  obtained  by  either  experimental  or 
theoretical  method.  Some  activity  is  underway  in  the  direction  of  obtaining 
such  Information  in  the  authors'  laboratory.  When  the  distribution  becomes 
available,  it  Is  most  convenient  to  express  It  in  a  coordinate  system  centered 
at  the  fuel  injection  nozzle. 

Although,  as  explained  above,  the  species  distribution  in  plumes  of 
combustors  equipped  with  the  fuel  injection  nozzle  is  not  readily  available,  a 
new  approach  explained  In  the  following  may  shed  light  in  the  search  for  a 
viable  solution  to  the  problem.  It  is  found  from  some  experimental  and 
theoretical  results  in  literature  that  the  species  distribution  in 
non-ax i  syirmetric  plume  of  such  systems  [14,15]  may  conveniently  be  described  by 
the  following  equation  In  a  cylindrical  coordinate: 

F  *  Fq  exp(-ao-b«^-cz^ )  (8) 

where  F  may  represent  either  fuel/air  ratio,  or  soot  concentration,  or  CO^  or 
H^O  concentrations;  Fq,  a,  b  and  c  are  constants  to  be  separately  determined  by 
either  experimental  or  theoretical  means  for  individual  combustion  units.  This 
expression  of  plumes,  however,  deemed  rather  simple  and  probably  less  realistic, 
is  considered  to  serve  a  suitable  method  for  some  parametric  analysis  of  the 
confronting  problem.  Among  the  advantages  of  this  approach  is  that  the  deformation 
of  plume,  e.g.  in  spiral  direction,  due  to  swirl  motion. while  the  fuel  spray  is 
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formed,  may  be  reasonably  described  for  both  transient  and  steady  combustion 
systems,  i.e,  by  including  it  In  $-di rectlon,  as  explained  later. 

The  next  task  to  be  implemented  for  the  final  solution  may  be  a  proper 
coupling  of,  with  Eq.(5),  the  details  of  in-combustor  species  distribution  that 
may  be  obtained  from  the  above  method  or  others.  Several  techniques  have  been 
employed  In  finding  a  solution  of  the  equation  of  radiation  heat  transfer  (Eq. 
(5),  namely,  the  Monte-Carlo  method  [10,16],  the  zonal  method  [17,18],  the  use 
of  geometric  factors  [6,19],  etc.  Those  methods  did  not  rigorously  implement 
the  solution  of  Eq.  (5),  but  sought  for  approximated  solution  and  gray-body 
computation.  Among  the  main  reasons  for  this  is  that  it  is  difficult  to  find 
species  distribution  along  individual  directions  of  integration  in  the 
hemispherical  volume  faced  by  the  location  where  heat  transfer  is  considered. 
The  difficulty  is  compounded  when  the  distribution  varies  with  time  as  in 
transient-flow  combustors. 

COORDINATE  TRANSFORMATION 

In  view  of  the  discussion  above  on  the  overall  problem  in  finding  space- 
and  spectral -resol ved  radiation  heat  transfer  in  combustors,  it  is  highly 
desirable  to  find  a  more  versatile  method  enabling  us  to  find  a  better  solution 
of  Eq.  (5).  A  new  technique  for  the  goal  is  presented  in  the  following.  The 
basic  idea  of  the  technique  Is  to  employ  coordinate  transformations  for  finding 
the  species  distribution  along  the  individual  optical  direction,  r,  (Fig.  1)  in 
an  equational  form  by  using  the  species  distribution  given  in  a  coordinate  with 
respect  to  the  fuel  injector.  When  such  an  equation  in  a  suitable  coordinate 
centered  at  detectors  on  the  combustor  wall  becomes  available,  i.e.  along 
individual  optical  direction  with  respect  to  the  detector,  a  better  computation 
of  the  space-resol ved  radiation  heat  transfer  can  be  attempted. 


As  explained  earlier,  it  was  found  to  be  convenient  to  express,  in  a 
cylindrical  coordinate,  the  spray  plumes  of  which  geometry  and  others  vary  with 
time.  The  following  discussion  is  concerned  with  description  of  the  same  plumes 
in  a  new  coordinate  system  with  respect  to  a  detector  of  interest,  consequently 
to  make  available  the  species  distribution  along  individual  optical  paths 
centered  at  the  detector.  This  goal  is  achieved  by  using  a  new  coordinate 
transformati on  technique  Involving  two  main  sets  of  coordinate  systems,  i.e.,  a 
cylindrical  coordinate  centered  at  the  injection  nozzle  hole,  0,  for  describing 
the  spray  plumes  (<>,<&, z)  and  a  spherical  coordinate  (r, a,c)  for  expressing  the 
same  plumes  with  respect  to  the  detector  located  at  D.  Two  transformation 
techniques  are  presented:  (1)  point-to-point  transformati on ;  and  (2)  distribution 
equation-to-distributlon  equation  transformation.  The  former  is  considered  to 
be  particularly  useful  when  some  discrete  point  data  are  available  in  a  coordinate 
with  respect  to  the  Injection  nozzle  for  expressing  them  in  a  new  coordinate 
centered  at  a  detector.  The  latter  is  developed  for  transforming  a  functional 
expression  of  the  spray  plume  in  a  nozzle  coordinate,  e.g.  Eq.  (8),  to  find  an 
equational  spray  description  in  a  new  detector  coordinate. 

( 1 )  Point-to-point  Transformati ons 

MUTUALLY  PARALLEL  COORDINATES  -  When  in-combustor  species  data  is 
experimentally  acquired  [20],  It  is  always  likely  that  they  become  available 
in  a  set  of  discrete  data  points  as  opposed  to  in  a  form  of  continuous  variation 
In  species  concentration  of  plumes.  When  one  uses  such  discrete  data  points  in 
computing  radiation  characteri sties  In  a  direction  of  optical  path  with  respect 
to  a  detector,  It  is  necessary  to  know  a  species  distribution  along  individual 
optical  paths  for  solving  Eq.  (5).  For  such  needs,  a  new  approach  of  using 
point-to-point  coordinate  transformat i on  is  introduced.  To  perform  a  coordinate 
transformation  between  these  systems,  two  more  intermediate  coordinate  systems. 


X  and  Y,  are  introduced  for  the  convenience  of  derivation.  The  axes  of  the 
coordinates  under  consideration  are  either  parallel  or  perpendicul ar  to  each 
other  in  order  to  avoid  the  complexity  of  the  derivation  although  cases  with 
the  mutually  non-parallel  coordinates  are  discussed  later.  The  interim 
coordinate  systems  are  listed  in  the  following  and  shown  in  Fig.  2- ( A ) . 


a. 

Spray  plume  coordinate  (cylindrical), 

C 

(  0 

b. 

X  coordinate, 

X 

(xi’x2,x3 

c. 

Y  coordinate. 

Y 

(yl,y2’y3 

d. 

Detection  coordinate  (spherical) 

S 

(r,e,0 

Identifying  the  location  of  a  detector  at  X|d<  x^.  a  series  of 

transformations  are  carried  out  to  obtain 


2  2  2  2 
o  =  r  sin  9  +  2  r  sine  cos  (<?d+0  +  , 

®  =-tan~^[(od  sin®d  -  r  sine  si nf<)/(od  cos  ®d  +  r  sine  cose)], 


Z  = 

Z  .  -  r  cose 
d 

1 

(9) 

2 

2  2 

•d  =  tan_1( 

-1 

where. 

°d 

=  X1  +  x2  , 

x  x  ^ )  > 

and  zd  =  Xy 

Note  that  ® 

=  tan  *(x7/x- 

w  1 

if  X, 

*  0 

r — * 

1 

c 

4-> 

II 

© 

( x2/xl )  +  * 

if  x1  < 

0  and  X2>  0 

and  ®  =  tan 

‘  x 2/x i '  "  '  i# 

<  0  and  ^  <  0. 

MUTUALLY  UNPARALLEL  COORDINATES  -  However,  since  in  many  practical  systems, 
e.g.,  diesel  engines,  the  fuel  is  injected  into  the  combustion  chamber  at  an  ang 
9 .  ,  as  shown  in  Figs.  1  and  2— ( 8 ) ,  this  has  to  be  taken  into  consideration  in 
the  coordination  transformati on .  In  order  to  meet  this  need,  an  additional 
coordinate  system  having  an  angle  of  e.^  with  respect  to  the  cylinder 
coordinate  (x^.x^.x-j)  is  introduced.  The  new  coordinate  (u, ,u9,u,) 
is  not  parallel  to  the  detector  coordinate  (y^,y-,,y^).  The'r  mutual 
+el at i onshi p  may  be  listed  as  follows: 


*  o  cos*,  %2  -  0  Sin*,  X3  -  L  . 

U1  =  X1  cosdin  '  x3  Sin  9 in  ; 
u2  =  x2  ’ 

u3  =  X1  Sin9in  +  x3  COS91n  ; 


2  2  2 
=  U1  +  u2  *» 


8  =  tan  (u-/u. )  ;  and 

u  v  2  1 ' 


Z  -  u_ 
u  3 


Further,  an  expression  similar  to  £q.  (9)  can  be  obtained  as 

2  .222 

0  -  to  cose,  cos*  -  2  sin  9.  )  +0  sin  * 

u  In  in 

*  =  tan_i[(o  sin*)/(p  cose,  cos*  -  z  sine.  )] 

u  in  i  n  J 

z  =  0  sine,  cos*  +  z  cose .  . 

u  in  in 

Accordingly,  the  location  of  detector,  initially  identified  by  using  the 

parallel  coordinate  (p,e,z)  can  finally  be  rewritten  in  terms  of  the  unparallel 

system  (o^.r^zj  (see  Fig.  2-(8)  as, 

2  2  2  2 
o,,  =  (pH  cose.  cos*.  -  z,  sine.  )  +  p,  sin  0  , 

u  d  in  a  d  in'  d  d 

$  =  tan  ^[(o,  sin(j>.)/(0.  cose.  cos<t>,  -  z.  sine.  )]  and 

u  L  d  d  d  in  d  d  in/J 

z  =  0  .  sine.  cos<t>.  +  2.  cose.  .  (10) 

u  d  i  n  d  d  i  n 


Since  it  is  readily  found  that  the  geometric  configurations  of  the  inclined 
plume  vs.  the  coordinate  (y^y^.y^)  is  equivalent  to  those  of  the  similar 
plume  with  axis  in  x^  vs.  an  equally  inclined  detector  coordinate 
(uj.u^u^),  some  pertinent  conversion  may  be  made.  Converting  the 
geometric  details  of  a  given  solid  angle  initially  represented  by  using  :a,q) 
into  an  expression  in  terms  of  correspond! ng  unparalle'  coordinate  ( 0  }, 

the  results  obtained  for  the  case  of  mutually  parallel  coordinate,  Eq.  ;?>,  c. 
directly  be  used.  For  this,  the  conversion  for  the  sol'd  anqie  is  acnieveo  o- 


(11) 


9  =  cos”  f-sine  cos^  sin  9.  +  cose  cose.  1 

u  1  in  l n J 

c  =  tan'^fsine  sinc/(sine  coss  cose.  +  cose  sine.  )] 

The  location  of  a  data  point  originally  identified  in  terms  of  the  coordinate 
plume  of  which  axis  is  inclined  by  an  angle  of  e.^  with  respect  to  a 
perpendi cul ar  coordinate  can,  therefore,  be  expressed  in  terms  of  a  spherical 
coordinate  centered  at  a  detector  on  the  combustor  wall  by  introducing  Eqs.  ( 
and  (10)  i nto  Eq .  ( 11 ) . 

( 2 )  01 stri bution-to-Qi stribution  Trans  forma tl ons 

WITHOUT  SWIRL  MOTION  -  As  explained  earlier,  the  details  of  in-combustor 
spray  plume  may  be  obtained  in  an  equational  form,  e.g.  Eq.  (9),  expressed  w; 
respect  to  the  nozzle  hole.  The  transformati on  of  such  distribution  to 
individual  distributions  along  an  individual  optical  path  at  a  detector  is 
considered  for  a  combustor  without  swirl  motion.  .The  logical  step  for  this  goa 
would  have  been  introducing  the  correspondi ng  terms  of  Eq.  (9)  into  Eq.  (3). 
Such  transformati on ,  however,  results  in  an  extremely  complex  form  of  Ec.  (31 
expressed  in  terms  of  the  spherical  coordinate  at  the  detector.  Furthermore, 
the  new  transformed  distribution  equation  becomes  too  difficult  to  be  usee  in 
Eq.  (5).  In  order  to  alleviate  the  problem,  an  approx imati on  technique  was 
sought  in  such  a  way  as  to  conveniently  use  its  results  in  obtaining  the 
solution  of  Eq.  (5).  It  was  found  that  use  of  the  Taylor  expansion  metnoo  'o 
the  transformati on  not  only  greatly  facilities  the  steps  toward  tne  f'naT  gca 
but  also  lead  to  very  accurate  results. 


o- terms : 


Letting  \  *  r  sins/pd,  A  *  cos(«d  +  ?),  with  f  *  o/od 

?  1/2 

f(x)  =  (1  +  2Ax  +  x  )  , 

f‘(x)  =  (A+x)/f  ,  and 
■  fM(x)  =  [f2(x)  -  ( A+x ) 2 ]/f3  ■ 

For  x  =  x^  +  c  and  e<<1  the  above  may  be  rewritten  as, 

fix)  =  f(x1)  +  f’U,)  £  +  f"(xi)E2/2 

Further,  since  e  =  x-x.,.  fix)  m ay  be  expressed  as 

f(x)  s  o/od  =  +  fl  +  f2  x2 

where,  fQ  =  f(Xl)  -  ^ '  ( x1 )  Xj  +  ^"(x^  x2/2, 

=  f ' (x1 )  -  f"(x1 )  x1 ,  and 
f 2  s  ^"(Xi )/2- 

For  X>>1,  fi»)  -  o/od  =  x  +  A 

<t»  terms: 

From  En.  (9),  having  f(x)  =  tan(«) 

=  (sin*  .  -  x  sin<;)/(cos<t>i  + 


one  finds 


(12) 


x  cose  )  , 


one  obtains 


2 

f'(x)  =  -sin{<j>d  +  ?)/(cos«d  +  x  cose)  » 
f"(x)  =  -2  f'(x)  cose/(cos«j>d  +  X  cose), 

*(x)  =  tan“1(f) 

♦'(x)  =  ^'/(1+f2),  and 

♦"(x)  =  [-2  ff'2  +  (1+f2)  f"]/(l+f2)2. 

2 

Note,  for  cos<s>d  +  x  cose  =  0,  <t>  ( x )  =  */2,  <t> '  ( x )  =  -cos  e/sin(<t>d+e)  and 

<t>"(x)  =  0.  For  x  =  xi+s  and  e<<l,  $(x)  may  be  written  as  <s(x)  =  + 

<t> '  ( X  ^ )  e  +  <>',(xi  )  e2/2. 

2 

In  view  of  E q.  (8),  one  defines  <t(x)  -  <9  (x)  to  rewrite 

♦  (X)  =  *  <  X 1 )  +  *'(xi)  E  +  <9"(xi)  e2/2 

where  ♦'  (x)  =  2<J» ( xn- )  * '  ( X1 ) 

*"(x)  =  2*(Xl)  ♦  ,,(x1)  +  20'2(Xi) 

Further,  since  c  =  x-Xii  ♦  may  be  expressed 

♦(x)  *  ♦  +  ♦.j  X  +  X2  (13) 

where,  4>g  =  ♦{x^  -  ♦  ’  ( )  xi  *  x2/2 

*2  s  ♦(x^  -  ( Xi )  xi 

*2  *  ♦,,(x1)/2 

-  -  2 

For  x>>1»  *(x)  =  -e,  or  <9(x)  =  e  • 

z-terms : 

Since  z  Is  a  single  first  order  polynomial  function  of  r,  the  following  is 


Obtained  from  Eq.  (9) . 


2  :  2  .  ■  2  2,  r  COS  8  +  r  COS  8 

d  d 


Finally,  by  substituting  terms  in  Eq.  (8)  with  Eqs.  (12),  (13)  and  (14),  one 
obtains 

F  =  f  Q-Kr-ro)/rJ2  (1 

Q 


where,  f  .  =  f  e"(arbl  /4Cl\ 


bx/2  Cr 


a  °o  +  b  *0  +  C  zd  » 


=  -(a  +  b  sin8/od  +  2  C2d  coss,  and 

=  [a  p.,  +  b  <>2]  sin28/od2  +  C-cos2e. 


When  the  above  results  are  used  in  solving  Eq.  (5),  it  is  found  that  the  above 
approximation  method  generates  better  accuracy  when  the  integral  interval  ;s 
divided  into  smaller  subranges. 


WITH  SWIRL  MOTION  -  The  present  coordinate  transformation  method  is  also 
considered  for  a  combustor  with  swirl  motion.  With  swirl  motion,  prior  to  the 
transformation,  the  spray  plume  distribution  equation,  Eq.  (8),  needs 
adjustments  to  take  into  account  the  bending  of  spray  plume.  Although  the  exact 
formation  mechanism  of  bent  sprays  in  the  presence  of  gas  swirl  in  the 
combustor  is  not  properly  understood  at  present,  the  following  two  probable 


cases  are  considered. 


t  (■  p «» ■«*i  » i « i » l ' rr^TrTTrT*Tr' 


Droplet  with  Constant  Axial  Velocity: 

It  is  assumed  that  (1)  fuel  droplets  ejected  from  the  nozzle  hole  to  move  at 
a  constant  axial  velocity,  v^.;  (2)  the  angular  velocity  of  the  swirl,  Q,  is 
constant  along  the  axis  of  the  plume;  (3)  the  angular  velocity  of  fuel  droplets 
along  the  plume  axis  is  the  same  as  Q.  With  the  above,  a  point  character i sti c 
at  a  location  (o,<s,Z)  of  no-swirl  condition  is  found  at  a  location  (p,o+®  ,Z), 
due  to  the  swirl  during  a  period  of  time  at,  according  to 
4>s  =  (ft/vf)p, 

•  =  fl/vf  and 


7  7 

In  this  case  Eq.  (8)  becomes  F  =  Fq  exp[-ao  -  b(d»-n>  )  -Cz  ].  The 
transformation  of  the  above  equation  into  the  detector  coordinate  expression 
Is  exactly  identical  except  for  a  new  set  of  equations  in  d>- terms.  They  are: 

♦(x^  -  *2(x1 ), 

4  '  (x1 )  =  2'i(xi  )  u» '  ( x  1 )  *  and 

♦"(x1 )  =  2 -4» ( x i )'!»"( x i )  +  (16) 


where,  *(xf)  =  a»(  x  i )  -  ♦jfx,), 

’>>'  (xi )  =  o'  (xi )  -  4>s  (  x  i )  and 
o"(xi )  =  o"(xi )  -  os(xi )  • 


Droplets  with  Specified  Axial  Velocity: 

Melton  [21]  analyzed  a  case  for  the  plume  formation  as 


v,  =  v  /( 1+ao ) 


to 


where, v$  is  the  exit  velocity  of  the  fuel  droplet  at  the  nozzle  hole,  a  is 
about  0.085/Rs  ahd  Rs  is  the  radius  of  the  nozzle  hole.  Assumptions  of  the 

previous  case  are  also  similarly  taken  here  for  finding  =  Q- 1 ,  since 

2 

v^  =dp/dt  =  v  / ( 1+ap )  and  p  +  ap  /2  =  v$  t. 

Combining  the  above,  one  finds 

$s  s  (ilp/vs)  (l+ap/2). 

Letting  <o s ^  a  (Q/v  )  (l+a/2),  the  above  equation  may  be  expressed  in 
terms  of  <<>s^  to  obtain 

<t>s  =  ^sl  D(^+at>/2)/(l+«/2) , 

4>s  =  't>s^(l+ap)/(l+a/2)  and 
<»s  3  <i» s ^  <*/ ( l+a/2 ) . 

Substituting  respective  terms  in  Eq.  (16)  by  the  above  results,  the  coordinate 
transformation  of  bending  spray  plumes,  as  specified  here,  can  be  achieved. 

Regarding  the  plumes  with  their  axis  having  an  inclined  angle  of  e.  ,  the 
distribution-to-distribution  transformation  can  be  obtained  by  introducing 
Eqs .  (11)  and  (12)  into  Eq.  (15). 

COMPUTATION  OF  RAOIATION  IN  OIESEL  C0M8USTI0N 

The  present  method  may  readily  be  used  for  analysis  of  radiation  heat 
transfer  processes  in  combustors  having  jet  flames.  An  analysis  may  be  made 
over  the  entire  wall  of  the  reactor  to  introduce  the  results  to  implement  a 
global  analysis  of  the  system  heat  transfer  Including  convective  and  conductio 
heat  transfers,  which  will  enable  prediction  of  the  temperature  distribution 
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into  and  over  the  wall.  When  an  unsteady  combustion  reactor  is  considered,  e.g., 
a  diesel  engine,  the  analysis  will  have  to  be  made  over  the  entire  cycle  period 
along  with  energy  conversion  process,  as  presently  attempted  in  our  laboratory. 
Tne  present  section,  however,  considers  the  computation  of  radiation  heat 
transfer  at  a  particular  engine  crank  angle.  Since  more  extensive  results  are 
reported  elsewhere  [22],  only  a  few  computations  are  described  here.  Among  the 
assumptions  employed  due  to  the  absence  of  either  sufficient  engine  data  or 
suitable  engine  modeling  are  (1)  the  values  of  constants  in  Eq.  (8)  are  for 
soot,  H^O ,  CO^  and  burned  fuel/air  ratio,  8xl0~°,  0.01,  0.01,  and  1.0, 
respectively;  the  distribution  constants,  a,  b  and  c  are  0.6,  2.36,  and  2  for 
the  above  cases  accordingly.  (For  example,  the  CO^  distribution  in  the  plume 
may  be  written  as  F  =  {0.01  Exp  -  [0.6(o/R)  +  2. A3 (®/$  )  +  2 - 0( z/R ) ^ ] } ,  where 
R-is  the  piston  radius  and  <t>r  is  */(number  of  plumes),  (2)  other  engine 
details  are,  the  surface  temperature,  500°K;  the  number  of  plumes,  4,  the  piston 
radius,  4.92  cm;  the  combustion  chamber-  bowl  radius,  3.44  cm,  and  the  fuel 
compos i tion  c16H34’  (3  )  the  computation  of  the  temperature  distribution  is 
computed  by  using  the  equilibrium  adiabatic  temperature  relationship  to  the 
fuel/air  ratio  (e.g.,  using  a  new  method  [23]),  etc.  The  computed  radiation 
heat  flux  is  presented  presented  in  terms  of  the  following  nondimens i ona 1 i zed 
terms:  the  normalized  spectral  energy,  Q^/Q^.  and  the  normalized  directional 
radiation  intensity,  I  (  o  ,  0  ,  z ) /I  b ,  where,  Qfar  =  =  2,400°K,  I  =  Q  /*, 

and  a  is  the  Stefan-3ol tzman  constant  ( 5 . 670xl0~^W/m3  ) . 

Shown  here  are  the  following  case  analyses:  (1)  the  spectral  radiation  heat 
incident  upon  locations  on  the  cylinder  head  along  the  axis  of  a  spray  plume  at 
r  =  0.0R,  0.5R  and  0.9R  (Fig.  3),  (2)  the  normalized  directional  radiation 
intensity  streaming  into  the  Injection  nozzle  tip  {o  =  0.0,  <i  =  0.0)  along  tne 
direction  or  zenith  angle  in  the  hemispherical  volume  faced  Dy  the  tip  '~;o.  4. 
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and  (3)  the  normalized  directional  radiation  intensity  along  the  locations  on 

the  cylinder  head  right  above  the  axis  of  the  spray  plume  in  directions  of 

zenith  angle,  9  =  17.6,  60.0  and  76.7  in  the  hemispherical  volume  faced  by  each 

location  (Fig.  5).  Without  discussing  the  details  of  the  results,  one  may  readily  flnd 

from  the  above  illustration  that  improved  radiation  heat  transfer  computations 

can  be  achieved  by  employing  the  present  new  coordinate  transformati on  method 

in  its  model Ing. 
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NOMENCLATURE 

A  band  absorptance 

a  constant  determining  species  distribution 
8  radiosity 

b  constant  determining  species  distribution 
c  speed  of  light  or  constant  determining  species  distribution 
e  Plank’s  emission  power 
F  species  distribution  function 

f  soot  volume  fraction  or  transformation  function 
h  Plank's  constant 
I  radiation  intensity 

k  8oltzman  constant  or  imaginary  part  of  reflective  index 
n  real  part  of  the  reflective  index 
Q  black  body  radiation 


q  radiation  neat  flux 


R  radius  of  cylinder 
T  temperature 


t  time 
v  velocity 

x  component  of  cartesian  coordinate  X 
y  component  of  cartesian  coordinate  V 

a  constant  for  plume  injection  speed 
s  reflectivity  of  cylinder  wall 
r,0,t  components  in  spherical  coordinate,  C 
c  emittance 

e  angle  of  optical  path  with  respect  to  normal 
<  volume  absorptance 
\  wavelength 

o,<t,Z  components  in  cylindrical  coordinate,  S 
x  extinction  coefficient 
<&  transformation  function 
Q  angular  speed 

Subscripts 

b  blackbody 
d  detector 
g  gas 

in  inclined 
o  constant 
s  soot 


t  total 


20 


u  unparallel 
\  spectral 

9s  optical  path  of  solid  angle 
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ABSTRACT 

A  new  Integral  function  has  been  discovered.  The  present  function 
la  found  to  be  useful  in  computing  the  spectral  emlssivlty  of  an 
isothermal  volume  containing  either  soot  or  gaseous  species,  or  both. 
Examples  of  its  application  are  discussed  herein. 


Introduction 


In  the  application  of  the  Planck  formula  for  computing  the  energy  radi¬ 
ated  from  an  isothermal  source,  it  is  necessary  to  find  the  emlssivlty  of  the 
source.  The  emlssivlty,  c,  is  expressed  in  terms  of  its  spectral  amissivity, 
e^,  as  follows 


e 


(T)dA/<JT4 


(1) 


where  E^  (T)  is  the  Planck  radiation  function 
,  .  2 

2  she _ 


E  (T)  -  -T- 

X3  [exp  (hc/lKT)  -1] 

c,  X  and  h,  constants. 


In  the  past,  the  determination  of  has  been  made  for  various  radiating 
sources  by  either  experimental  or  theoretical  means,  e.g.,  for  solid  surface 
[1]  ,  gases  (2,3],  soot  suspension  in  transparent  gases  [4, SI,  and  combustion 
products  (6,8].  For  the  computation  of  of  an  isothermal  volume  with  an 
optical  length,  L,  containing  radiating  gases  and/or  soot,  the  following 
equation  was  alloyed  (1] 
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where  k is  the  spectral  volume  absorption  coefficient. 

Regarding  the  emiasivity  by  soot,  tj(  one  may  introduce  eg.  (2)  into  eq. (1) 
to  find ' 


e 

s 


1 


E^mdl/oT4 


(3) 


Although  significant  progress  has  b«en  made  on  the  determination  of  [3-8]  , 
researchers  have  resorted  to  numerical  methods  in  finding  a  solution  of  eq. (3) 
with  additional  assumptions.  A  coemon  assumption  made  in  the  past  [6,91  is 
the  relationship,  ■  Kjf^/X^,  where,  <9,  the  constant)  f^,  the  soot  volume 

fraction;  and  d,  the  number  chosen  close  to  1.  This  allowed  one  to  obtain  the 
final  equation  of  emiasivity  including  the  hexagantna  function 


H(L)  -  -  t4e-Lt/(l-e-t)dt  (4) 

■  « 

where  L  represents  the  space  variable.  Consequently,  the  final  equation  was 
numerically  solved,  e.g.,  using  Simpson's  rule  [6,9]. 

The  New  Integral  Function  and  Applications 

The  present  paper  offers  an  exact  solution  of  eq.(3).  The  spectral  ab¬ 
sorption  coefficient  of  an  isothermal  volume  containing  soot  that  was  found, 
by  using  the  dispersion  theory  [4,5,7]  and  the  Ray leigh- limit ,  was  introduced 
into  eq. (3) .  It  was  followed  by  a  rearrangement  to  obtain  the  following 
tidy  equation 


c 

s 


where 


x 

a 


-ax 


a 


3 

x 

X 

e  - 


1 


dx 


hc/XKT  , 

36im2V 

[n2-nk2*2]-Mn4)t2 


f  tS 

v  ch 


and 


(5) 


n+lnk,  the  complex  refractive  index. 

An  exact  solution  of  eq. (5]  has  been  discovered  by  the  authors  in  a  closed 
form  by  defining  a  more  general  expression  of  indefinite  integration,  F(m,a,x) 


’  -ax  m 

r(m,a,x)  »  -  -  dx  (a  >  0,  a,  positive  integer) 

ex-l 


m  -<a+l)x 

x  e 


(1-e  ) 


dx 
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x»#“(»+l)x  r  -nx 


r  -nx  . 

i  «  Ax 

n«0 


*  I 

n»l 


-(n+a)x  m  . 

•  x  dx  . 


By  using  the  Heaviside  linaar  diffarantial  inversion  operator,  D ,  the  authors 
found 


F(m,a,x) 


r  -(n-ta)x  n, 
x  ] 

-  (n+a)  x  1 

B-  (n+a) 


31, 
X  1 


l 

J 


l 

n*  1 


-<n+a)x  *  i 

a  *  *“> 


(n+a) 

Finally,  the  indefinite  solution  was  obtained  as 


F(m,a,x)  » 


l  {' 

n-1  *■ 


-<n+a)x 


rtm- 


rtm+l)  ,-i-l  m-i, 

-  -  (n+a)  x  1 


i-0 


i+1) 


where  i"  is  the  gamma  function  defined  as 
(■(m+1)  ■  ar  (m)  and  T(l)  »  1. 

The  definite  integral  in  eg. (S) ,  therefore,  becomes 

x- 


/  i  -ax  a 
I  a  x 

x  , 

x^  e  -1 


dx  ■  F(m,a,x^)  -  F(m,a,x.,). 


(6) 


(7) 


In  fact,  the  solution  for  the  hexagamma  function,  eg.  (4)  ,  can  be  readily 
found  as 


H(L)  -  -  [F (4 ,L,0)  -  F(4,L,»)1 . 


The  new  solution  can  be  conveniently  used  even  in  the  computation  of  the 
gas  emissivity,  that  employs  the  use  of  the  band  model  [2,6]  to  have 


c 


g.i 


IS 

4 


J 


band  i 


(8' 


Although  the  solution  of  eg. (8)  may  be  approximated  by  using  tabulated  values 
of  Debye  functions  (101  or  other  approximations  [11],  the  new  solution, 
eg.  (61  ,  offers  a  snore  convenient  and  accurate  result,  l.e., 
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Hare,  u  and  1  Indicate  the  upper  and  lower  bound  of  the  band. 

In  addition,  it  was  found  that  the  new  solution  enables  the  authors  to 
construct  an  extensive  three-dimensional  model  of  radiation  in  diesel  combus¬ 
tion  as  long  as  in  eq.  (3)  is  known  in  terms  of  L. 
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Nomenclature 

a  function  defined  by  eq. (5) 

c  speed  of  light 

D  Heaviside  linear  differential  inversion  operator 
d  number  chosen  close  to  1 

E^  Planck's  spectral  black  body  radiation  function 

F  indifinite  integral  function  defined  by  eq. (6) 

f  soot  volume  fraction 

v 

g  subscript  denoting  gas 

H  hexagamma  function  defined  by  eq.  (4) 

h  Planck 1 s  constant 

K  Boltzmann ■ s  constant 

L  space  variable 

m  positive  integer 

n+ink  complex  refractive  index 
s  subscript  denoting  soot 

T  temperature 

t  variable  used  in  eq. (4) 

x  variable  defined  by  eq. (S) 

c  emissivity 

r  ganma  function 

c  Stef an-Boltzman  constant 

spectral  volume  absorption  coefficient 
< ,  constant 
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ABSTRACT 

The  present  paper  reports  an  exact  solution  of  the  fractional 
function  of  blackbody  radiation.  The  new  equation  Is  expressed 

as 


F 


o-AT 


/IT 


VT> 


d(AT) 


oT 


where  x  =  hC  /XkT. 

o 


Introduction 

It  Is  often  necessary  to  calculate  the  fraction  of  the  total  radiation 
emission  from  a  blackbody  chat  Is  emitted  In  a  given  wave  length  lncerval  or 
band.  This  fraction,  for  a  given  temperature,  T  and  the  wavelength  Interval 
from  o  to  \,  may  be  expressed  as 
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where,  E^  Is  the  Plank  radiation  function 

„  2i»h  C2 

v  »  _ — - - 

0  X  t  exp (hC  /XkT)-l] 
o 

C,  C  ,  k  and  h  are  constants, 
o 

Due  to  the  unavailability  of  an  exact  solution  of  eq.  (1)  in  the  past 

and  Its  frequent  uses  In  radiation  heat  transfer  computations,  F  has 

o- AT 

been  numerically  obtained  and  Its  tabulated  values  for  different  IT  are  found 
in  most  of  the  fundamental  heat  transfer  textbooks.  In  ac dition,  for  the  com¬ 
puter  solution  of  various  types  of  radiation  problems,  Plvovonsky  and  Nagel  [2] 
and  Wlebelt  [3]  have  presented  polynomial  curves  fitted  to  the  function,  e.g.. 


15  -®v 

F0_xt  "  — —  l  — { [  (mv+3)mv>+6  lmvN-6) 

»  m«l  ,2,  . .  m 


for  v  >  1, 


o-XT 


2  4  6  8 

1  v  _1_  *  v  a  V  V  x  t  ^ 

ff4  ‘  8  60  ■  5040  272160  "  13305600)  for  v  <  2 


where  v  i  —  and  -  14388  um.X  . 


The  New  Blackbodv  Radiation  Functions 
The  present  paper  reports  an  exact  solution  of  eq.  (1).  The  solution  has 
been  derived  by  a  direct  application  of  the  new  Integral  function  thac  the 
authors  have  recently  presented  demonstrating  its  usefulness  in  thermal  radia¬ 
tion  calculations.  The  reported  integral  is  expressed  as  [4] 
rX 2  -  ax  m 

J  ~  'x~*-  dx  -  F(ra,s,x1)  -  F(m,a,x2)  (2) 


where,  F(o,a,x) 


l  (« 

i-l 


■(n+a)x 


l  l 


<.-1 


r (©♦•I)  ,  .  ,-c-l  m~t, , 

rT£=3T)(n+a)  x  » 


r(art-l)  »  or(m)  and  r ( 1 )  -  1. 
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Sine*  eh«  speed  of  light,  C0  ■  vX,  where  v  Is  the  wave  number,  eq.  (1) 
con  be  rewritten  as 


fAT  hx  .  fv  Sv 

— r  d(XT)  -  1  -  — J 

aTJ  J A  flT 


2*hv 

where-  ^V  “  TT. - TTT-- 

Co  (exp  (^)-ll 


And  since  E,  dv/<7T4  »  1,  by  using  eq.  (2)  with  m  -  3  and  a  «  0,  we  find 
o  ov 


,  ,  4  (X  3 
2  irk  x 


_  2,  3  J  x 
oC  h  ■'o  e  -1 
0 


Hence, 


„  2.  3 
oC  h 
o 


2irk4 

oC  2h3 
o 


2  ,V 
15  c  V 

o 


3!  •  I  4 


3!  •  W  *  1  * 


The  Integral  In  the  righthand  side  of  eq.  (3),  therefore,  can  be  expressed  by 
using  eqs.  (2)  and  (4),  with  x  i  hCWXkT,  as 

fV  j  15  f*  x3  . 

— r  dv  •  —  -  dx 

'  o  oT  it  o  eX-l 


,  e  -  -n* 

3  -  H  nv* 

it  n-1 


n  n 


Finally,  from  the  relationships  In  eqs.  (3)  and  (5),  we  find  the  exact  form 
of  the  fractional  function  of  the  blackbody  radlacion  as 
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Fig.  1  Blackbody  Radiation  Functions 
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Nomenclature 

function  used  in  eq.  (2) 
speed  of  light 

Plank's  spectral  blackbody  radiation  function 
integral  function  defined  by  eq.  (2) 
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Fq  fractional  function  of  the  blackbody  radiation 

h  Flank's  constant 

k  Boltzman's  constant 

m  positive  Integer 

n  variable  used  In  eq.  (5) 

T  temperature 

x  variable  defined  In  eq.  (6) 

r  gamma  function 

o  Stef an-Boltzman  constant 

X  wavelength 

v  wavenumber 
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Abstract — This  paper  reports  new  formulas  of  aitiabauc  flame  temperature.  They  arc  functionally 
expressed  in  terms  of  fuel'air  ratio,  reaction  pressure,  and  the  number  of  carbon  atoms  m  me 
individual  fuel.  Among  the  fuels  presently  considered  for  the  formulas  are  members  of  parailin. 
aromatic  and  olefin  families,  acetylene,  alcohols,  and  hydrogen. 


INTRODUCTION 

The  combustion  reactions  in  many  air-breathing  engines  take  place  at  high  tem¬ 
peratures.  When  the  energy  released  in  such  reactions  is  completely  used  to  heat 
the  combustion  products,  the  product  temperature  may  attain  the  highest  level.  At 
such  a  temperature,  however,  the  combustion  products  dissociate  into  fractions  of 
molecule,  c.g.,  atom  and  free  radicals,  by  up  to  several  percentage  points.  Since  the 
dissociation  reaction  is  highly  endothermic,  the  processes  substantially  lower  the 
flame  temperature,  consequently  achieving  the  equilibrium  adiabatic  flame  tem¬ 
perature,  Ta. 

NASA-Lewis  pioneered  the  computation  of  T„  and  thermodynamic  equilibrium 
composition  of  combustion  products  of  various  fuels  (Gordon  and  McBride.  197]). 
Among  the  same  computational  approaches  taken  by  others  for  the  equilibrium  data 
of  gaseous  hydrocarbon-air  systems  are  those  for  modeling  the  combustion  processes 
in  internal  combustion  engines  (Olikara  and  Borman.  I975)  and  for  computing 
radiation  heat  transfer  in  diesel  combustion  (Chung  and  Rliee.  1983).  In  computing 
the  equilibrium  data  of  a  specific  fuel-air  system,  eleven  gaseous  species  were  con¬ 
sidered  beside  the  fuel,  including  seven  dissociation  reactions  with  known  individual 
equilibrium  constants,  a  stoichiometric  equation  and  an  energy  equation.  Due  to 
the  highly  nonlinear  nature  of  the  problem,  the  solution  was  numerically  found  by 
using  Newton's  iteration  method  (Gordon  and  McBride.  1 97 1  ;  Olikara  and  Borman. 
1975).  Although  the  authors  employed  a  new  numerical  method  for  the  solution 
that  substantially  reduced  the  computation  time,  one  of  the  more  time-consuming 
steps  for  a  recent  report  (Chang  and  Rhce.  1983)  was  still  that  for  finding  T„  at  each 
nodal  point  with  a  burned  fuel/air  ratio  in  the  combustion  volume.  In  order  to 
simplify  the  computational  process  of  T„  in  their  more  involved  on-going  compu¬ 
tational  model  of  diesel  combustion,  the  authors  have  developed  simple  formulas  of 
Ti  for  various  fucl-air  systems. 
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Tim  paper  reports  the  new  formulas  of  T„  for  an  individual  family  of  fuels  expressed 
in  terms  of  the  fuel  air  ratio  and  the  reaction  pressure.  The  formulas  were  determined 
by  using  the  least-squares  method  for  computed  results  obtained  from  the  above 
mentioned  method  with  variation  of  fuel  (gas  slate),  fuel  air  ratio,  and  pressure 
Note  that  the  present  formulas  may  not  be  useful  for  those  beyond  the  lean  llam 
mahility  because  the  flame  does  not  virtually  exist  in  the  mixtures,  and  that  the 
formulas  for  those  richer  than  stoichiometric  are  not  considered  due  to  the  unavail 
ability  of  the  equilibrium  dissociation  equations  for  such  mixtures.  The  new  formula 
of  r„  of  a  fuel-air  system  leaner  than  stoichiometric  under  reaction  pressure.  P.  is 
expressed  as. 

T„  =  «(l  -  />  •  In </>  -  c  ■  (In./.)-]  (  K) 

where  a  =  </i  -  o-_»  •  ln(  P  P,>) 

h  =  />,  -  /».,  In (/’  /’„) 

t  =  <  i  •  ln(  P'Pn) 

Fuel  volume  fraction 

4  - - 

(Fuel  volume  fraction )m.„.-u 

Pn  =  I  atm 

iii.  h\. .  .  .  arc  coefficients  of  each  family  of  fuels  given  in  the  following. 

1)  Paraffin  family.  C„Hj„  :  (I  ^  n  $  16) 

in  =  2295.0  +  9.5  ln</»  -  0.99) 

<h  =  13.2  +  0.7  In («  -  0.97) 

/>,  =  0.544  -  0.0019  ln</»  -  0.99) 
h*  =  0.018  +  0.0009  ln</»  -  0.96) 
r,  =  0.070  -  0.003  ln<«  -  0.98) 
r.  =  0.0169  -  0.0008  ln(/»  -  0.97) 

2)  Olefin  family.  (2  ^  n  ^  16) 

ay  =  2377.0  -  17  9  Info  -  1.85) 
a>  =  19  9  -  1.4  InOr  -  1.89) 
hi  =  0.512  +  0.0085  ln(«  -  1.89) 
hi  =  0.025  -  0.0014  I n(/f  -  1.86) 
d  =  0.035  +  0.0083  I n(/i  -  1.89) 

<•.  =  0.023  -  0.0013  ln(ji  -  1.86) 


3)  Aromatic  family.  CiH-:,, 


(6  $  $  22) 


cn  =  2378.0  -  15.3  !n<//  -  5.36) 
o3  -  21.0  -  1.6  Into  -  5.29) 
hi  =  0.518  *  0.006  I n(/r  -  5.51) 
h2  =  0.027  -  0.002  ln(n  -  5.24) 
n  =  0.034  -  0.008  ln(/t  -  5.38) 
c2  =  0.025  -  0.0016  Info  -  5.27) 
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4)  Alcohol  family,  C„H>  tOH  (2  $  «  $  IS) 

«,  =  2153.5  -  60.2  I n(//  -  1.51) 

«...  =  -0.5  -  5.2  ln(/i  +  0.11) 

/>i  =  0.641  -  0.030  I n(//  -  3.4) 
h-i  =  0.0004  -  0.0066  ln(/t  0.02) 
c,  =  0.193  -  0.041  ln(ir  -  3.5) 
t  .  =  -0.0002  -  0  0066  ln(/r  -  0.09) 


5)  Hydrogen.  H- 

«,  =  2412.2 
a-  =  17.9 
/>i  =  0.473 
/>,.  =  0.021 
n  =  0.0217 
=  0.0191 


6)  Acetylene.  C:H: 

,/,  =  2562.9 
«...  =  42.3 
/»,  =  0.390 
h.  =  0.0426 
r,  =  -0.0X19 
,  .  =  0.0379 


A  comparison  of  computed  T«  and  approximated  T„  hv  the  present  formula  oi 
para fTin-air  systems  at  reaction  pressure  of  100  atm  can  be  seen  in  Figure  I.  In 
almost  all  the  cases  with  the  formula,  the  discrepancy  between  the  two  sources  of 
value  was  less  than  one  percent.  The  formula  enables  one  to  find  the  adiabatic  flame 
temperature  of  a  wide  range  of  fuels,  fuel  air  ratio,  and  reaction  pressures  even  b\ 
using  a  simple  pocket  calculator. 


Computed,  n«l 
Approximated,  n«l 
Computed,  n«8 
Approximated,  n*8 
Computed,  n«  16 
Approximated,  n*l6 


Fuel  Volume  Fraction  of  Stoichiometric 

FIGURE  I  Adiabatic  llame  temperature  of  paralTm  fjiiul>.  C.Hj.,  :.  as  function  ol  tucl  an 
ratio  and  number  of  carbon  aioms  of  fuel  for  reaction  pressure  ot  100  atm. 
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Erratum 


1  ffrt  curt/  4ftn/  2  * 

Chany,  b  i.  ..mi  Rhcc.  K  I  (I9K3)  Adiabatic  flame  temperature  estimates  of 
lean  tuel  air  mixtures  Sliort  CommunietiMon.  C»mhtisimn  Science  anil  Tech- 
Hum” >'  35.  ’(Is 


In  onr  resent  papei  e  reported  ness  formulas  lor  adiabatic  flame  temperature  (AFT) 
’.■'pressed  m  terms  ol  the  fuel  air  rain)  (entering  at  29fc  K  ).  the  reaction  pressure,  and 
•'v  main  he:  <>  .  arbon  atoms  in  the  individual  fuel  Recently,  Gulaer  (It  has  per- 
'ormeJ  .'Insc  -tud>  f>t  the  formulas,  which  lead  us  to  diw.'.'Ci  a  numerical  error  in 
.ops  me  :;u  ncislcnc  heal  of  formulation  in  the  computei  .ode  In  addition,  there 
a  as  .it:  es!i  i  or 'a  I  error  in  our  computer  program  of  the  A  f  T  tor  alcohol  -air  mixtures 
The  nrre'.t  Vrmul.ts  lor  hoth  cases  are  eivr.ii  helow.  In  the  meantime,  the  authors 
nave  .omplcted  additional  work  i'i  id  find  new  improved  formulas  of  the  AFT 
similar  to  the  previous  ones  \vi<h  inclusion  of  the  initial  temperature  as  a  variable 

\celvlene.  (  IF  Alcohol  family  C„FI>„OH  (2  s-  n  X) 

-  -523  1  at  •  2276.X  +  12.6  ln(/t  -  1.5 1) 

•  i-  44  5  a-:  -  1 2.65 -i-  1 . 147  ln(/i  +  0. 1 1  > 

•  0  374  5,  =0.4677  +  0.0261  ln(/i  t- 3.4) 

«439  />•_•  -  0.01 7.3  +  0  0013  ln(x  •+■  0.02) 

0.094*  .  ,  ;  0.0272 -Ml  01 22  ln(n  t  3.5) 

-0  03X6  ,  .  -  ft  0163  e  0.001  7  ln<»  i  0  09) 
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Abstract — A  new  simple  formula  for  calculanng  the  adiabatic  flame  temperature  of  fuel  air 
mixtures  has  been  developed.  The  formula  is  usable  In  a  wide  range  of  both  fuel  and  reaction 
variables.  Among  the  fuels  considered  for  the  formula  are  members  of  the  paraffin,  aromatic, 
olefin  and  alcohol  families,  and  hydrogen  and  acetylene.  The  formula  is  functionally  expressed 
in  terms  of  the  fuel  air  ratio,  the  reaction  pressure,  the  initial  nixiure  temperature  and  the  number 
of  carbon  atoms  in  the  individual  fuel.  By  using  the  formula,  the  adiabatic  llame  temperature  for 
a  designated  fuel  and  reaction  condition  can  be  found  within  an  accuracy  of  one  percent. 


INTRODUCTION 

The  computation  of  the  equilibrium  adiabatic  flame  temperature  is  often  performed 
in  the  analysis  of  various  airbreathmg  combustion  devices.  Due  to  its  complexity, 
however,  researchers  have  resorted  to  computer  operation  for  ns  computation  in 
the  past  (Gordon  and  McBride,  1971 ;  Olikara  and  Borman.  1975;  Chang  and  Rhec. 
1983).  The  computation  for  a  typical  fuel-air  system  considers  eleven  gaseous 
species  besides  fuel,  in  simultaneously  solving  the  following  equations-,  seven  dis¬ 
sociation  reaction  equations  with  the  respective  temperature-dependent  equilibrium 
constants;  a  stoichiometric  equation  and  an  energy  equation. 

In  spite  of  a  considerable  reduction  in  computer  time  by  the  use  of  a  new  numerical 
method  for  the  solution,  one  of  the  more  time-consuming  steps  for  a  recent  paper 
(Chang  and  Rhce,  1983)  was  still  that  for  finding  the  adiabatic  llame  temperature  in 
each  nodal  point  with  a  burned  fuel-air  ratio  in  the  combustion  volume.  In  order 
to  facilitate  its  involved  computation,  the  authors  recently  introduced  simple  formu¬ 
las  for  adiabatic  llame  temperature  estimates  (Chang  and  Rhee.  1983).  They  arc 
expressed  in  terms  of  fuel  air  ratio  (entering  at  298  K).  reaction  pressure,  and  the 
number  of  carbon  atoms  in  the  individual  fuels.  For  a  similar  goal,  graphs  for  adia¬ 
batic  flame  temperature  of  hydrocarbon  fuel -air  mixtures  were  introduced  by  others 
(Glassman  and  Clark.  1983).  As  a  continuation  of  the  previous  endeavor,  the  authors 
have  developed  a  new  improved  formula  for  adiabatic  llame  temperature  which 
additionally  includes  the  initial  reactant  temperature  as  a  variable  and  generate 
more  accurate  results.  In  addition,  as  an  exploratory  procedure  in  the  development 
of  more  accurate  and  consistent  formulas  of  adiabatic  llame  temperature,  factors 
controlling  the  temperature  in  the  First  Law  relation  were  reviewed.  Sonic  findings 
seem  to  be  worth  reporting,  along  with  the  newly  developed  formula. 


Adiabatic  Flame  Temperature  of  Fuel- Air  Combustion 


|t  is  generally  expected  that  the  energy  released  in  a  combustion  reaction  is  best 
represented  by  the  adiabatic  flame  temperature  tGaydon  and  Wolfhard.  1970)  In 
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the  calculation  of  the  heat  of  reaction,  as  a  convenient  and  logical  procedure,  the 
heat  of  formation  at  the  standard  state.  A H'  -w.  of  chemical  substances  arc  used. 
To  a  good  approximation,  in  many  cases,  it  is  possible  to  determine  A//  joh  of  a 
molecule  by  adding  the  energy  of  bonds  forming  the  molecule.  The  heat  of  reaction 
(or  heat  of  combustion),  Q„ ,  of  a  fuel  may  be  determined  by  using  A/Z'jsg  of  sub¬ 
stances  involved  in  the  reaction  as  follows 

-  QP  =  I  .»■<( A//'  »os)i  -  £  iv(A//  .98)/.  ( I ) 

/-products  ./-reactants 

For  a  closer  look  at  the  above  relationship  of  common  fuels,  the  molar  quantities 
of  A H'ioa  and  Qp  at  a  reference  temperature  of  298 :  K  arc  plotted  againsr  the  number 
of  carbon  atoms,  ,Vf,  in  individual  fuels  as  shown  in  Figures  I  and  2.  respectively 
A  brief  look  at  the  Figures  and  Eq.  (I)  shows,  although 'not  new,  that  A/f  ;<»  of 
fuel  is  much  smaller  than  the  corresponding  Q,,  in  magnitude  and  that  for  the  respect¬ 
ive  fuel  family  there  is  no  mutually  comparable  trend  between  AF/'mim-V,  and 
QP-Sr  relationships.  This  indicates  that  A/f  for  fuels  may  not  play  a  great  role 
in  determining  the  adiabatic  tlame  temperature.  T„  For  the  purpose  of  comparison 


FIGURE  I  Molar  heat  of  formation  of  IlicK  i,.  number  of  .jihnn  junn- 
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77 


Number  of  Carbon  Atoms  in  Fuel 

FIGURE  2  Molar  hear  of  combustion  of  fucl-air  vs.  number  of  carbon  aioms  in  fuel. 

jnd  then  finding  the  possibility  of  using  Q,,  ns  a  parameter  in  the  new  formulas  of 
T,  the  relationship  of  T«  to  yr  is  reviewed  for  fuel-air  mixtures  (initially  at  298  K. 
and  reacting  at  I  atml  at  a  fuel  volume  fraction  of  stoichiometric.  <h=  I  (in  addition, 
for  alcohol.  n»0.7Sl.  as  shown  in  Figure  3.  Putting  aside  the  mutual  irrcspectivity 
between  the  2 \H  .^g-  V,  and  7",,-  Vr  relationships,  the  similar  irrcspcctivity  between 
the  'Jj-''  -  and  77,- V-  relationship  could  be  surprising.  But  one  soon  finds  that  Q,, 
nav  not  he  a  proper  indicator  for  77,  by  recognicmg  that  Q „  is  used  up.  in  a  highly 
nonlinear  manner,  m  heating  the  products  and  the  endothermic  dissociation  reaction 
as 

I  -\\H  u- H~ :,»)-■  r,-H  ..ps) -  A//  =  -  O,,.  (2) 

,-proOucts  ./-reactants 

wnere  T  ndicatcs  the  initial  temperature  of  reactants,  r  is  the  number  of  moles  of 
each  species,  and  Q ,,  becomes  zero  in  an  equilibrium.  Further,  since  the  analysis  of 
the  open  thermodynamic  systems  is  more  conveniently  performed  by  using  mass 
thermodynamic  quantities  than  molar  basis,  the  mass  quantities  of  ±11  and  O ,, 
may  be  worth  investigating  as  shown  in  Figures  4  and  5.  Unlike  the  hitherto-held 
notion  that  ±H  may  not  well  represent  77,.  it  is  rather  unexpected  to  lind  a  mutual 
similarity  m  the  trend  existing  between  the  ±11  „••!*( mass)-  V,  and  7’,,-  V..  relation¬ 
ships  except  for  the  alcohol  fuel  family;  no  such  trend  exists  between  the  O,,  (massi  V 
and  77,- V,  relationships.  Although  there  is  no  explanation  that  can  be  given  for  the 
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unexpected  mutual  similarity,  this  trend  seems  to  support  choosing  A//  as  a 
parameter  in  the  construction  of  7',,-graphs  introduced  by  Classman  and  Clark  ( 198.'). 
In  search  of  a  proper  formula  of  T„,  however,  the  close  relationship  of  A//  to 
7"„  found  in  the  present  study  is  not  used  due  to  the  inability  to  include  the  alcohol 
family. 

In  parallel  with  the  above-mentioned  approach,  an  extensive  parametric  study  of 
T«  has  been  carried  out  for  the  following  major  variables  by  the  numerical  method 
for  T«  mentioned  earlier:  the  fuel/air  ratio;  the  reaction  pressure  and  the  initial 
temperature  of  reactants;  and  the  number  of  carbon  atoms  in  individual  fuels. 
Some  of  the  results  are  summarized  in  the  following  forms:  adiabatic  flame  tem¬ 
perature  vs.  fuel  volume  fraction  of  stoichiometric  at  varied  reaction  pressure 
(Figure  6);  adiabatic  flame  temperature  vs.  initial  mixture  temperature  for  rep¬ 
resentative  fuels  at  reaction  pressure  of  I  atm  (Figure  7);  and  adiabatic  flame  tem¬ 
perature  vs.  reaction  pressure  with  varied  initial  mixture  temperature  (Figure  8). 


Initial  Mixture  Temperature,  K 


FIGURE  7  Adiabauc  flame  temperature  of  various  Fuels  iuuh  s  I  0  rc.iumu  .u  1  .nun 
initial  mixture  temperature. 


Among  the  important  characteristics  of  the  results,  already  well  known  arc  the 
following. 

I)  The  higher  the  reaction  pressure,  the  higher  the  adiabatic  flame  'cmpci.iunc 
Tilts  is  due  to  the  more  reverse  of  dissociation  reactions  iChane  and  Race  . 


Adiabatic  Flame  Temperature, 
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Number  of  Carbon  Atoms  in  Fuel 


5  Mass  heal  of  combustion  of  fuels  vs.  number  of  carbon  atoms  in  fuel 


Ertcct  of  fuel  air  ratio  and  reaction  procure  on  ailtahaiK 
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FIGURE  8  EfTeci  of  reaction  pressure  and  initial  mixture  temperature  on  adiabatic  flame 
temperature. 


higher  reaction  pressure,  in  particular. 

CO>  -CO-iOi 
H20  —  H2-402 
*H2  —  H 
l 0  >  -»  0,  and 
Jn2  —  N. 

21  The  higher  the  adiabatic  flame  temperature  of  the  mixture,  the  greater  the  effect 
of  reaction  pressure  on  the  temperature.  The  effect  is  greater  because  of  the  greater 
tendency  of  dissociation  at  higher  flame  temperature. 

3)  An  increment  in  the  initial  mixture  temperature  by  AT  is  reflected  in  the  adia¬ 
batic  flame  temperature  by  an  amount  much  smaller  titan  A 7".  this  tendency  is  greater 
in  higher  initial  mixture  temperature  and  for  mixtures  with  higher  llanie  temperature 
This  is  explained  by  the  greater  tendency  ol  dissociation  jt  higher  flame  temperatures, 
and  the  change  in  heat  capacity  with  temperature  and  composition 

A  careful  review  of  results  (refer  to  Figures  3.  '  and  X)  revealed  that  the  adiabatic 
(lame  temperature  could  be  expressed  in  terms  of  the  above-mentioned  maior  variables 
by  using  the  least-squares  fit  method,  as  similarly  attempted  by  the  present  authors 
(Chang  and  Rhee,  19X3).  The  mam  task  for  this  was  to  idennfy  us  equauonal  lorm 

Veu  Formulas  for  -I diabolic  Flume  Temperature 

The  search  for  new  formulas  for  T„  of  fuel-air  systems  was  made  in  such  a  wav  .i- 
to  meet  the  following  requirements,  ilia  simple  form  of  equation  >2)  convenient 
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and  easy  to  use;  (3)  accuracy  within  one  percent;  (4)  applicable  within  the  flammability 
limit  of  individual  fuels  but  leaner  than  stoichiometric;  (5)  applicable  in  the  range  of 
initial  mixture  temperature  from  298CK  to  the  self-ignition  temperature  of  individual 
fuels,  (6)  applicable,  at  least,  in  the  range  of  reaction  pressure  from  I  atm  to  100  atm, 
and  (7)  mixtures  in  gas  state. 

Upon  performing  a  considerable  amount  of  computational  work  in  correlating  7"a 
to  the  previously  described  main  variables,  largely  through  the  trial  and  error 
approach,  the  authors  have  found  the  following  equational  form  which  oilers  a 
considerable  flexibility  in  use  and  to  meet  most  of  the  abovementioned  requirements. 
The  new  formula  of  Ta  is 


Ta  =  ,4i[l  4-  He  If*  i  ■+■  4j(ln  'A)~], 

(31 

A, 

=  Hn  I-.- 4,2  In  P  +  At s(\n  P)2  (/  = 

1,2.  3) 

(41 

An 

=  An i  +  An  ’  In  Ft  .4,,3<ln  T)-  (j  = 

1.  2.  3) 

(51 

Ain 

=  Atjk i  —  In  ,Ve  —  4 A3(  1  n  ,VC)J- 

II 

1  J 

(6) 

Ta.  the  estimated  adiabatic  flame  temperature  |'K.) 

■A,  the  burned  fuel  volume  fraction  of  stoichiometric 

P.  the  nondtmensionul  form  of  reaction  pressure  (/>)  defined  as  P  =/»■  I  atm,  and 
7",  the  nondimensional  form  of  initial  mixture  temperature  ( 7*, )  defined  as 
r-Ti/ 298°  K. 

The  correlation  coefficients  in  the  above  equations  were  determined  for  the 
following  cases  to  offer  more  convenience  in  use:  (I)  for  individual  fuel  families 
(0.5  <,4<  1.01;  (2)  for  hydrogen  with  fuel  air  ratio  ranges.  0.|4<iA<0.6  and 
0.6  <  6  <  1 .0,  respectively ;  (3)  for  acetylene  with  fuel  air  ratio  ranges.  0.3I<A<0.6 
and  0.6  <4<  1.0,  respectively,  and  (4)  for  individual  fuels  with  respective  fuel/air. 
ratio  ranging  from  lean  flammability  to  stoichiometric.  For  the  sake  of  convenience 
one  could  use  one  of  the  coefficient  matrices  for  individual  fuel  families  in  Table  I 
instead  of  one  for  individual  fuels  in  Tabic  II,  or  vice  versa.  Note  that  if  the  coef¬ 
ficient  matrix  in  Table  II  is  used,  Eq.  (6),  winch  implements  the  variation  in  the 
number  of  carbon  atoms  in  the  fuel,  is  not  required  in  the  computation,  making  it 
simpler.  In  addition,  the  coefficients  for  lean  mixtures  of  hydrogen  and  acetylene 
fuels  are  zero  since  there  is  almost  no  effect  of  the  reaction  pressure  on  Ta  for  the 
cases. 

Use  of  New  Formulas  for  Ta  Estimate 

Complex  as  it  may  look  at  first  glance,  the  use  of  the  coefficient  matrix  for  com¬ 
puting  r„  is  relatively  simple  on  a  pocket  calculator  as  well  as  in  computer  coding 
as  a  part  of  combustion  modeling.  In  order  to  familiarize  oneself  with  its  use.  the 
typical  combinations  of  fuel  jir  ratio  and  reaction  condition  are  summarized  in 
Table  111  by  listing  the  coefficients  needed  for  each  case 

I'  is  expected  that  the  use  of  the  present  formulas  for  the  7\  estimate  will  help 
facilitate  various  analyses  of  air  breathing  combustion  systems.  It  will  be  useful  in 
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a  wide  range  of  applications  since  the  formula  is  expressed  in  terms  of  the  major 
variables  determining  T„\  the  fuel/air  ratio,  the  reaction  pressure,  the  initial  mixture 
temperature,  and  the  number  of  carbon  atoms  in  individual  fuels.  The  accuracy  of 
the  formula  may  be  as  great  as  one  percent  at  either  of  the  lower  or  upper  limits  of 
fuel/air  ratio  considered  for  the  formula  but  much  smaller  otherwise. 


TABLE  m 

Coefficient  matrix  needed  in  individual  cases 


Mixture  Use  of 

and  matrix 

reaction  cond. 

Coefficients  of  fuel 
family  in  (IMIV) 

Coefficients  of  individual 
fuels  in  (VMVII) 
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A  Parametric  Analysis  of  Radiation  Heat  Transfer  in 
Direct  Injection  Diesel  Combustion 


S.L.  Chang,  X.L.  Yang  and  K.T.  Rhee 
Department  of  Mechanical  and  Aerospace  Engineering 
Rutgers.  The  State  University  of  New  Jersey 
New  Brunswick,  New  Jersey  08903 


ABSTRACT 

Radiation  heat  transfer  In  a  direct 
inject  Ion-type  diesel  engine  Is  investigated  by  a 
computational  parametric  analysis  for  various 
engine  and  combustion  variables.  A  new  model  has 
been  developed  to  compute  the  spectral  radiation 
heat  flux  incident  on  various  locations  In  the 
combustion  chamber  and  to  Incorporate  it  with  the 
spectral  emlsslvity  of  the  chamber  wall  surface. 
The  model  mainly  differs  from  previous  ones  in 
several  ways:  (1)  It  uses  the  in-cylinder  species 
distribution;  (2)  a  new  coordinate  transformation 
method  is  introduced,  instead  of  using  geometric 
factors  or  the  zonal  method;  (3)  it  computes  the 
spectral  volume  absorptance  of  the  combustion 
products  through  optical  paths;  (4)  a  new 
formulation  and  Integration  method  are  employed 
for  the  governing  equation  of  radiation  heat  flux; 
etc.  The  main  purpose  of  the  present  paper  Is  to 
report  some  of  results  from  a  parametric  analysis 
carried  out  by  using  the  new  computational  model. 


A  NEW  MODEL 

A  better  modeling  of  radiation  heat  transfer 
is  desired  for  a  more  comprehensive  analysis  of 
thermal  loading  of  a  diesel  engine.  A  rigorous 
way  of  computing  the  radiation  heat  flux  Incident 
on  a  particular  location  of  the  combustion  chamber 
wall  will  be  implementing  the  equation  of 
radiation  heat  transfer  as  accurately  as  possible. 
The  equation  calls  for  the  detailed  distributions 
of  optical  and  thermal  properties  along  individual 
optical  paths  in  the  hemispherical  volume  faced  by 
the  location.  Several  difficult  issues  remain 
to  be  overcome  in  achieving  the  above  goals  for  a 
long  time.  They  are  the  proper  use  of  the 
In-cylinder  species  and  temperature  distributions 
In  the  modeling,  the  computation  of  optical  and 
thermal  properties  of  the  species  along  the 
Individual  optical  paths,  the  accurate 
Implementation  of  the  equation  of  monochromatic 
radiation  heat  transfer,  to  name  a  few.  Each 
problem  is  discussed  elsewhere  by  proposing  the 
authors'  solution.  However,  a  brief  description 
of  the  authors'  methods  for  the  Individual 
problems  are  made  in  the  following.  Rote  that  the 
main  purpose  of  the  present  paper  is  to  report 
results  from  a  parametric  analysis  made  in  the 
present  numerical  method,  which  are  not  included 
In  others. 


In  order  to  properly  evaluate  the  space- 
resolved  nature  of  the  radiation  heat  transfer  in 
diesel  combustion.  It  Is  important  to  determine 
the  Instantaneous  in-cylinder  species  and 
temperature  distributions.  It  cannot  be 
overemphasl zed  that  an  assumption  of  uniform 
species  distribution  in  diesel  combustion  is 
incapable  of  properly  evaluating  the  spatial 
dependence  of  heat  transfer.  Such  space-resolved 
distribution  data,  which  may  be  found  by  using 
either  experimental  methods  [1-3]  or  computational 
means,  presently  are  not  sufficient  for  an 
extensive  analysis  of  the  radiation  process,  in 
particular,  for  a  parametric  study  of  the  process. 
As  an  alternative,  a  plume  equation  which  seems  to 
reasonably  express  the  instantaneous  species 
distribution  has  been  Introduced  for  radiation 
heat  transfer  analysis  as  [4], 

f  -  f 0  exp(-ao-b*2-cz2)  (1) 

where  f  represents  the  local  burned  fuel/air  ratio 
or  COj  or  HgO  concentrations,  etc.,  and  fg,  a,  b, 
C  are  the  ones  to  be  respectively  determined  at 
each  engine  crank  angle  by  using  the  above 
mentioned  methods.  For  simplicity  but  perhaps  as 
one  of  the  only  presently  available  alternatives, 
the  burned  fuel  air  ratio  may  be  used  for 
computing  the  Instantaneous  local  temperature 
distribution  [5],  The  above  plume  equation  enables 
one  to  choose  plume  geometries  of  varied  species 
distributions.  Therefore,  the  above  equation  can 
be  used  In  a  parametric  analysis  of  radiation  heat 
transfer  In  a  combustion  chamber  when  one  chooses 
plume  details  as  independent  variables,  i.e.,  by 
varying  values  of  the  constants  in  the  equation. 

Even  though  the  in-cylinder  species 
distribution  is  given,  very  likely  in  a  coordinate 
system  with  respect  to  the  injection  nozzle,  it  is 
needed  to  find  the  species  distribution  along  the 
optical  paths  centered  at  the  location  of  which 
incident  heat  flux  is  sought.  Because  of  its 
complexity,  simplified  methods  have  been  employed 
in  previous  radiation  computations;  uses  of 
geometric  factors  [5,6],  zonal  method  [7],  etc.  A 
new  coordinate  transformation  method  is  Introduced 
for  finding  an  accurate  species  distribution  along 
each  optical  path  centered  at  individual  locations 
over  the  combustion  chamber,  when  the  in-cylinder 
distribution  Is  given  with  respect  to  the  nozzle 
hole  [8].  A  short  summary  of  the  method  is  given 
here:  Referring  to  Fig.  1  which  shows  plumes 


distributed  In  a  cylindrical  coordinate  centered 
at  the  nozzle  hole,  R,  one  wants  to  find  the 
species  distribution  along  an  optical  path  faced 
by  a  particular  location  on  the  wall  (called  a 
detector  for  simplicity's  sake),  0. 

A  new  equation  was  derived  for  this  purpose 
as  shown  in  the  following  (Fig.  1)  or  for  a  given 
property  distribution  with  respect  to  the  nozzle 
as  shown  in  Eq.  (1),  the  property  distribution 
along  an  optical  path  (r,a,c),  s,  is  expressed  as. 


Fig.  1.  A  Schematics  of  Flame  Plumes 
in  a  Combustion  Chamber 
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ci  *  [ao2+b«2]sin2e/od  ♦  c  cos2  and  oi  and 
(1  *  0,1,  and  2)  are  second  order 
polynomials  expressed  by  using  Taylor's 
expansion  in  terms  of  x  *  rsin  s/o^,  r, 
a,  and  the  plume  angle  (refer  to 

Fig.  1).  " 


Upon  the  determination  of  the  species 
distribution  along  the  individual  optical  paths, 
the  'radiation  heat  flux  incident  on  the  detector 
through  each  path  is  computed.  In  order  to  achieve 
this,  the  following  two  steps  have  to  be  taken, 
i.e.,  to  compute  the  thermal  and  optical 
properties  of  the  species  and  to  implement 
integration  of  the  equation  of  radiation  by 
combining  all  the  above  results.  The  computations 
of  optical  properties  along  each  path  were  made  by 


finding  radiation  from  both  soot  and  gases 
respectively:  Raylei gh- 1 imi t  expression  and  the 
dispersion  equation  are  used  for  finding  the 
monochromatic  volume  absorptance  of  soots  [5,9] 
and  the  non-gray  semiempirical  band  model  proposed 
by  Edwards  and  Balakrishnan  [5,10]  is  employed  for 
computing  gas  volume  absorption.  The  summation  of 
radiations  from  both  contributions,  £t.\«  are  made 
equivalent  to  *  1- ( 1  -cg  x ) ( 1  -. c*s . a  ) »  where 
Eg.x  and  Es.x  ar®  the  spectra)  emissivity  due  to 
the  presence  of  gaseous  species  and  soot, 
respectively.  Notice  that  the  monochromatic 
property  variation  is  a  function  of  the  distance 
along  each  optical  path.  Regarding  the  equation 
of  radiation  heat  transfer,  the  following  equation 
was  derived.  In  reference  to  Fig.  1,  when  the 
distribution  of  the  total  volume  absorptance, 
<v(r),  is  found,  the  governing  equation  of  local 
directional  spectral  radiation  intensity,  I,(r), 
is  written  as 
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: r0 
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ebx(r')exp(-j 
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The  heat  flux  along  an  optical  path  of  a  solid 
angle  (e,c)  In  the  thermal  radiation  range  may  be, 
then,  found  from 
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(3) 


where  <\  is  the  spectral  volume  absorptance  at 
(r,9,  c) .  Combining  Eq.  (2)  with  the  above 
equation,  a  better  implementation  of  its 
computation  has  been  made  by  using  Gaussian 
integration  method  [11]. 


Due  to  the  repeated  computation  of  the 
adiabatic  flame  temperature  based  on  the  local 
burned  fuel/air  ratio  along  the  individual  optical 
paths,  a  convenient  way  was  sought  to  find  the 
temperature  values.  A  simple  form  of  equation  was 
developed  in  a  direction  of  its  proper  coupling 
with  Eq.  (3)  for  facilitating  the  subsequent 
solution.  The  new  equation  of  the  adiabatic  flame 
temperature,  Ta,  is  functionally  related  to 
several  variables,  i.e.,  the  fuel  volume  fraction 
of  stoichiometric,  the  reaction  pressure  p,  the 
initial  temperature  of  mixture,  T,  and  the  number 
of  carbon  atoms  in  the  fuel,  I,  as  shown  in  the 
following  [12]: 

T^  •  Aj[l+  A?inA  ♦  A3(tnA)2]  0.5  <  A  <  1.0 

(4) 

where  Aj,  A2  and  A3  are  functionally  expressed  in 
terms  of  p,  T  and  I.  The  constants  in  those 
expressions  for  most  of  practially  used  fuel/air 
systems  are  given  in  the  above  reference. 

In  addition,  the  spectral  surface  reflectivity 
is  incorporated  with  the  spectral  radiation  heat 
flux  incident  on  each  detector  over  the  combustion 
chamber  wall.  By  defining  the  extinction 
coefficient,  t,  along  the  optical  path  and  giving 
the  spectral  surface  reflectivities,  and  =2 
(see  Fig.  1),  the  radiosity  of  the  surface  1,  Bj, 
is  expressed  as 
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where  I(  *  e^T*  O  *  1*2)  and  ax  ’  1*£\. 


(5) 


Location  of  Detector  (  0(j,  jg.zg).  4 
computat i on  was  carr i ed  out  for  heat  fluxes 
through  the  locations  on  the  cylinder  head  across 
a  plume  at  =  0.5  and  along  tne  direction  of  a 
plume  axis  ($g  *  0)  and  their  results  are  shown  in 
Figs.  2  and  3  respectively.  The  trend  shows  the 
radiation  flux  quite  similarly  follows  the  species 


PARAMETRIC  ANALYSIS  OF  0IESEL  RADIATION  HEAT 
TRANSFER 

Since  the  details  of  instantaneous  fn- 
cylinder  events  are  not  known  due  to  either  the 
insufficient  experimental  data  or  the  limitations 
In  combustion  model  ing,  several  assumptions  are 
taken  in  the  present  computational  analysis.  The 
main  assumptions  are:  (1)  the  system  for  each 
computed  entity  is  Isolated;  (2)  the  species  in 
the  system  attains  an  equilibrium  at  successive 
moments;  (3)  the  temperature  distribution  is  that 
of  the  adiabatic  flame  temperature  determined  by 
the  burned  local  fuel/air  ratio  and  cylinder 
pressure;  (4)  the  fuel  is  injected  Into  the 
combustion  chamber  at  an  inclined  angle,  * 
150*;  (5)  the  scattering  in  the  radiation  is 
negligible;  and  (6)  the  standard  values  of 
constants  for  convenience  of  comparison  in  £q. 
(1)  are,  fo  for  soot,  H2O,  CO2  and  fuel/air  ratio, 
3  x  10*®,  0.01,  0.01  and  1.0,  respectively.  In 
addition,  the  distribution  constants  a,  b,  and  c 
are  0.6,  2.36  and  2  accordingly.  For  example,  the 
fuel/air  ratio  distribution,  may  be  written  as 
•  {1.0  exp  -[0 . 6  ( 0  /R  )  ♦  2.43(9/3r)  ♦  2.0(z/R)2]), 
where  R  is  the  piston  radius  and  9r  is  Vfthe 
number  of  plumes).  Further,  the  analysis  was  made 
for  the  following  engine  details;  the  surface 
temperature.  500*K;  number  of  spray  plumes,  4; 
piston  radius,  4.92  cm;  combustion  chamber  bowl 
radius,  3.44  cm  and  the  fuel  composition,  CigH34. 

The  computation  was  carried  out  to  evaluate 
the  following:  the  apparent  emissivity  (or 
normalized  radiation  heat  transfer),  ea  *  q/qr; 
the  spectral  emissivity  (or  normalized  spectral 
radiation  heat  transfer),  •  q  /qr  ;  the 

spectral  apparent  emissivity  (or  normalized 
relative  spectral  radiation  heat  transfer),  exa  * 
d  \/lr  x  .where  qr  «  sTjl.  Tr  •  2400*K:  \r  •  6.0  um; 
and  a  •  2 n®kVl5c2h3 ,  c  •  2.998x10^0  cm/sec,  h  • 
6.625x10*27  erg-sec,  and  k  •  1.380x10**®  erg/K. 
Additionally,  in  order  to  evaluate  the  effect  of 
the  presence  of  the  combustion  chamber  wall  on  the 
radiation,  three  different  computations  of  heat 
flux  were  evaluated  (as  shown  In  results):  Net, 
computation  of  radiation  on  heat  flux  through  the 
wall  with  inclusion  of  the  emission  of  the  spray 
plumes,  the  absorption  and  emission  of  the  wall 
and  the  reflection  from  the  surroundings;  No 
Reflection  ,  computation  Including  the  emission  of 
the  spray  plixnes  plus  the  absorption  and  emission 
of  the  wall  but  excluding  the  reflection  and 
emission  from  the  opposite  walls;  and  Spray  Plume, 
computation  of  radiation  streaming  out  of  the 
s&ray  plumes  only. 

The  results  obtained  by  using  the  above 
conditions  and  assumptions  are  shown  in  the 
following  for  various  parameters  with  relevant 
expl anatlons. 


Delector  Location  in  Anguior  Direction,  4>, 


Fig.  2.  Radiation  Heat  Flux  Through  a 
Cylinder  Wall  Along  the  Axial 
Direction  of  a  Plume 


Fig.  3.  Radiation  Heat  Flux  Through  a 
Cylinder  Wall  Along  the 
Circumferential  Direction 
Across  a  Plume 

distribution  in  the  corresponding  plume.  This 
suggests  the  greatest  portion  of  heat  flux 
Incident  on  a  detector  comes  from  the  radiation 
soiree  right  beneath  the  detector,  as  indicated  by 
the  directional  cosine  in  Eq.  (3).  A  similar 
conclusion  was  made  in  the  authors'  preliminary 
radiation  modeling  of  diesel  combustion  [5],  The 
local  minimum  around  og  *  0.1  may  be  of  surprise. 
In  order  to  find  its  reason,  the  net  directional 
radiation  intensity  was  calculated  at  different 
locations  along  the  axis  {*  g  •  0)  by  varying  the 
azimuthal  angle  ( ;)  of  integration  with  an  zineith 


!t>9 


angle  of  ^  *  39.7*,  as  shown  in  Fig.  4:  The  reason 
for  the  minimum  at  od  *  0.1  is  readily  explained 
form  the  figure  by  finding  that  while  the 
radiation  heat  flux  incident  on  the  detector  at 
the  nozzle  (  *  0)  is  from  all  of  the  four 
plumes,  the  heat  flux  on  the  detector  at  »  0.1 
is  from  the  species  from  one  plume,  i.e.,  the 
nearest  plume. 
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was  not  found  for  the  straight  end  of  the  plume 
(*(j  *  0)  but  for  »(j  »  10.  The  humps  in  the  curve 
were  found  to  occur  due  to  the  geometric 
configuration  of  other  plumes  from  a  similar 
analysis  to  results  in  fig.  4. 
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Location  of  Detector  on  Side  Wall  of  Piston  Bowl, 
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fig.  4.  Net  Directional  Intensity  Incident 
on  Zineith  Angle  of  39.7  deg  in  the 
Hemispheric  Volume  (r,a,  c)  faced  by 
Oetectors  Along  a  Plume  Axis 

The  radiation  heat  flux  Incident  on  the 
piston  surface  (with  the  cup  depth  of  0.74  cm)  for 
s  *  0  was  calculated  as  shown  in  fig.  5.  The 
trend  of  results  Is  predictable;  the  local  minimum 
at  a  d  «  0  is  caused  by  the  weak  emission  species 
right  underneath  the  nozzle;  the  abrupt  change 
around  34  *  0.7  is  due  to  the  presence  of  the 
piston  bowl  bank. 
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fig.  5.  Radiation  Heat  flux  Through  a  Piston 

Surface  Along  the  Axial  Oirection  of  a 
Plume 

A  similar  computation  was  conducted  for  the 
heat  flux  incident  on  the  side  wall  of  the  piston 
cup  as  shown  in  fig.  6.  The  strongest  heat  flux 


fig.  6.  Radiation  Heat  flux  Through  a 
Piston-Cup  Side  Wal 1 

The  heat  flux  on  the  cylinder  linear  (3^  =  1) 
may  be  of  interest:  Computation  was  made  for  zd  * 
-0.044  as  shown  in  fig.  7.  The  decrease  in  heat 
flux  with  »d  was  expected  because  the  species 
concentration  near  the  detector  becomes  leaner 
with  increase  in  The  relatively  flat  curve 
between  » d  •  10*  and  3^  *  35*  indicates  the 
increasing  effect  of  the  neighboring  plume  on  the 
heat  flux.  Note  that  the  heat  flux  striking  the 
cylinder  wal  1  at  •  1  is  stronger  than  that  at 
the  similar  location  on  the  cylinder  head.  Also 
note  that  the  heat  flux,  in  general,  is  much 
Stronger  on  the  cylinder  head  (fig.  3)  than  the 
cylinder  liner  (fig.  7). 
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fig.  7.  Radiation  Heat  flux  Through  a 
Cylinder  Liner  Along  the 
Circumferential  Direction 

Surf  ace  Emi ss 1 v 1 1  y .  In  figs.  2  and  3,  the 
energy-spray  plume  is  greatest  and  the  energy-net 
is  greater  than  the  energy-no  reflection.  This 
result  may  be  conveniently  explained  by  the 
following  example.  At  ;d  *  0,  since  the  detection 


surface  adsorbs  95X  of  the  incident  energy  (with 
reflection  of  5X)  and  emits  an  amount  of  energy 
(Tw  »  500*K)  equivalent  to  1.4X  of  radiation 
energy  from  the  spray  plume  (qsp),  the  energy-no 
reflection  is  93. 6X  of  the  energy-spray  plume. 
Computational  results  show  the  energy-net  is  96. 3X 
of  q$p;  the  difference  of  the  two,  3.2X 
(96.8X-93.6X)  of  qsp  is,  therefore,  the  energy 
reflected  from  the  opposite  surface  walls.  Since, 
from  this  example,  5X  of  reflection  and  1.4X  of 
emission  from  the  chamber  surfaces  result  in  a 
3.2X  of  increment  of  the  energy  detected  by  a 
location  of  the  chamber,  it  may  be  concluded  the 
following:  A  net  energy  equivalent  to  about  one 
half  of  the  energy  leaving  a  surface  will  arrive 
at  the  same  surface  when  incorporated  the 
emission  from  the  opposite  wall  with  the 
extinction  along  the  optical  paths. 

The  effect  of  the  above  trend  may  be  inferred 
for  different  surface  emissivitias  from  the 
results  shown  in  Fig.  8.  The  decrease  of  heat 
flux-net  with  decrease  in  emissivity  is  caused  by 
the  increase  in  reflectivity.  At  the  present,  the 
spectral  emissivity  of  the  wall  is  being 
incorporatedinto  the  analysis  for  further 
investigation  as  a  continuing  activity  at  Rutgers. 
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Fig.  8.  Effect  of  Surface  Emissivity  on 
Radiation  Heat  Flux  Through  a 
Cylinder  Head 

Soot  Concentration.  The  radiation  heat  flux 
Increases  with  the  soot  concentration  in  the  plume 
as  shown  in  Fig.  9.  The  Increase  is  greater  in 
the  lower  concentration  range  than  that  high 
concentration.  The  reason  for  the  convergence  of 
the  emissivity  at  high  soot  concentration  is  due 
to  the  exponential  function  relationship  of  the 
emissivity  to  the  volume  absorptance.  Another 
loqk  at  Fig.  9  will  lead  one  to  find  the  local 
minimum  in  heat  flux  does  not  occur  for  low  soot 
concentrations.  The  result  is  more  evidence  to 
explain  the  unique  receiving  geometric 
configuration  at  the  injection  nozzle  for  some 
cases. 
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Fig.  9.  Soot  Concentration  Variation  Effect 
on  Radiation  Heat  Flux  Through  a 
Cylinder  Head 

Soot  and  Gas  Radiation,  As  expected  from  the 
authors'  previous  study  [Sj,  it  is  found  that  the 
soot  radiation  is  much  stronger  than  the  gas 
radiation  (Fig.  10).  The  direct  summation  of  soot 
radiation  and  gas  radiation  exceeds  the  computed 
combined  radiation.  While  its  reason  may  be 
explained  by  looking  into  the  method  of 
combination  [4,8],  it  can  also  be  seen  from  the 
fact  that  some  gas  emission  is  absorbed  by  soots 
and  that  soot  radiation  is  similarly  attenuated  by 
the  presence  of  the  gas  species. 
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Fig.  10.  Comparison  of  Radiation  Heat  Fluxes 
Due  to  Soot  and  Gaseous  Species 

In  order  to  investigate  predominant  wave 
bands  of  radiation,  a  wave  number-resolved 
computation  was  made  as  shown  in  Fig.  11.  There 
are  three  distinctively  strong  bands  of  gas 
radiations  at  v  *  1,  3  and  10  jm  (recall  that  the 
normalized  wave  number  *  '/6  j).  In  fact,  each 
band  consists  of  many  small  bands  of  COj  and  HjO 
in  the  plume.  Over  the  spectra  of  radiation,  one 
finds  that  soot  radiation  is  dominant  in  short 
wavelength,  ><  2  jm,  while  the  gas  radiation 


becomes  important  in  longer  wavelengths,  *  >  4  urn 
(Fig.  12). 


Normalized  Wavenumber,  6«/A 


Fig.  11.  Spectral  Radiation  Heat  Flux 

Oue  to  Soot  and  Gaseous  Species 
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Fig.  12.  Relative  Spectral  Radiation  Heat 
Flux  Oue  to  Soot  and  Gaseous 
Species 

Pressure .  The  effect  of  pressure  on 
radiation  heat  transfer  was  studied  assuming  the 
injection  processes  and  the  gas  motions  are  not 
affected.  Figure  13  shows  the  computational 
results  of  the  pressure  effect  on  spectrum- 
resolved  radiation  heat  flux  at  *  0.  The 
variation  of  the  spectral  radiation  with  pressure 
is  somewhat  noticeable  at  the  long  wavelength 
region  (  •  >  3  -m)  where  the  gas  radiation  is 
Important.  Since  the  flame  temperature  is 
relatively  insensitive  to  the  pressure  variation, 
the  results  in  Fig.  13  only  show  the  pressure 
broadening  effect  on  the  gas  radiation  in  diesel 
combustion. 
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Fig.  13.  Pressure  Broadening  Effect  on 
Diesel  Radiation  Heat  Transfer 

Surface  Temperature.  The  weak  effect  of  the 
chamber  surface  temperature  variation  on  the 
radiation  heat  transfer  is  shown  in  Fig.  14.  There 
is  some  reduction  in  the  energy-net  due  to  the 
stronger  surface  emission  at  higher  surface 
temperature.  Such  an  effect  will  be  greatly 
increased  in  some  local  places  where  the  surface 
temperature,  however,  is  extremely  high  in  the 
engine  combustion  chamber. 


Fig.  14.  Surface  Temperature  Effect  on 
Diesel  Radiation  Heat  Transfer 


8urned  Fuel/Air  Ratio  vs.  Soot  Concentrat ion 
vs.  Flame  Temperature.  In  order  to  mutual ly 
compare  parameters  greatly  affecting  the  radiation 
heat  transfer,  Figs.  14,  15  and  16  are  offered. 
Several  important  points  found  from  the  results 
may  be  in  order:  the  effect  on  the  radiation  heat 
transfer  is  much  greater  with  the  relative 
variation  of  the  fuel/air  ratio  than  that  of  the 
soot  concentration.  The  greater  effect  of  the 
fuel/air  variation  is  caused  by  the  correspond i ng 
changes  in  C0£  and  H2O  concentrations  and  the 
adiabatic  flame  temperature.  Note  that  the 
temperature  effect  on  volume  absorptance  is  the 
quadruple  power  dependency  of  temperature  and  the 
soot  concentration  effect  on  it  is  the  linear 


relationship  of  soot  concentration. 
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fig.  15.  Comparison  of  Effects  on  Diesel 
Radiation  Heat  Transfer  Due  to 
Variations  in  the  Soot 
Concentration  and  the  fuel/Air 
Ratio 


fig.  16.  Comparison  of  Effects  on  Relative 
Spectral  Radiation  Heat  Transfer 
Due  to  Variations  in  the  Soot 
Concentration  and  the  fuel/Air  Ratio 

Since  the  adiaPatic  flame  temperature  is  not 
realistic,  at  all,  in  diesel  combustion,  the 
sensitivity  of  radiation  heat  transfer  to  the 
temperature  variation  is  worth  investigating.  Note 
that  the  actual  flame  temperature  in  diesel 
combustion  may  either  higher  or  lower  than  the 
adiabatic  flame  temperature  due  to  the  post¬ 
reaction  compression,  heat  transfer,  etc.  figure 
15  shows  the  effect  of  temperature  on  the 

radiation  heat  flux.  Note  that  there  Is  an 

assumption  of  no  change  in  gaseous  species 

concentrations  In  this  variation,  unlike  in  the 
fuel/air  ratio  variation.  As  mentioned  above, 
because  of  the  quadruple  dependence  of  temperature 
on  the  volume  absorptance  and  due  to  the  weak 
effect  of  the  variation  in  gaseous  species 

concentrations  on  the  radiation,  the  impact  of 


temperature  variation  on  the  process  is  concluded 
to  be  very  high  as  shown  in  the  results,  figure 
16  shows  the  stronger  effect  of  the  temperature 
variation  on  the  radiation  heat  transfer  at 
shorter  wavelengths.  This  may  be  well  explained  by 
Plank's  radiation  function. 


SUMMARY 

A  new  radiation  heat  transfer  model  of  direct 
injection  diesel  combustion  has  been  developed. 
The  model  has  been  used  for  a  parametric  analysis 
of  radiation  transmission  in  a  diesel  engine.  Some 
of  the  results  may  be  misleading  due  to  the  lack 
of  sufficient  combustion  details  of  the  engine. 
Nevertheless,  the  present  results  clearly  show 
that  the  engine  radiation  processes  may  be  well 
analyzed  by  using  the  new  model  and  that  some  new 
findings  are  worth  noting: 

(1)  The  most  portion  of  radiation  heat  flux 
incident  on  a  location  in  the  chamber  wall  comes 
from  the  radiation  source  right  near  the  location. 

(2)  The  effect  of  surface  emissivity  on 
diesel  radiation  may  not  be  negligible,  in 
particular,  when  the  surface  temperature  is  high, 
as  in  new  uncooled  diesel  engines. 

(3)  The  diesel  radiation  heat  transfer  is 
highly  dependent  on  the  soot  concentration  and 
weakly  sensitive  to  the  gaseous  species 
concentrations. 

(4)  The  pressure  effect  on  the  radiation  is 
almost  negligible  in  diesel  combustion. 

(5)  The  surface  temperature  of  the  chamber 
wall  is  not  a  sensitive  parameter  on  the  diesel 
radiation  heat  transfer,  but  it  will  not  be 
negligible  where  the  local  temperature  is 
exceedingly  high. 

(6)  The  flame  temperature  or  the  fuel/air 
ratio  variation  produces  very  high  impact  on  the 
radiation  process,  more  than  any  other  parameter 
studied. 
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NOMENCLATURE 

a  species  distribution  constant 

A  constants  in  adiabatic  flame  temperature 

equation 

b  species  distribution  constant 

3  radiosity  of  surface 

c  speed  of  light 

f  species  distribution  with  respect  to 

injection  nozzle 

fo  species  distribution  constant 

fy  soot  volume  fraction 

h  Planck's  constant 

I  radiation  intensity  or  number  of  carbon 

atoms  in  a  fuel 


k  8oltzmann  constant 

p  pressure 

q  radiation  heat  flux 

r  optical  path 

R  radius  of  cyl inder 

s  species  distribution  with  respect  to 

detector 

T  temperature 

components  in  cylindrical  coordinate 
r.s,;  components  in  spherical  coordinate 

•  surface  emissivity 

<  volume  absorptance 

wavelength 

•  fuel  volune  fraction  of  stoichiometric 

extinction  coefficient 
Stefan-Soltzman  constant 

Subscripts 

a  adiabatic 

b  black  body 

d  detector 

9  gas 

i  index  number 

op  optical  path 

r  reference 

s  soot 

t  total 

spectral 
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